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Abstract 



Using the gauge-gravity duality we study strongly coupled dynamics of fundamental 
flavours in large Nc gauge theories in a constant external field. We primarily focus on 
the effects of an external magnetic field. We use two holographic models realized in 
the Type IIB and Type HA supergravity and present a comparative case study. In both 
these models, by studying the dynamics of probe branes, we explicitly demonstrate and 
discuss the magnetically induced chiral symmetry breaking effect ("magnetic cataly- 
sis") in the flavour sector. We also study the associated thermodynamics and the meson 
spectrum and realize e.g. Zeeman splitting, stability enhancement of the mesons in the 
presence of an external magnetic field etc. By studying the quasinormal modes of the 
probe brane fluctuation in the hydrodynamic limit we also obtain an analytic disper- 
sion relation in the presence of a magnetic field in the Type IIA model. This dispersion 
relation consists of a propagating sound mode in the otherwise diffusive channel and 
is sourced by the quantum anomaly of the global U(l) current. We briefly discuss the 
effects of an external electric field and observe that the flavour bound states dissociate 
for sufficiently high electric fields and an electric current is induced. 
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Chapter 1: Introduction 



1.1 Motivation 

It is well-known that elementary particles interact by four fundamental interactions: 
gravitation, electromagnetism, the weak interaction and the strong interaction. The 
framework of quantum field theory has been extremely successful in describing many 
aspects of the last three interactions within remarkable experimental precision. There are 
however numerous aspects which remain poorly understood. One such long-standing 
puzzle comes from the strong interaction. Conventional quantum field theory meth- 
ods are mainly based on (but not limited to) perturbative techniques and are therefore 
generally inadequate to address various aspects of strong coupling dynamics. 

The quantum field theory describing the strong interaction is called Quantum Chro- 
modynamics (QCD). The elementary degrees of freedom of QCD are massive spin- 1/2 
particles called "quarks" and massless spin-1 particles called "gluons". There are sev- 
eral different quarks (named "up", "down", "charm", "strange", "top" and "bottom"). 
However depending on their masses the quarks can be broadly divided into two cate- 
gories: the light quarks (including the up, down and the strange) and the heavy quarks 
(including the charm, bottom and the top). As a first approximation we can set the 
masses of the light quarks to zero and the masses of the heavy quarks to infinity. In this 
approximation the only parameter in QCD is the dimensionless coupling constant. 

This coupling constant is a function of the energy scale, which is encoded in the 
QCD beta function. The QCD beta function was calculated a long time ago in refs. 
||62l 1 13211 and it was shown that it is an asymptotically free theory. This means that the 
coupling constant becomes small as the energy scale increases (in the UV) and becomes 
large when the energy scale becomes small (in the IR). This notion of UV or IR energy 
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scale in QCD is decided in terms of the only relevant scale Aqcd generated by dimen- 
sional transmutation. Physically, Aqcd (around 200 million electron volt (MeV)) cor- 
responds to the energy scale where the QCD coupling constant becomes of order one, 
hence the theory becomes non-perturbative. 

Thus perturbative field theory methods based on expansion in the small coupling 
constant are applicable only in the energy regime much larger than Aqcd- Indeed a lot 
of efforts have been made to address aspects of QCD within such perturbative regime, 
e.g. in ref. [I133II aspects of dense quark matter have been addressed in a perturbative 
framework, refs. [|8T1 11561 144| have perturbatively computed thermodynamic pressure 
of hot QCD. There is a vast literature on similar aspects. 

Although such perturbative approach has led us to remarkable insights about phases 
of QCD, this is only a small part of what we would want to understand. It is of vital the- 
oretical interest to build up useful and efficient non-perturbative techniques to address 
strong coupling dynamics in general. QCD is one such strongly coupled theory which 
posses (at least) two remarkable features in the low energy regime: confinement and 
chiral symmetry breaking. Both these effects have been observed in experiments, how- 
ever their dynamical origin is not well- understood. In a confining theory, the elementary 
degrees of freedom exist only in the form of bound states. In QCD, we can only observe 
the quark bound states such as mesons or baryons (collectively known as the hadrons) 
in the low energy regime, but not free quarks. Chiral symmetry breaking, on the other 
hand, occurs via the formation of a condensate consisting of fermions of opposite chi- 
ralities. The understanding of both these features remains an open problem till today. 

There exist powerful non-perturbative techniques. Lattice QCD is one such exam- 
ple. In this framework, spacetime is represented by discrete lattice points and quantum 
fields are represented by degrees of freedom living on each lattice site. Such a discrete 
formulation of spacetime provides a natural cut-off for the energy scale given by the 
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inverse lattice spacing. The QCD path integral, and therefore any observable, can in 
principle be computed by numerical simulations on a computer. 

Although lattice QCD is a very useful computational tool, it may not shed any light 
on the conceptual understanding. Moreover, the lattice computations have technical lim- 
itations of their own. Specifically, this framework is not (yet) equipped well-enough to 
address non-equilibrium time dependent dynamics. These computations are also limited 
to vanishing chemical potential due to the so called "sign problem". However various 
methods have been developed over the past decade to circumvent the sign problem; we 
refer to ref . [40] for a recent review. Moreover, lattice QCD remains a very useful tool in 
studying thermal properties of QCD at vanishing chemical potential; for a recent review 
see e.g. ref. PTTl . Some pedagogical introduction on lattice QCD can be found in e.g. 
refs. [l35l[TT4l[39ll. 

The effort towards understanding the strong coupling dynamics of QCD is not purely 
theoretically motivated. In recent years, the Relativistic Heavy Ion Collision (RHIC) 
experiments at Brookhaven National Laboratory have probed the strong coupling regime 
of QCD and thus pose a challenge for a theoretical understanding of the data. The 
quark-gluon plasma (QGP) produced at RHIC is at the energy scale roughly same as 
the Aqcd, thus it is expected that the dynamics is presumably strongly coupled. This 
intuition has further been supported by the fact that RHIC data is best described or fitted 
in the framework of ideal fluid dvnamics [|143l[T44l . 

Fluid dynamic^can be viewed as an effective description of a system when the mean 
free path of the constituent microscopic degrees of freedom is much smaller than any 
characteristic length scale of the system. Thus, by definition it implies the presence of 
strong coupling in the microscopic description of the corresponding system. For the case 
of RHIC, some direct evidences of strong coupling further exist and can be measured in 

'The standard pedagogical textbook for this is e.g. ref. IIIOII . 
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e.g., the large elliptic flow, jet quenching etc. For a pedagogical introduction of heavy 
ion physics and discussions on related phenomena see e.g. ref. [TTOII . The tools of lattice 
QCD are not very helpful to understand the RHIC results since many of the observed 
phenomena are time dependent non-equilibrium dynamics. 

Surprisingly string theory provides an unique toolkit to study strongly coupled 
aspects of QCD, at least in some qualitative sense. This is rooted on the idea of hologra- 
phy or gauge-gravity duality. At present, the holographic principle is a conjecture. The 
progress in understanding of black hole entropy and its dependence on its area rather 
than its volume led to the key idea [|154[|151[|28l . The precise statement of holography 
is as follows: any theory of quantum gravity in a given background is entirely equivalent 
to a non-gravitational theory describing the degrees of freedom living at the boundary of 
the spacetime and vice versa. This is also referred as the gauge-gravity duality. Further- 
more in the known examples of this duality the non-gravitational theory turns out to be 
a certain class of conventional quantum field theories. String theory provides us with a 
natural and well-defined framework where most of such examples of gauge-gravity dual- 
ity can be understood and used. However, much of the effective and working knowledge 
of the gauge-gravity duality can be learnt with minimal involvement of string theory. We 
refer the reader to ref. [|112|| for a recent review from this point of view. We will take a 
more conventional route and briefly review how this duality emerges in string theory. 

In ref. [I153II 't Hooft made the crucial step in recognizing the rank of the non-Abelian 
gauge group as an expansion parameter for the non-Abelian gauge theory. Namely, he 
proposed to consider an SU(A'^c) Yang-Mills theory in the limit Nc ^ oo keeping the 
't Hooft coupling A = gyy^Nc constant. In this so called large A^^^ limit many simplifi- 
cations happen. It can be shown that the partition function of the large Nc Yang-Mills 
theory is dominated by the planar diagrams only and the next corrections occur at the 
order of l/N^. Viewed as an approximation scheme for QCD (where Nc = 3), this 
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gives a result correct up to about 10% error margin and thus not completely irrelevant. 
Furthermore, the 1 /Nc expansion of the Yang-Mills partition function takes an identical 
form of the genus expansion of the string partition function, thus suggesting a possible 
connection between string theory and large N^. Yang-Mills theory. 

It was realized that string theory naturally contains various extended objects called 
"branes". A certain kind of such "branes" called the D-branes was originally realized in 
ref. [I125II : however now it is well-known that there are several kinds of branes in string 
theory. We refer to the book in ref. [,78. 1 for an extensive discussion on the properties of 
various kinds of branes. Here we will focus entirely on the D-branes. These extended 
objects are the places where an open string can end. Such an extended object with p 
spatial dimensions is called a Dp-brane. It was further realized that these Dp-branes arise 
as solutions of supergravity which is obtained as the low energy limit of string theory. 
Now the crucial step was to realize that the low energy dynamics of Nc coincident Dp- 
branes can be equivalently described by the gravitational background of the Dp-brane 
and an SU(A''c) Yang-Mills theory in {p + l)-dimensions living at the boundary of the 
Dp-brane background. 

String theory provides a large class of such examples. The existence of Dp-branes 
in both Type IIA and Type IIB supergravity (obtained as the low energy limits of Type 
IIA and Type IIB string theory respectively) provides a rich arena to explore and realize 
numerous examples of gauge-gravity duality. The most popular example of this is the 
duality between M = A Super Yang-Mills (SYM) theory in (3 + l)-dimensions and the 
gravitational background of anti de-Sitter (AdS) space in 5-dimensions. 

However none of the gauge theories mentioned above resemble QCD in precise 
details. It is still an open problem to construct the precise gravitational dual to QCD. The 
idea of using, e.g. N = A SYM to understand features of QCD is therefore completely 
based on the concept of universality. This idea was initially explored in refs. Il32l l33ll 
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to learn about the phase diagram of QCD using the dual gravitational background. See 
also refs. [l63l l36l l37l for discussions on related physics. For an early review exploring 
this idea we refer to the Chapters 1 8 and 19 of ref. [1781 . A more recent and less technical 
review can be found in ref. [l64ll . 

The notion of universality class is quite old in statistical mechanical systems. Phase 
transitions in physical systems are described by critical exponents which capture the 
behavior of physical quantities near the phase transition. Various physical systems 
have been observed to have the same critical exponents based on which many macro- 
scopic systems can be grouped into one set of universahty class. We refer the reader 
to ref. HI 6011 for a more detailed discussion of the physics of phase transitions and the 
notion of universality class. However, here we will consider a much relaxed notion of 
universality. We will merely concentrate on any statement (qualitative or quantitative) 
that holds universally for such large Nc gauge theories. 

That such universal statements exist is further supported by the following observa- 
tion. At zero temperature, the physics of QCD differs considerably from the physics of 
A/" = 4 SYM (or any such large Nc gauge theory for which we know the gravitational 
dual). Firstly, the degrees of freedom do not match. Moreover QCD is confining, with 
a discrete spectrum and running coupling constant whereas N = A SYM does not share 
any of these features. However at finite temperature both these descriptions consist of 
deconfined strongly coupled plasma of gluons and various matter fields. Furthermore as 
it has been emphasized before, the best model describing the dynamics of quark gluon 
plasma (QGP) at RHIC is fluid dynamics which is an effective long-wavelength descrip- 
tion and does not necessarily rely on the microscopic details. 

Interestingly RHIC provides us with a remarkable quantitative evidence of such pos- 
sible universality. Recent investigations, e.g. ref. [|148 l| seem to suggest that the ratio 
of viscosity to entropy density for the QGP plasma r^/s ~ h/{A7ikB), where h is the 
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Planck constant and A;^ is the Boltzmann constant. Although presently we do not have 
any QCD calculation for this ratio, using AdS/CFT refs. [fT29l fBOl |95l |96|| calculated 
the ratio to be r^/s = h/ {AixkB). It was further argued in refs. [|96ll29ll30l that the ratio 
rj/ s = h/{A'KkB) is universal for any gauge theory with a classical gravity dual. This 
remarkable qualitative and quantitative (in terms of orders of magnitude) resemblance 
to the observed data at RHIC is suggestive of an universality we alluded to. 

Let us briefly comment on a further interesting observation. Generically, the finite 
temperature strong coupling physics seems to have universal properties which can be 
obtained using the methods of gauge-gravity duality and ideal fluid dynamics. Thus a 
natural question arises: whether there is a connection between the dual gravitational and 
the fluid dynamic descriptions. During the past few years this possible connection has 
been analyzed in detail and it has been concluded that there is a precise map between a 
class of inhomogeneous, dynamical black hole solutions in asymptotically AdS space- 
time and arbitrary fluid flows in the strongly coupled field theory living at the boundary. 
This is now known as the fluid-gravity correspondence[26J. We refer to the ref. HI 3411 
for an extensive review and further references on this topic. 

Thus we can boldly step forward in using techniques of gauge-gravity duality to fur- 
ther understand the strong coupling dynamics of QCD in some universal sense. How- 
ever we will not attempt to compare any of our results with any experimental data, at 
most we will merely comment or compare with the results and expectations from QCD 
whereever necessary. Our approach should thus be viewed as an attempt to use the 
tool of the AdS/CFT correspondence (or gauge-gravity duality) to learn about strong 
coupling dynamics of large A^^^ non-abelian gauge theories in as much generality as pos- 
sible. 
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1.2 Gauge-gravity duality from string theory 

String theory is a vast subject and we will just refer the standard pedagogical text- 
books [|60l [6T1 11261 11271 1781 for interested readers. Here we will briefly review some 
basic facts about string theory. The basic degrees of freedom of string theory are one 
dimensional objects called strings. The fundamental parameter of the theory is the char- 
acteristic length scale is = which sets the string tension Tg = l/(27ra'). The 
interaction strength between the strings is given by the dimensionless string coupling 
constant denoted by gg. Strings can be both open and closed. Open strings are one- 
dimensional line segments with a boundary and two end points. On the other hand, 
closed strings are loops with no boundary or end points. 

Based on the presence of supersymmetry in the underlying action, string theory can 
be divided into two broad categories: the bosonic string theory and the superstring the- 
ory. There are various kinds of superstring theories, which go under the technical names 
of Type I, Type IIA, Type IIB and heterotic string theories etc. These apparently differ- 
ent theories are actually related by certain duality transformations and are connected to 
what is known as M-theory. 

Although presently we lack a complete understanding of string theory (or M-theory), 
the low energy limit of this theory is rather well-studied. In the low energy limit, string 
theory /M-theory reduces to ten/eleven dimensional supergravity. It is in the framework 
of supergravity where many aspects of the gauge-gravity duality have been discovered 
and used. Note however that many aspects (e.g. later when we consider branes and 
fundamental strings with gauge fields on them) of the duality depends on the full string 
theory and not just its low energy limit. Nonetheless, all our computations in the later 
chapters are performed within the supergravity approximation. In this thesis we will 
discuss two different (the Type IIB and the Type IIA) supergravity frameworks in par- 
ticular. 
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1.2.1 The Type IIB theory: low energy dynamics on D3-branes 

In Type IIB theory, we can have Dl, D3, D5 and D7 (stable) branes. Here we will restrict 
ourselves to the low energy limit of D3-branes only. 

To lead to the conjecture, let us start with a stack of Nc parallel D3-branes. In this 
set-up there are two kinds of excitations: the open strings that start and end on the D3- 
branes and the closed strings that propagate through the ten dimensional bulk. At low 
energies (i.e. energies much smaller than £j^), we can integrate out the massive modes 
to obtain an effective action given by -S" = -S'buik + -S'brane + -S'int- The bulk action S'buik 
is the action for Type IIB supergravity up to corrections in a' — (which appear as 
higher curvature corrections to the Type IIB supergravity action). On the other hand, 
the massless spectrum of open strings starting and ending on the D3-branes constitute 
the vector supermultiplet of J\f — A SYM in (3 + l)-dimensions with U(A^c) gauge 
group. The presence of the gauge group can be understood as follows. The massless 
open string spectrum ending on a stack of Nc parallel D3-branes is labeled by an index 
known as the Chan-Paton factor. In case of oriented strings, it can be shown that these 
indices transform in the adjoint representation of an U(A^c) gauge group. This gives the 
brane action 5'branc- There are also interactions between the bulk modes and the brane 
modes which is denoted by 5'int. In the limit a' 0, these interactions are suppressed 
reflecting that gravity becomes free at large distance scales. Thus in the low energy 
limit, we get a complete decoupling between the open and the closed string modes. Let 
us now discuss this decoupling in some details. 

1.2.1a The decoupling limit and the correspondence 

In supergravity D-branes are massive charged objects sourcing the supergravity fields. 
To illustrate the decoupling limit, let us start from the extremal D3 -brane background in 
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flat space obtained as a solution of Type IIB supergravity equations of motion. Follow- 
ing the notation of ref. iTTSl this background is given by 



ds^ = H'^^^ri^^dx^'dx'' + Hl'^dx'dx' , 
C(4) = H^^g^^dx^ A... Adx^ , 

e' = g,, H, = l + ^, = ATTg,N,a'\ (1) 

where n = 0, ... 3 and i = 4, ... 9 and ?7^j, is the flat Lorentz metric; r is the radial 
coordinate and is a function harmonic in the six directions transverse to the D3- 
branes. C4 is the four-form RR potential sourced by the D3-brane and Nc is an integer 
measuring the flux quanta of the RR field strength. This supergravity background has 
a horizon at r = and the boundary is located at r = 00. From the point of view 
of a boundary observer, the low energy Type IIB string modes near the horizon will 
be infinitely redshifted by the factor H^^^^ and thus these low energy excitations near 
the horizon should describe the low energy massless degrees of freedom seen by the 
asymptotic observer. On the other hand, the low energy massless excitations of Type 
IIB string theory is described by the Type IIB supergravity as seen by the boundary 
observer. 

Now any bulk excitation in the full background interacts with the near-horizon exci- 
tations, which is captured by the absorption cross-section of graviton scattered by the 
near-horizon regime, a ~ co^R^ as obtained in refs. [|67ll90ll , where tu denotes the energy 
of the graviton. However, in the low energy limit (denoted by w — t- 0), cr — )■ since the 
bulk excitations have a wavelength much larger than the gravitational size of the brane. 
Thus it is appropriate to conclude that bulk excitations decouple from the near-horizon 
excitations in the low energy limit. 
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Now we will precisely define what we mean by the "near-horizon" limit of the back- 
ground in ([T|. We achieve the near-horizon limit of ([T]) by sending a' — t- and keeping 
u = r joL fixed. In leading order in ol we get the following background !! 10511 

C4 = jdx^ A . . . Adx^ , 

gsU^ 

6"= 9s. (2) 

The metric in ^ is the metric of AdSs x background. Thus we learn that the near- 
horizon geometry of D3-branes is identical to an anti-de Sitter space. 

Therefore we arrive at two physical descriptions for the low energy physics of a 
stack of Nc D3-branes: the A/" = 4 SYM theory and the Type IIB string theory in 
AdSs X background. The essence of Maldacena's AdS/CFT correspondence is to 
conjecture that these two descriptions are entirely equivalent. Let us now formally state 
the correspondence. 

Type IIB string theory on AdS^ x background is equivalent to Af = 4 SYM theory 
in (3 + 1)- dimensions. 

One hint in favor of this correspondence comes from the non-trivial matching of the 
global symmetry group in the gauge theory with the isometry group of the geometry. Let 
us therefore discuss this in some more details. We will not explicitly prove or demon- 
strate any of the following assertions, for our purposes it would suffice to just review 
them. 

A/" = 4 SYM is a maximally supersymmetric gauge theory with SU(A^c) gauge 
group. The degrees of freedom are contained in the gauge multiplet of this theory which 
consists of gauge fields, four left-handed Weyl fermions and six real scalar fields. All 
these fields transform in the adjoint representation of the SU(Ai'c) gauge group. This 
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theory has only one gauge coupling constant usually denoted by gym- A remarkable 
feature of this theory is the exact conformal invariance implying that the coupling con- 
stant gYM does not run at all. The global symmetry group of this theory is the conformal 
group in (3 + 1) -dimensions, which is given by S0(4, 2). J\f = A SYM also has sixteen 
real supercharges which can be rotated into each other by an SU(4)r group called the 
R-symmetry group. Thus the global symmetry group is given by S0(4, 2) x SU(4)r. 
For future reference, let us also note that SU(4) is locally isomorphic to S0(6), i.e. 
SU(4) ~ S0(6). 

Now we focus on the geometric side. Anti-de Sitter spaces are defined as spaces 
with constant negative curvature. These spaces are obtained as a solution to Einstein's 
equations with a negative cosmological constant. The (4 + 1) -dimensional anti-de Sitter 
space can be represented as the following hyperboloid in MP 

4 

X',+X',-J2xf = R\ (3) 

i=l 

where R is the radius of curvature of the hyperboloid. It is evident from this description 
that AdSs space has S0(4, 2) isometry acting as the rotation group in (4+2)-dimensions. 
A simple way to parametrize the solution of eqn. ([3]) is to set 




and u > 0. These coordinates are called the Poincare patch which cover only half of the 
hyperboloid. There exists also a global patch which cover the entire hyperboloid, but 
we refrain from discussing it since it will not be relevant for out purposes. For interested 
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readers, we refer to the review in ref. [|3l. In the Poincare patch, the induced metric on 
the AdSs space is given by 




which is what we have met in eqn. ([2]). In eqn. ([2]) we obtained a background in a product 
form AdSs x S^. The isometry group of is S0(6) corresponding to the rotation 
group. Thus the geometric background has an isometry group SO (4, 2) x SO (6), which 
is exactly identical to the global symmetry group of the J\f = A S YM theory. The relation 
between the gauge coupling and the string coupling is simply given by Qs = gy^- For 
further details we refer the review [|3l. 

We have, so far, stated the strong form of the correspondence. However, it is not 
known how to quantize string theory on the curved AdSs background. Thus the strong 
form of the conjecture is hard to confirm. A weaker form of the conjecture can be based 
on the 't Hooft limit with a fixed 't Hooft coupling defined by A = gyM^c, while taking 
the limit Nc oo. In this large Nc expansion the planar A/" = 4 SYM is conjectured to 
be equivalent to Type IIB string theory perturbatively expanded in powers of the string 
coupling gs = X/N^. A further weaker form of the conjecture is in the large 't Hooft 
coupling limit, i.e. A ^ 1. Using the definition of R from eqn. ([T]), we can conclude 
that in this limit R^/a''^ ^ 1, meaning that the curvature scale is much larger than the 
string scale, where classical supergravity is a valid description. Thus the weakest form 
of the conjecture states that classical Type IIB supergravity is dual to the = 4 SYM 
in the strong 't Hooft coupling limit. In this thesis we will use the weakest form of the 
conjecture and use it to learn strong coupling aspects in a purely classical (super) gravity 
set-up. 
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1.2.1b The dictionary of the correspondence 

Let us now review the precise map between the gauge invariant operators in A/" = 4 
SYM and the supergravity fields. To demonstrate this let us consider a simple massive 
scalar field (p propagating in the background AdSs space. The equation of motion for 
the scalar is given by the massive Klein-Gordon equation 

V^cj) - = , (6) 



which has two independent asymptotic solutions given by 

~ Vo + u-^{0^) , A = 2 + V4 + m2i?2 , (7) 

where A is the conformal dimension of the field. The first term with coefficient 0o 
is the non-normalizable solution and the second term with the coefficient (O^) is the 
normalizable solution. It was suggested in ref. [I158II that the non-normalizable mode 
should be identified with the source and the normalizable mode should be identified 
with the expectation value corresponding to the source. Furthermore it was suggested 
in refs. I1158[ [66l that the precise form of the dictionary is given by equating the gen- 
erating functional of the conformal field theory with the partition function of classical 
supergravity 

^exp {^j d^xcly^O^ ^ = Z^^^,^ [(/)o] = exp (i^sugra [<Po]) , (8) 

where 0o is the boundary value of the bulk scalar field. The proposal of refs. ||158[ |66ll 
holds true for any bulk supergravity field irrespective of whether it is a scalar, vector 
or tensor. This proposal can therefore be viewed as the operational definition of the 
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AdS/CFT correspondence in its weak form. Furthermore, it is noteworthy that this pro- 
posal can be generally extended to any AdSd+i/CFT^ correspondence. 

However the precise correspondence in ([8]) does not give a prescription of how to 
determine the dual bulk field corresponding to a given operator in the gauge theory and 
vice versa. Usually this identification is uniquely fixed based on the scaling dimension 
and the symmetry property of the operator or the dual bulk field. Another caveat with 
this precise statement is that on the face of it both the field theory generating functional 
and the on-shell supergravity action are divergent quantities. Thus we need to introduce 
a consistent regularization and renormalization scheme. The divergence in the super- 
gravity action comes from the infinite volume of AdS space. Thus while evaluating the 
supergravity action integral we introduce a cut-off at some u = const hypersurface and 
add appropriate counterterms on the hypersurface to cancel the divergences. This is for- 
mally known as the "holographic renormalization" scheme and we refer the review II145II 
for further details on this. 

It is interesting to note that the regularization scheme we just described is very sim- 
ilar to the familiar regularization scheme in field theory if we identify the radial direc- 
tion u of the AdS background to be the renormalization scale of the dual field theory. 
This identification is indeed true and can be understood simply as follows. If we scale 
Xfj, — )■ ax^, where a is some real positive number; then for the AdS-metric in ([5]) to be 
invariant under the scaling transformation we must require u a^^u, which is how 
the energy must transform. In the AdS background, m — > is the region deep inside 
and u — oo corresponds to the region near the boundary. Motivated from the above 
observation we will sometimes refer to these as "the IR" and "the UV" respectively. 
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1.2.1c The correspondence at finite temperature 

Let us review the key aspects of the correspondence at finite temperature. As we will see 
later, for the most of the cases we will be interested in finite temperature physics. The 
notion of a temperature can be introduced very easily in any quantum field theory. This 
is simply achieved by Euclideanizing the time-coordinate by an Wick rotation t — )• —it^ 
and compactifying this Euclidean direction on a circle with period (3 = 1/T, where T 
is the corresponding temperature. The frequency modes corresponding to this compact 
direction naturally becomes discrete, which are known as the "Matsubara frequencies". 
This is known as the "imaginary time formalism". It is to be noted that in this formalism, 
once temperature is introduced we loose the notion of time and therefore we cannot 
describe any time-dependent phenomenon. 

The dual gravity background of finite temperature J\f = A SYM is given by the 
AdSs-SchwarzschildxS^ background, where the Hawking temperature and the entropy 
of the black hole is identified with the temperature and entropy of the field theory ll65l 
I158L In view of the precise map in ([8]), the generating functional of the quantum field 
theory becomes the thermal partition function. Thus we get that the free energy of the 
field theory is given by the Euclidean on-shell supergravity action 

F = T^,„gra [M ■ (9) 

Therefore the thermodynamics oi M = A SYM can be calculated from the Euclidean 
gravity action. Finite temperature thermodynamic quantities such as the free energy or 
the entropy density, thus computed, scales as A'^^ as observed in ref. [|65l . This factor 
of Nl is understood as counting the degrees of freedom for the SU(Ac) gauge theory. 
Since all the degrees of freedom contribute to thermodynamic observables, we define 
this as the "deconfined" phase of the N = A SYM theory (as opposed to the "confined" 
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phase where degrees of freedom exist only in the form of bound states and thus provide a 
factor of in front of thermodynamic observables). All our analysis in the subsequent 
chapters will mostly be restricted to this "deconfined" phase only. 

As clear from the formahsm itself, this is not just limited to the finite temperature 
Af = 4 SYM. 

1.2.2 The Type IIB theory: introducing flavours 

As we have already reviewed, J\f = A SYM theory contains only adjoint matter sector 
whereas QCD contains both adjoint (the gluons) and fundamental matters (the quarks). 
We will focus in particular to understand the dynamics of these fundamental flavours 
and thus we need to introduce the flavour fields in the gravity set-up. 

So far we have focussed on the low energy dynamics of D3-branes only. In this 
case, open strings that start on the D3-branes must end on the D3-branes itself. These 
two end-points of the open string correspond to point charge in the fundamental and anti- 
fundamental representation of the SU(Ai'c) gauge group. Together they transform under 
the adjoint representation. It is therefore clear that in order to introduce fundamental 
flavour degrees of freedom we need to introduce open strings with only one end on the 
D3-brane. 

The simplest way in which this is achieved is by considering an additional stack of 
Nf D7-branes [[83l . The relative orientation of the Nc D3-branes and Nf D7-branes is 
given by the following array. Here " — " represents that the brane is extended along 
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Table 1: The brane construction including fundamental flavours. 

the corresponding direction whereas " • " represents that the brane is point-like in the 

17 



corresponding direction. It is to be noted that the D3-branes' directions are aligned with 
four of the D7-branes' directions. 

The low energy dynamics of the above configuration is now given by open strings 
that start and end on the D3-branes, open strings that start and end on the D7-branes 
and finally open strings that start on the D3/D7-branes and end on the D7/D3-branes. 
The low energy dynamics of the open strings that start and end on the D3-branes give 
the — A SYM as before. The D7-branes are extended in the non-compact direction 
transverse to the D3-branes and thus have infinite volume. Therefore the strings that 
start and end on the D7-branes yield a global symmetry as opposed to a gauge symmetry 
since the gauge coupling vanishes. 

The strings that are stretched between the D3 and the D7-branes give rise to the so 
called M = 2 hypermultiplet. The N — 2 hypermultiplet consists of two complex 
scalars and two Weyl spinors of opposite chiraUty. One scalar and Weyl fermion trans- 
form under the fundamental of SU( A^c) and the other pair transform under the conjugate 
representation. In an abuse of language, we will call that the hypermultiplet transforms 
under the fundamental representation. We will not discuss the explicit Lagrangian of 
this theory in details. 

A few words about the R-symmetry group are in order. As seen from the explicit ori- 
entation of the D7-branes the full S0(6)r symmetry is broken down to an S0(4)r sym- 
metry corresponding to the group rotation in the {4, 5, 6, 7}-directions and an S0(2)r 
symmetry corresponding to the group rotation in the {8, 9}-directions. 

This S0(2)r ~ U(1)r is identified as the chiral symmetry under which the two Weyl 
fermions are oppositely charged. Note that the D3 and the D7-branes can be separated 
in the {8, 9}-plane in which case the strings stretched between these branes will have 
a finite length and finite energy. This corresponds to giving mass to the hypermultiplet 
fields, which explicitly breaks the U(1)r symmetry. This is analogous to the chiral 
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symmetry breaking in QCD when the quarks become massive. Note that in principle 
we can give different mass to each flavour sector by separating each of the Nj different 
D7-branes from the D3-branes by different amounts in the {8, 9}-plane. However, we 
will content ourselves with the simple case when all the Nf D7-branes are coincident. 

In what follows we will restrict ourselves to the so called "probe limit" which means 
that we will take the large Nc limit keeping Nj fixed such that Nc ^ Nj. In this limit, 
the quantum corrections leading to the U(1)r anomaly are suppressed lll57L thus the 
U(1)r symmetry becomes an honest global symmetry. 

This probe limit has another interesting consequence. In the case when the mass 
of the hypermultiplet is zero and there is no background temperature, the theory is still 
conformal. It can be shown that introducing Nf flavours modifies the beta function of 
this theory by adding a positive contribution at the order Nj /Nc and is thus negligible 
in the probe limit [[88l . However in the parametric regime Nf ~ N^ the coupling of 
the theory increases with increasing energy and finally acquires a Landau pole where 
the gauge coupling constant becomes infinitely strong. This behaviour is completely 
different from QCD and requires an UV completion. 

The gravity dual picture is rather simple in this probe approximation. We will replace 
the near-horizon geometry of the D3-branes by the AdSs x background (or with the 
AdS-Schwarszchild background in the finite temperature case) and treat the D7-branes 
as probes in this gravity background which only "sense" the background geometry but 
do not backreact on it. This is similar to the "quenched approximation" used in the 
Lattice QCD studies to analyze the effect of a test quark in the background of gluons. In 
this probe limit, the action of the probe brane is given by the Dirac-Bom-Infeld (DBI) 
action. This action will be our main tool and we will discuss some general aspects of 
the DBI action in a later section. 
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1.2.3 The Type IIA theory: low energy dynamics on D4-branes 

Immediately after Maldacena's AdS/CFT correspondence, the low energy dynamics of 
a general Dp-brane was investigated in ref. [1741 to realize more examples of gauge- 
gravity duality. Here we will not review the details of Type IIA string theory or Type 
IIA supergravity. However, we will just recall that in Type IIA supergravity Dp-branes 
exist only with even values for p. In particular we will focus on the low energy dynamics 
of D4-branes. 

As clear from the discussions in the previous sections, the gauge theory represent- 
ing the low energy dynamics of a stack of D4-branes is (4 + 1) -dimensional. It 
was illustrated in ref. [1741 that the decoupling limit is now taken by sending the string 
scale a' —7- and the asymptotic string coupling Qs ^ oo while keeping the (4 + 1)- 
dimensional gauge coupling = (2txY ggy/o' fixed. The (4 + 1) -dimensional SU(Ai'c) 
gauge theory in this case is not a conformal theory. Moreover this is not a renormaliz- 
able theory and thus requires an UV completion at high energies. This UV completion 
is given by the (5 + 1) -dimensional N = (2,0) superconformal theory defined on the 
worldvolume of M5-branes compactified on a circle]^ However, we will not elaborate 
on this issue since it will not be relevant for our subsequent discussions and refer the 
interested reader to ref. [74] for further discussions on this. 

1.2.3a The decoupling limit 

As we have motivated ourselves to learn aspects of QCD (or at least (3 + 1) -dimensional 
gauge theory) using the gauge-gravity duality, we will restrict ourselves to the construc- 
tion described by Witten in ref. I1 159II . Let us discuss this set-up in some detail. Witten 
proposed to consider the near-horizon limit of a stack of N^. D4-branes compactified 

^M-branes are objects like D-branes in M-theory. For an extensive list of the examples of gauge- 
gravity duality arising from various D-branes and M-branes we refer to the classic review [31 • 
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on a spatial circle (an S^) with a spin structure on the that breaks supersymmetry 
(achieved by imposing an anti-periodic boundary condition for fermions in the adjoint 
of the SU(A^c) gauge theory). The near-horizon gravity background is given by 

3 



u 



J > ^(4) = -^e4, /H = !-( — ) , (10) 



where t and x\i = 1, 2, 3 are the four gauge theory directions and is the compact 
S^. The radial variable is denoted by u and Ukk is a constant. The S0(5) invariant line 
element is denoted by dQl, e(4) denotes the corresponding volume form and V4 is the 
volume of an unit 4-sphere. Also, $ is the dilaton and F(4) is the RR four-form field 
strength. Since is a compact direction, to avoid a conical singularity in the {x^, u} 
plane one should make periodic identification: 

Sx' = ^(^) =2vri?4. (11) 

This endows the background with a smooth cigar geometry in the {x'^, u} plane. When 
Ukk — 0, there is a naked singularity at = with no restriction on the period of the 
direction. 

The radius of curvature of the background is given by 

R'^TrgsN.el. (12) 

The gauge couphng of the five-dimensional field theory is given by 95 = {2-Kygs\/a' as 
alluded to above. The corresponding five-dimensional 't Hooft coupling is given by 

h^glK. (13) 
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Since the five-dimensional gauge coupling and the 't Hooft coupling are dimensionful 
quantities, we find a power law running of the dimensionless effective coupling given 
by [l74ll gig = g1 Ncu/ {2tx a'). There is a tower of Kaluza-Klein (KK) modes arising 
from the compactification along the direction, which sets the typical glueball mass 



scale as obtained by studying excitations around the background in ( 10): Mkk = 1/-R4. 
Below this scale, the effective low energy gauge coupling is given by gl = g\l (27ri?4) 
which defines a (3 + 1) -dimensional 't Hooft coupling A4 = g\ N^. 

Let us now comment on the regime of validity of supergravity approximation. In 
order for this to be true, we impose that the spacetime curvature must be much smaller 
than the string length scale. The maximum value of the Ricci scalar of the background 
in ( 10) goes as (UkkR^) hence we impose the condition 



>1 ^ glN,:^l ^ -^>1, (14) 



£2 ^^^-^^ R 



4 



which corresponds to the large (3 + 1) -dimensional 't Hooft coupling limit. 

Note that the dilaton is not constant in this background and string loop effects are 
not suppressed for arbitrary large values of u. However, the critical value of u at which 
the dilaton is order one can be obtained to be 

Mcutoflf ~ , p , (15) 
A5lt4 

which does not play any important role in the large Nc limit (with fixed A5). 

It was shown in ref. [| 15911 that the theory develops a mass gap. In the regime A5 <^ 
i?4, the scale of the mass gap is exponentially suppressed compared to I/-R4. Therefore 
the gauge theory is effectively (3 + 1) -dimensional below the KK-scale set by l/R^. 
However in the strongly coupled regime A5 ^ R^, the dual gauge theory develops a 
mass gap of the order I/-R4 and therefore cannot be separated from the KK scale. Thus 
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the gravity calculation (which is valid at the strong coupling regime) corresponds to a 
(4 + 1) -dimensional gauge theory which is non-renormalizable. The presence of the 
cut-off in the radial coordinate is a reflection of this fact. The scheme of holographic 
renormalization in this set-up is therefore somewhat subtle. The systematics of the 
holographic renormalization of such non-conformal systems (including the example of 
the near-horizon D4-brane geometry) has been worked out recently in ref. If82ll . 

A few words about the dual gauge theory are in order. Let us briefly mention the 
various degrees of freedom that constitute the low energy description of the D4-branes 
below the KK-scale. The massless degrees of freedom that arise from the open strings 
that start and end on the D4-branes are fermions, gauge fields and six real scalars trans- 
forming in the adjoint representation of the SU(Ai'c) gauge group. The component of 
the gauge field along the compact direction also appears to be an adjoint scalar as 
seen from the effective (3 + 1) -dimensional perspective. Since there is no supersymme- 
try, these scalars obtain mass via one-loop corrections and thus decouple from the low 
energy physics. 

The dictionary of this correspondence works in a similar manner as we discussed 
earlier. Any bulk supergravity field will have two types of modes near the boundary 
(which is now placed at some cut-off radial value of u): the non-normalizable and the 
normalizable, which are identified with the source and the vev respectively. The iden- 
tification of a bulk field with an appropriate boundary gauge theory operator is again 
based on its symmetry properties and engineering dimension. 

1.2.3b The finite temperature background 

We can discuss finite temperature aspect of the correspondence by Euclideanizing the 



time coordinate t — )■ —He in the gravity background ( 10) and compactifying along ^e- 



Now we have two directions which are compact. In this case, the x'^ circle smoothly 
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shrinks away at m = Ukk but the circle is fixed. One can also construct a finite 
temperature version by interchanging the role of the t and circles so that now time 
circle shrinks away at some value u = Ut but the circle remains fixed. It turns out that 
these are the only known Euclidean continuations of the background in equation ([10]) 
with the right asymptote. This second solution is given by (in Euclidean signature) 

ds' = (1)'^' [dxdx' + f{u)dtl + {dxf) + (1)"'^' (^-^ + u'dnl^ , 



where Ut is a new constant and T is the temperature. 

Therefore at finite temperature we have two candidate gravity backgrounds: the 



Euclidean version of (10) and the background in (16) and hence we have to resort to 
energetics to determine the favoured background. This can be done by comparing the 
on-shell actions for the two gravity backgrounds, which in the light of our earlier dis- 
cussions on finite temperature gauge-gravity dictionary corresponds to different phases 
of the dual gauge theory. It can be shown iT59ll that at temperature = 1/(277/24), 
there is a first order phase transition between these two backgrounds. At low tempera- 



tures Td < 1/ (27ri?4) the Euclidean version of ( 10 ) is energetically favoured, whereas at 
high temperature Td > 1/ {271 Ri) the background in ( 16 ) is favoured. Geometrically the 
energetics of this transition is fairly intuitive: it is a competition between the two circles 
along t and x^ and the smaller circles prefers to shrink. 

The gauge theory interpretation of this transition is fairly straightforward. A com- 
parison of the appropriately renormalized free energies in these two phases shows that 
in the low temperature phase the free energy scales as A'^^ whereas in the high temper- 
ature phase it scales as A"^. This is, according to our broad characterization, a confine- 
ment/deconfinement transition. Thus this model, as suggested by Witten in ref. [11591 , 
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is a holographic model dual to a non-supersymmetric gauge theory which undergoes a 
confinement/deconfinement transition, which is similar to QCD. 

1.2.4 The Type IIA theory: introducing flavours 

We will now discuss introducing flavour degrees of freedom in this model following 
refs. [11371 1138II which is known as the Sakai-Sugimoto model. We will call this the 
Witten-Sakai-Sugimoto model. Our approach will be very similar to the Type IIB con- 
struction discussed before. In this case flavour degrees of freedom are introduced by 
intersecting the Nc D4-branes by Nf D8-branes at = — ^ and Nf D8-branes at 
= ^ with the obvious constraint L < ttR^. This configuration can be represented by 
the following array. 








1 


2 


3 


4 


5 


6 


7 


8 


9 


D4 






















D8, D8 










• 













Table 2: The brane construction including fundamental flavours. 

From the open strings with one end attached to the D4-branes and the other end 
on the D8 or the D8-brane we obtain Nj fundamental fermions of opposite chirality. 
Since the D4-branes are intersected by the D8/D8-branes these fermions are necessarily 
massless. As seen from the D4-brane worldvolume perspective the flavour branes intro- 
duce a global flavour symmetry group U(A^/)l x \J{Nf)ii, which is exactly analogous 
to the flavour symmetry group in massless QCD. Henceforth we will call this the chiral 
symmetry. 

Since the brane construction here involves both branes and anti-branes, we may 
worry about the stability of the construction. Usually a brane-anti-brane pair is unstable 
since oppositely charged branes annihilate each other, which is manifest as the presence 
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of a negative mass squared tachyon in the spectrum of the open string stretching between 
the brane-anti-brane pair. The mass of this "tachyon" in this case is explicitly given by 



which is positive definite for L > n\'2a'. Therefore as long as the brane-anti-brane pair 

is separated by a length much larger than the string length scale, this system should be 
stable. We will again treat the flavour branes as probes in the background of D4-branes. 

1.2.5 The Type IIA theory: an instructive limit 

Before proceeding further, we will comment on a particularly instructive limit of the 
previous model. This has been realized recently in ref. [11 J and further generalized in 
refs. [l57l[T2l[T3l . Let us consider sending the radius of the compact direction to infin- 
ity making it a flat direction. The flavour brane-anti-brane pair are aligned (separated 
by a distance L in the x^'-direction) in exactly the same way as illustrated in table 2. We 
will focus on the low energy dynamics from the point of view of a (3 + 1) -dimensional 
observer living in the intersection region. We will concentrate in the limit gs — > 0, 
Nc oo with gsNc and Nf held fixed. 

In this special limit (R^ — )■ oo) we have three dimensionful parameters: the brane- 
anti-brane separation L, the five-dimensional 't Hooft coupling A5 and the string length 
is- The chiral fermions are located at the (3 + 1) -dimensional intersections but they 
interact via a (4 + 1) -dimensional gauge field. This interaction strength is represented 
by the dimensionless ratio A5/L. There is another dimensionless ratio given by X^/ig 
which is relevant in this model. In the regime of the following hierarchy of scales 




(17) 




A5 < 4 < L 



(18) 
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we can neglect all stringy effects. We can also neglect all gauge field dynamics on 
the worldvolume of the D4 and the D8-branes. The only relevant interaction is the 
exchange of (five-dimensional) gluons between the chiral fermions located at the brane- 
intersections. In this regime the theory is weakly coupled with a natural UV cut-off 
given by is- In the regime is ^ ^5 ^ L, the gauge theory is still weakly coupled, but 
the natural UV cut-off is given by A5. Finally in the regime ^ L ^ A5 the gauge 
theory becomes strongly coupled. In this limit we can reliably use the supergravity 
description (as demonstrated in ref. L74J) to address the strong coupling dynamics. 

It was further shown in ref. [[TT| that in the weak coupling regime the effective 
description of the gauge theory is given by a Nambu-Jona-Lasinio (NIL) model llll6B 
with a non-local four fermi interaction between the chiral fermions. The corresponding 
Lagrangian was obtained by writing down an effective action for a Weyl fermion located 
at the (3 + 1) -dimensional intersection and a (4 + 1) -dimensional SU(Ai'c) gauge field. 
After integrating out the gauge field, one arrives at the non-local NIL model with a non- 
local four fermi interaction term. This non-locality arises from the (4 + 1) -dimensional 
gauge field exchange between the chiral fermions. We will not elaborate more on this 
topic here, but refer the interested reader to ref. IfTTll . 

The NIL model was first proposed in ref. ] 116 ||. This model is an useful effective 
low energy description of QCD as well. It has been successfully used in several phe- 
nomenological studies of hadrons and therefore provide an interesting avenue to under- 
stand many aspects of QCD. For a review of some such applications see e.g. ref. [|93l . 
We will not use this model explicitly to provide quantitative estimates or comparisons 
on the topics we will study using the gauge-gravity duality, however, we will comment 
on the qualitative aspects where necessary. It is very exciting that the gravity dual of a 
certain kind of NJL model can be realized in a precise sense. 
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1.3 The probe brane action: the tool of the game 

Our general strategy will be to study the dynamics of a probe brane (of various dimen- 
sions) in an appropriate background. In subsequent chapters we will make heavy use 
of the probe action, so let us review some general aspects of this action. Here we will 
restrict ourselves only to the bosonic part of the probe action. 

The bosonic action consists of two parts: the Dirac-Born-Infeld (DBI) action and 
the Wess-Zumino (WZ) term. In supergravity we have two kinds of fields to which a 
D-brane can couple to, namely the Neveu-Schwarz (NS) fields (containing the metric, 
dilaton and the B-field) and the Ramond-Ramond (RR) fields (containing the various 
potentials). The DBI action describes the coupling of the D-brane with the background 
NS fields and the WZ action describes its coupling to the background RR fields. 

The DBI action was formulated by Born and Infeld[|27l to avoid various singularities 
of Maxwell's theory of electrodynamics and later studied in more detail by Dirac [|43l . 
This action was obtained by requiring Lorentz and gauge invariance and also the condi- 
tion that for small electromagnetic fields it should reduce to the familiar Maxwell action. 
The DBI action is also naturally equipped to provide a cut-off for the maximum electric 
field. This fact will play an interesting role in the later chapters. 

The DBI action was derived explicitly from string theory in refs. [|54l I102[ 11551 . 
Furthermore, the fermionic part of the DBI action has been analyzed in refs. Ill06[|107] 
I108II : however we will not need it. The bosonic part of the DBI action for a Dp-brane is 
given by 



where P [Gab + Bab] is the pull-back of the background metric and the B-field, Fab is 
the field strength corresponding to the gauge field Aa defined on the worldvolume of 





(19) 
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the probe brane. The constant fip is related to the tension of the D-brane and is given 
by[[l2l,/i, = (27r)-^a'-(P+i)/2. 

The Wess-Zumino action describing the coupling of the D-brane to the background 
RR potentials (denoted by C(p+i)) is a topological term. This action in complete gener- 
ality is given by (see e.g. ref. [[781 for a review) 



'S'wz — f'p 

JMp+i 



^ C'(P+1) 
. P 



(20) 



For most of the cases that we will consider in later chapters this term will not play any 
role in determining the brane's classical profile. It will, however, play an important role 
in the study of the fluctuation spectrum around this classical profile of the brane. 



1.4 Outline 

This thesis is organized as follows. We begin with a brief review of chiral symmetry 
breaking in QCD in chapter 2. As an effective model for QCD, we also review the so 
called NIL model and the pattern of chiral symmetry breaking within this model. Our 
focus is to obtain a qualitative understanding of the physics of chiral symmetry breaking 
under the influence of several external parameters such as temperature or an external 
electromagnetic field in the NJL model. 

Then we start discussing these aspects in the holographic set-up. In chapter 3 we 
first discuss introducing fundamental flavour in finite temperature A/" = 4 SYM theory. 
At finite temperature the flavour sector undergoes a first order meson melting transition. 
Now we introduce an external constant magnetic field and discuss its effect on the phase 
diagram of the fundamental flavours in finite temperature A/" = 4 SYM theory. We find 
that temperature and magnetic field are two competing parameters which give rise to an 
interesting phase diagram. We further realize that an external magnetic field can induce 
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spontaneous breaking of chiral symmetry which was first observed in ref. [53 1 in the 
zero temperature set-up. We begin by establishing the non-trivial phase structure that 
results from finite temperature. We observe, for example, that above the critical value of 
the field that generates a chiral condensate spontaneously, the meson melting transition 
disappears, leaving only a discrete spectrum of mesons at any temperature. We also 
compute several thermodynamic properties of the plasma. 

Following the same logical path, we investigate the effect of an external electric 
field in chapter 4. At zero temperature, we observe that the electric field induces a 
phase transition associated with the dissociation of the mesons into their constituent 
quarks. This is an analogue of an insulator-metal transition, since the system goes 
from being an insulator with zero current (in the applied field) to a conductor with 
free charge carriers (the quarks). At finite temperature this phenomenon persists, and 
the dissociation transition becomes subsumed into the more familiar meson melting 
transition. Here, the dissociation phenomenon reduces the critical melting temperature. 

We change gears in chapters 5 and 6 and work in a different model. In chapter 5 using 
the Witten-Sakai-Sugimoto model we study the effect of an external magnetic field on 
the dynamics of fundamental flavours in both the confined and deconfined phases of a 
large Nc gauge theory. We find that an external magnetic field promotes chiral sym- 
metry breaking, consistent with the "magnetic catalysis" observed in the field theory 
literature, and seen in the Type IIB set-up in chapter 3. The external field increases 
the separation between the deconfinement temperature and the chiral symmetry restor- 
ing temperature. In the deconfined phase we investigate the temperature-magnetic field 
phase diagram and observe, for example, there exists a maximum critical temperature 
(at which symmetry is restored) for very large magnetic field. We find that this and 
certain other phenomena persist for the Witten-Sakai-Sugimoto type models with probe 
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branes of diverse dimensions. We comment briefly on the dynamics in the presence of 
an external electric field. 

In chapter 6 we continue our study of the dynamics of the flavour sector of the 
Witten-Sakai-Sugimoto model in the presence of an external magnetic field, uncovering 
several features of the meson spectrum at high and low temperatures. We employ both 
analytical and numerical methods to study the coupled non-linear equations that result 
from the gravity dual. Finally in chapter 7 we conclude with open problems and future 
directions. 
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Chapter 2: Review of Chiral 
Symmetry Breaking 

In this chapter we briefly review the phenomenon of chiral symmetry breaking in QCD 
and the NIL model which is an effective model of QCD. We also discuss the effect of 
an external parameter (such as temperature and an external electromagnetic field) on the 
chiral symmetry breaking in the latter model. Our discussions here are entirely based 
on the more detailed reviews in ref. [11391 [681 . 

2.1 QCD and Symmetries 

We begin with a brief review of the symmetries of QCD. We will closely follow the 
notations and discussions in ref. [[139L The fundamental degrees of freedom of QCD 
are the fermionic quark fields (denoted by ip^) and the non-abelian gauge fields called 
gluons (denoted by Ap. We denote flavour indices by the letter i and color indices by 
the letter a whereas jj, labels the spacetime indices. The quarks transform under the 
fundamental representation of the SU(Ai'c) gauge group (where Nc = 3) and the gluons 
transform under the adjoint representation. The QCD Lagrangian is given by 

£qcd = iJ{iip-m)i;- ^F^^F'^'^^ , (21) 
where {d, - igAp^) , F;, = d,At - d.A^ + gr'^A^Al , (22) 

where g is the gauge coupling, are the generators of the SU(3) gauge group and /"'"^ 
are the structure constants. In writing the fermionic term in the Lagrangian we have 
implicitly assumed a summation over the flavour indices. 
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To begin with, the QCD Lagrangian is invariant under a local SU(3) gauge symmetry 

i/j{x) [/(x)V(x) , ^ U{x)Ap''U^{x) + iU{x)d^JJ^{x) . (23) 

If the masses of the three lightest quarks are the same then the Lagrangian is invariant 
under the following global flavour rotations 

= Viji^j , V e SU(3) , (24) 

which is known as the flavour isospin symmetry. Furthermore if the masses are equal to 
zero, then this global symmetry is enlarged. Let us define the chiral quark fields as 

V'L,i? = ^(l±75)V', (25) 
and rewrite the fermionic part of the QCD Lagrangian as 

>Cfermion = i^L^Ip^L + ^'l^^'R + ilRlIptpR + ^'uMlpL , (26) 

where M is the diagonal mass matrix. Clearly if M = 0, the Lagrangian is invariant 
under the following independent flavour transformations 

^L,^ ^ VL,^J^/JL,j , ^R,i ^ V^i^,^,^R,, where {Vl, Vn} G SU(3)^ x SU(3)r . (27) 

This is known as the chiral symmetry. In the real world the quark masses are not van- 
ishing. However, for the three lightest quarks the masses are negligible compared to the 
scale Aqcd and therefore it is an approximate symmetry of the theory. 
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Lastly, the QCD Lagrangian has the following two global U(l) symmetries as well 



U(1)b : i'L^R ^ e"f'ijL,R , (28) 
U(1)a: ^L^e^>L, ^fi^e-^>^, (29) 

where cp and a are constant angles. The U(1)b is an exact symmetry for massive quarks 
as well. On the other hand, the axial U(1)a is exact for the massless quarks only at the 
classical level and is broken by quantum anomalies. This can be represented by writing 
the divergence of the axial current 

d'jl = ^F^^F'"'' , F^''" = l^'^'^'^F^p . (30) 

At zero temperature and density, the QCD vacuum breaks chiral symmetry via the 
formation of a chiral condensate (ipip). The symmetry group breaks down to SU(3)l x 
SU(3)r —7- SU(3)y and an octet of approximately massless pseudoscalar Goldstone 
bosons appear. Chiral symmetry restricts the interactions of these Goldstone modes 
which can be obtained from the low energy effective chiral Lagrangian. The low energy 
effective chiral Lagrangian is given by 

Ahirai = ^Tr [a^Ea'^St] + [BTr (MS^) + h.c] + . . . , (31) 

where is the pion decay constant, S collectively denotes the Goldstone modes, 
M is the quark mass matrix and i? is a constant. The Goldstone modes associated 
with the chiral symmetry breaking satisfies the Gell-Mann-Oakes-Renner (GMOR) 
relationship dsn 



^l^-pii^i^), (32) 
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where m is the bare quark mass. 

At finite temperature, QCD undergoes the deconfinement and the chiral phase tran- 
sition. However, close to this critical temperature the essential dynamics remain non- 
perturbative. We will now turn to an effective model of QCD known as the Nambu- 
Jona-Lasinio (NJL) model. It turns out to be a useful toy model for understanding the 
physics of chiral symmetry breaking among other features of QCD. We will discuss the 
effect of external parameters on chiral symmetry breaking within this model. 

2.2 The Nambu-Jona-Lasinio model 

This model was originally proposed in ref. [I116II . The Lagrangian of the NJL model is 
given by 



(33) 



where G is the coupling corresponding to the four-Fermi interaction. For simplicity we 
will focus only on the single flavour case. Due to the presence of the four-Fermi terms, 
this is a non-renormalizable theory and therefore requires a UV cut-off. In writing the 
four-Fermi interaction terms we have also assumed a summation over the color indices. 

This bare Lagrangian is invariant under anU(l)/, x U(l)ij symmetry (which from 
now onwards we refer as the chiral symmetry) 



'4}L,R exp ( '^-eL,R j ^ , ^L,R exp { -^0L,R ] ^ , (34) 



where ipL,R is defined as in equation (25 1. Clearly the NJL model has an obvious gener- 
alization for Nf flavours. 
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There are two simple cases to consider: when Ol = 6^ = Oy, the chiral symmetry 
group is the U(l)y group; the other choice is 9l = —Or = 9a for which the invariance 
group is the axial U(1)a group. 

Now we want to couple the NJL model with an external electromagnetic field. The 
corresponding Lagrangian is then given by 

£njl = i^ilpif^ + G 

Note that we can rewrite this Lagrangian by introducing auxiliary fields (denoted by a 
and tt) 

£njl = i^iiPip -i^{(T + i7^7r) ^ - ^ (o^' + > (36) 
where a = -2G {^tp) , n = -2G {^ij^ip) . (37) 

To analyze the effect of the external field on the chiral symmetry we have to study what 
happens to the dynamically generated mass associated with the chiral symmetry break- 
ing. One way to achieve this is to compute the one-loop Coleman-Weinberg effective 
potential [|34l . We will not review the calculation of this potential, but point the inter- 
ested reader to the relevant reference, e.g. ref. [|68l . A pedagogical introduction to 
Coleman-Weinberg potential can be found in e.g. ref. 01241 . The effective potential for 
the NJL model, denoted by Vcs, depends on the U(1)l x U(1)/j invariant combination 
= (j^ + 7r^. Therefore without any loss of generality, we consider a configuration with 
TT = 0. We will now review the effective potential for this model under various external 
parameters, such as temperature and constant electromagnetic field. 



D, 



du, - ieA, 



(35) 
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2.2.1 The effective potential at finite temperature 

Let us first review the properties of the ground state by analyzing the effective potential 
at finite temperature. The effective potential at finite temperature is given by (see e.g. 

ref. ttnni) 

Kfr(p) = ^-^^^£ + 2Tlog (^1 + exp (^-1 v7+7) ) 

(38) 

where T is the temperature. 

The true vacuum of the system is given by dVes/dp = 0. This is popularly known 
as the Gap equation. At T = 0, the Gap equation gives [|68l 



pA' (- - 1 j = -P' [-) + ^( VA^) , 9 = ■ (39) 

Clearly in order for the above equation to have a real non-trivial solution one should 



have g > Qc = 1. Also, any solution of the Gap equation (39), p = m^y^ basically 
gives the dynamically generated quark mass (from the definition of a in the Lagrangian 
in ( [36] )) associated with the chiral symmetry breaking. In this case, the chiral symmetry 
is spontaneously broken if g > q^ [|68l . 

At finite temperature the resulting Gap equation is given by 

P N.P r" ^^^^ 



We will not explore this equation in detail, but we will try to estimate the effect of 
introducing a temperature. Specifically we want to estimate how the coupling constant 
must be tuned if we introduce an infinitesimally small temperature keeping the value of 
p fixed. 



37 



First let us note that a solution of equation (39), p ~ 0{JS). Therefore for small 
temperature, the exponential factor in the integrand in (|40]) can be well approximated by 
a small positive constant, which we denote by e. Further we assume that the change in 
the couphng is given hy g = qiq) + 5g, where g{Q) is the coupling at zero temperature. 
In this limit the Gap equation becomes 

pA^ Q-(l-e)) = V(l-6)log (^) 
^ 5g = tg^Q) > , (41) 

where we have used the fact that p satisfies the Gap equation at vanishing temperature. 
This suggests that in order to have a chiral symmetry breaking the four Fermi coupling 
needs to be stronger in finite temperature background. In other words, finite temperature 
tends to restore chiral symmetry. Physically this is intuitive to understand. Raising tem- 
perature will presumably melt the quark bound states and break the quark condensate, 
therefore tend to restore the chiral symmetry. For more detailed analysis we refer the 
reader to the review in ref. [|93l. 



2.2.2 The effective potential in a magnetic field 

In the presence of an external constant magnetic field the effective potential is given 
byil 

0^ N f°° ds 9 
Kfr(p) = j^ + ^ ^^"'^H coth (eHs) . (42) 
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The corresponding Gap equation in this case is given bv ll68l 



pA' 



= -p^ log 



2 

pH^ 



+ ip' + Y, log 



+ 0(1/A) 



47r 



(43) 



where 7 is the Euler constant and £^ = {eH) ^ is the magnetic length. 

The magnetic field modifies the Gap equation in a rather interesting way. To see this, 



let us assume that p = mdyn is a solution of (43) and p 
Further, let us \ 
can obtain ||68l 



.(0) 



m^^yn is a solution of (|39|). 

Further, let us work in the small magnetic field limit, i.e. m^^^ ^ i. Then from (|43|) we 



dyn 



m 



(0) 
dyn 



1 + 



{eHf 



3(miyl) log(A/mi;g' 



(44) 



which means that the magnetic field enhances the breaking of the chiral symmetry. 

Interestingly, the Gap equation in the presence of a magnetic field admits a solution 
even for[68J g <^ gc given by 



-dyn - efep J 



(45) 



It is therefore clear that any non-zero magnetic field in this model induces a spontaneous 
chiral symmetry breaking for arbitrarily small values of the coupling. This feature of 
a magnetic field inducing spontaneous breaking of chiral symmetry, however, is rather 
universal and commonly known as the magnetic catalysis in chiral symmetry breaking. 
In section 2.3, we will review the basic mechanism behind this catalysis. 
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2.2.3 The effective potential in an electric field 

At zero temperature, the unbroken Lorentz invariance particularly simplifies the expres- 
sion for the effective potential. Specifically, we can obtain the effective potential in a 
constant electric field from the potential in a constant magnetic field just by sending 



H — 7- iE, and vice versa, in e.g. equation (42) 



The consequences are simple to observe. Now, the dynamically generated mass 



mdyn < "^dyn vi^tue of equation (44). Hence we conclude that an electric field tends 



to restore the chiral symmetry. Physically this is expected, since an external electric 
field will tend to dissociate the quarks bound states into its constituents. 



2.3 Fermions in a constant magnetic field 

In this section we will review the basic properties of a relativistic fermion in a constant 
magnetic field in (3 + 1) -dimensions. We will closely follow the discussions in ref. [|68l . 
The corresponding Lagrangian is simply given by 

-Cfermion = ^ - m) ijj , D^ = d^- icA^, . (46) 

In the so called Landau gauge, the gauge potential corresponding to a constant external 
magnetic field (pointing in the -direction) is given by 

= -S^iHx^ . (47) 

The energy spectrum of the fermions is given by the Landau levels 

En ih) = ±Jm'^ + kl + 2\eH\n , where n G Z+ U {0} . (48) 
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Clearly, n = corresponds to the lowest Landau level (henceforth abbreviated as LLL). 
The LLL has degeneracy and the number of states is given by 



dNo = V 



\eH\ dks 
"2^ 2^ ' 



(49) 



where V is the spatial 3-volume. 

Let us first argue that the LLL plays a special role in the infrared dynamics. Clearly 
the energy of the LLL is given by a/ m'^ + fc| whereas the energy for the higher Landau 



levels goes as ^y\eH\. It is therefore obvious that the infrared dynamics in a strong 
magnetic field, i.e. when \/\eH\ ^ m, ^3, is insensitive to the higher Landau levels. 

We want to demonstrate that the dynamics of the LLL plays a crucial role in cat- 
alyzing the chiral symmetry breaking. To do so, we consider the fermion propagator in 



a background constant magnetic field computed in ref. 11141 



S (x, y) = exp 



{x-y fAf, (x + y) 



S{x-y) 



(50) 



where 



S{k) 



dsexp 



is ( kl 



2 ,2 tan [eHs) 



eHs 



— m 



[(fco7° - ^37^ + m) (1 + tan (eHs)) - k±-f± (l + tan^ (eHs))] 



(51) 



where k± = {ki,k2} and 'y± = {71,72}- Now using expression for the generating 
functional for Laguerre polynomials (denoted by L") 



xz 
z-l 



n=0 



[X)Z 



(52) 
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we can decompose the fermionic propagator over the Landau poles 



S{k) = iexp 



leH 



n=0 U 



Dn {eH, k) 



A;| — — 2\eH\n 



^^(")(A;), (53) 



n=0 



where 



Dn{eH,k) = (fco7°-A;37' + H 



(l-^7Vsign(eif))L0 



\eH\ 



-(l + z7Vsign(ei/))L0„i 



2kl 
\eH\ 



+ 4 {ka' + ^27') 



2kl 
\eH\ 



The contribution coming from the LLL is then given by 



(54) 



5(°)(A;) 



zexp 



\eH\ 



^o7° — ^37^ + ^ 



(l — 'i7"'^7^sign (eiJ)) 



(55) 



It is now clear from the LLL contribution to the propagator that in the infrared regime, 
i.e. when k'j_ <^ \eB\, S^'^^ is independent of k± completely and the dynamics is effec- 
tively (1 + 1) -dimensional. The physical reason behind this is simple: the motion of a 
charged particle is constrained in the plane perpendicular to the magnetic field. 

It is well-known that infrared dynamics is stronger in low dimensional field theo- 
ne ^ This can be qualitatively understood as follows: the measure in the momentum 
space goes as k^dk in d = 3 and kdk in d = 1. Clearly in the infrared, i.e. for small 
k, k^ suppresses the dynamics much stronger than k. Therefore it is expected that the 
infrared dynamics in (1 + 1) -dimensional theory will remain much strongly coupled as 
compared to the infrared dynamics of a (3 + 1) -dimensional theory. This therefore sug- 
gests an universal mechanism through which an external magnetic field facilitates the 



We thank V. Miransky for a discussion on this point. 
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spontaneous breaking of chiral symmetry by generating a dynamical quark mass. In the 
forthcoming chapters we will further observe that a similar universal mass generation 
mechanism in the presence of a magnetic field exists within the specific holographic 
models that we consider in this thesis. 
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Chapter 3: Flavoured AT = 4 



Yang-Mills theory at finite temperature 
and an external magnetic field 

3.1 Introductory remarks 

We begin with a holographic analysis of thermal properties of large Nc flavoured Yang- 
Mills theory. We will then consider the effect of both temperature and an external mag- 
netic field in later sections. Our holographic set-up is the so called Type IIB frame- 
work reviewed in the first chapter. The presented material about the thermal proper- 
ties is based on a close collaboration with Tameem Albash, Veselin Filev and Clifford 
Johnson[7]. Similar results were obtained in ref. UllOII and also in ref. [85] . The main 
result, however, was already known before in ref. ||89l . 

We study the geometry of AdSs-Schwarzschild xS^ which is the decoupled/near- 
horizon geometry of Nc D3-branes, where Nc is large and set by the (small, for relia- 
bility) curvature of the geometry. The physics of closed Type IIB string theory in this 
background is dual to the physics of A/" = 4 supersymmetric SU(Ai'c) gauge theory in 
four dimensions, with the supersymmetry broken by being at finite temperature !! 1 5911 . 
The temperature is set by the horizon radius of the Schwarzschild black hole, as we will 
recall below. 

We introduce a D7-brane probe into the background. Four of the brane's eight 
world-volume directions are parallel with those of the D3-branes, and three of them 
wrap an C S^. The remaining direction lies in the radial direction of the asymptoti- 
cally AdSs geometry. 
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Such a D3-D7 configuration controls tlie pliysics of tlie SU(Ai'c) gauge tlieory witli 
a dynamical quark in the fundamental representation [|83l . The configuration (at zero 
temperature) preserves J\f = 2 supersymmetry in D = A, and the quark is part of a 
hypermultiplet. Generically, we will be studying the physics at finite temperature, so 
supersymmetry will play no explicit role here. 

We are studying the D7-brane as a probe only, corresponding to taking the Nc S> 
Nf limit, and therefore there is no backreaction on the background geometry. This is 
roughly analogous to the quenched approximation in lattice QCD. The quark mass and 
other flavour physics-such as the vacuum expectation value (vev) of a condensate and 
the spectrum of mesons that can be constructed from the quarks-are all physics which 
are therefore invisible in the background geometry. We will learn nothing new from the 
background; our study is of the response of the probe D7-branes to the background, and 
this is where the new physics emerges from. 

We carefully study the physics of the probe itself as it moves in the background 
geometry. The coordinates of the probe in the background are fields in an effective D7- 
brane world-volume theory, and the geometry of the background enters as couplings 
controlling the dynamics of those fields. One such coupling in the effective model rep- 
resents the local separation, L(u), of the D7-brane probe from the D3-branes, where u 
is the radial AdSs-Schwarzschild coordinate. 

In fact, the asymptotic value of the separation between the D3-branes and D7-brane 
for large u yields the bare quark mass nig and the condensate vacuum expectation value 
(vev) (tptp) as follows [|99l[T28]l: 

lim L{u) = m H — - + . . . , (56) 
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where m = 27ia'mg and — c = {ijjil)) / {Sti^ a' N jt-j) , in this chapter we will set Nf = 1. 
The fundamental string tension is defined as T = l/{2na') here, and the D7-brane 
tension is = (2n)^'^(a')~^. The zero temperature behavior of the D7-branes in the 
geometry is simple. The D7-brane world-volume actually vanishes at finite u, corre- 
sponding to the part of the brane wrapped on the shrinking to zero size. The location 
in u where this vanishing happens encodes the mass of the quark, or equivalently, the 
separation of the probe from the D3-branes. In addition, in the zero temperature back- 
ground, the only value of c allowed is zero, meaning no condensate is allowed to form, 
as is expected from super symmetry. 

The finite temperature physics introduces an important new feature. As is standard 
Il59|, finite temperature is studied by Euclideanizing the geometry and identifying the 
temperature with the period of the time coordinate. The horizon of the background 
geometry is the place where that shrinks to zero size. The D7-brane is also wrapped 
on this S^, so it can vanish at the horizon, if it has not vanished due to the shrinking 
of the S^. For large quark mass compared to the temperature (horizon size), the 
shrinking will occur at some finite u > uh, and the physics will be similar to the zero 
temperature situation. However, for small quark mass, the world-volume will vanish 
due to the shrinking of the corresponding to the D7-branes going into the horizon. 
This is new physics of the flavour sector. 

The authors of ref. [ 15J explored some of the physics of this situation (the depen- 
dence of the condensate and of the meson mass on the bare quark mass), and predicted 
that a phase transition should occur when the topology of the probe D7-brane changes. 
However, they were not able to explicitly see this transition because of poor data res- 
olution in the transition region, coming from using UV boundary conditions on the 
scalar fields on the D7-brane world-volume. The origin of this phase transition, as 
we shall see, is as follows: The generic behavior of an allowed solution for L(u) as in 



46 



equation (56), is not enough to determine whether the behavior corresponds to an S^- 
vanishing D7-brane or an S^-vanishing D7-brane. The choices of branch of solutions 
have different values of c, generically. In other words, for a given value of m there 
can be more than one value of c. There are therefore two or more candidate solutions 
potentially controlling the physics. The actual physical solution is the one which has 
the lowest value for the D7-brane's free energy. The key point is that, at a certain value 
of the mass, the lowest energy solution may suddenly come from a different branch, 
and, as the corresponding value of the condensate changes discontinuously in moving 
between branches, we find that the system therefore undergoes a first order phase tran- 
sition. On the gauge theory side, we can imagine a similar situation occurring; two 
different branches of solution are competing, and the lowest energy branch is always 
picked. We are able to uncover this physics by doing a careful numerical analysis of 
the equations of motion for the probe dynamics on the gravity side of the AdS/CFT 
correspondence, by using IR boundary conditions instead of UV boundary conditions. 



3.2 Holographic meson melting 
3.2.1 General set up 

We begin by reviewing the physics of the D7-brane probe in the AdS5-Schwarzschild 
background solution [[151. The background metric is given by: 

ds^ = -^^df + -^du^ + ^dx-dx + R^dnl , f{u) =u^ ^ 
where d^l = dO^ + cos^ Odnj + sin^ ^^0^ , 

and dill = di/j"^ + cos^ ipd/d'^ + sin^ ipd'y'^ , (57) 
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where x is a three vector defining the three spacelike directions of the dual gauge theory, 
u G [0, oo) is the radial coordinate and we are using standard polar coordinate on the 
S^. The quantity is given by: 



where gs is the string coupling (which, with the inverse string tension a' sets for 
example, Newton's constant). The quantity b is related to the mass of the black-hole, 
= 8G'5mb.h./(37r). The temperature of the black hole can be extracted using the stan- 
dard Euclidean continuation and requiring regularity at the horizon. Doing this in the 



metric given by equation (57), we find that (3^^ = b/nR^. Therefore, by picking the 
value of b, we are choosing at what temperature we are holding the theory. 

Later we will use different coordinate system for convenience, so here we introduce 
the change of variables: 

=i(n2 + vV^)=p2 + L2, (58) 
with p = r cos 6 , L = r sin 6 . 

The expression for the metric now takes the form: 



•2 



+ ^ (dp"" + p^dnl + dL^ + L'dcl?) . (59) 



Note that with this variable change r — t- m as n — )• oo. 

We choose to embed the D7-brane probe transverse to and (p ala the background 



in (57). In order to study the embeddings with the lowest value of the on-shell action 



(and hence the lowest free energy), we choose an ansatz of the form = and 9 = 9{u). 
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The asymptotic separation of the D3- and D7-branes is given by L{u) = u sin 9. Given 
this particular choice of embedding, the world- volimie of the D7-brane is given by: 



= u'^ cos^ e{u)Vdet S^y/u^ + {u^ - b^) e'{uy , (60) 

where g is the determinant of the induced metric on the D7-brane given by the pull-back 
of the space-time metric G^i,. We are interested in two particular cases. First, there is 
the case where u goes to b, which corresponds to the D7-brane probe falUng into the 
event horizon. In the Euclidean section, this case corresponds to the shrinking of the 

of periodic time, we will refer to these solutions as the 'black hole' embeddings. 
Second, there is the case where 9 goes to 7r/2, which corresponds to the shrinking of the 
S^. We will name these Minkowski embeddings. It is this change in topology (S^ versus 

shrinking) between the different solutions that will correspond to a phase transition. 
The classical equation of motion for 9{u) is: 



d I u^u^ - b^) 9' cos^ 9 
du 



+ 3u'cos^9sm9Ju^ + {u'^ -b^)9" = 0. (61) 



When u goes to infinity and the background metric becomes asymptotically AdSs x S^, 
the equation of motion reduces to: 

-^(u^9'(u))+3u^9(u) = 0, (62) 
du 

which has solution: 

9{u) = - fm + 4) • (63) 
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These two terms are exactly the non-normalizable and normalizable terms correspond- 
ing to a dimension 3 operator (tp-ip) in the dual field theory with source m and vacuum 
expectation value (vev) c. 



Next we solve equation (61 ) numerically using a shooting technique. Shooting from 
infinity towards the horizon, physical solutions are those that have a finite value at the 
horizon. This will only be accomplished for a particular m and c value from equation 



(56) which would have to be delicately chosen by hand. Therefore, if we instead start 
from the horizon with a finite solution, it will shoot towards the physical asymptotic 
solutions we desire. In order to be able to analyze the phase transition between the 
black hole and Minkowski solutions, we shoot from the horizon for the condensate 
solutions and from 9 = 7r/2 for the Minkowski solutions. This technique avoids having 
to correctly choose the boundary conditions at infinity, which is a sensitive procedure, 
allowing us to have many more data points to analyze the phase transition. We impose 
the boundary condition that, at our starting point-the horizon, we have: 



We argue that this is the physical boundary condition to take; the first is simply a result 
of taking the limit of n — )• 6 in the equation of motion in equation ( [6T] ), whereas the 
second is a result of requiring no conical singularity as the shrinks to zero size [|85l . 

For the sake of visualization, we will present the profile of the probe branes in 
{p, L}-plane. Therefore we define the following dimensionless quantities 

p=-, L = - , "^ = y' (65) 
L = m + + . . . . (66) 
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In this coordinate system several D7-brane embedding solutions are shown in Figure [T] 
The red (solid) lines correspond to Minkowski solutions, and the blue (dashed) lines 
correspond to black hole solutions. From each of these solutions, we can extrapolate the 
bare quark mass and quark condensate vev. 

L{~p) 




A ' p 

Figure 1: Solutions for the D7-brane probe in the black hole background. 



3.2.2 The first order phase transition 

We plot the c values as a function of m in Figure |2(a) When enlarged, as shown in 



Figure 2(b), we find, as anticipated in the introductory remarks, the multi-valuedness 
in c for a given m. Physics will choose just one answer for c. There is therefore the 
possibility of a transition from one branch to another as one changes m. 

In order to determine exactly where the transition takes place, we have to calculate 
the free energy of the D7-brane. In the semi-classical limit that we are considering, the 
free energy is given by the on-shell action times 13'^. For our case, this is simply given 
by [193: 

jr = p-^TjNf j (fx dn-i du ^- det g , (67) 
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where here Nf = 1. Akematively we can determine the phase transition point by con- 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 

(a) 



0.918 0.920 0.922 0.924 0.926 0.928 
(b) 



Figure 2: The condensate vev at the phase transition. See text for more details. 



sidering an equivalent Maxwell construction. Since c and m are thermodynamically 
conjugate to each other, the free energy difference can be expressed as dj' ~ —cdm. 
Therefore the phase transition will happen exactly at the point where the areas under 



the curve separated by the vertical dashed line in Figure 2(b) are equal. Thus we find 
a first order phase transition — at m ^ 0.92345 — where the condensate's vev jumps 



discontinuously, which is shown with the dashed black line in figure 2(b) 



This first order phase transition was studied in details in ref. [73], where the authors 
developed an appropriate formalism to study the meson spectrum of the theory before 
and after the phase transition. They were able to show that after the phase transition 
there are no bound meson states, but only quasi-normal excitations corresponding to a 
melting mesons. This is why the present understanding of the observed first order phase 
transition is that it corresponds to the confinement/deconfinement phase transition of 
the fundamental matter which in the case of consideration is the meson melting phase 
transition. 
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3.3 Introducing an external magnetic field 

We will investigate the effect of an external constant magnetic field in the dynamics 
of probe flavours in the thermal N = A Yang-Mills theory. The material presented 
here is based on the work done in collaboration with Tameem Albash, Veselin Filev 
and Clifford JohnsonflH. Similar work was also done in ref. [|49l . Ref. [|53l analyzed the 
dynamics of probe flavours in a background constant magnetic field at zero temperature, 
which is the premise of the finite temperature physics we will discuss in this chapter. 

It was shown in ref. [[53l that the presence of an external magnetic field induces 
a spontaneous chiral symmetry breaking by generating a non-zero chiral condensate, 
where the chiral symmetry is an U(1)r symmetry as explained in chapter 1. However 
it was described in ref. [fTSfl that finite temperature melts the mesons and restores this 
symmetry at zero bare quark mass. Therefore we have two competing parameters and 
we will investigate their effects on the phase structure of the fundamental flavours. 

We will consider similar embedding of the probe brane represented by = 0, L = 
L(p). An external constant magnetic field can be introduced in this set-up by exciting 
an Abelian gauge field on the worldvolume of the probe brane [|53lj ^ 

^3 = Hx"^ , (68) 



where is a real constant. Exciting a non-trivial gauge field on the worldvolume of 
the probe brane (resulting in a non-trivial value for Fab) should not backreact on the 
background in the probe limit, however it will affect the probe brane. Also note that 



exciting a constant anti- symmetric field as in (68) on the worldvolume of the probe 



brane corresponds to exciting a non-normalizable mode in the boundary gauge theory. It 



"^It is easy to check a posteriori that this ansatz is consistent with the equations of motion and the 
Bianchi identity for the field strength of this gauge potential. 
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is fairly evident from the form of (68 1 that this anti- symmetric field should be identified 
with an external constant magnetic field (the strength of which is denoted by H). 

On the gauge theory perspective, this magnetic field therefore couples only to the 
flavour sector and does not affect the background adjoint fields. We will see that this 
magnetic field affects the finite temperature physics in many interesting ways. 

To study the effects on the probe, let us consider the general (Abelian) DBI action: 

Sdbi = -NjTdj J rf'e dee/\P[Gab + Bab] + 2Tra'Fab) , (69) 

Ma 

where Tor = /iy/fi's = [(2vr)'''a''^(7s]~^ is the D7-brane tension, P[Gab] and P[Bab] are 
the induced metric and induced fi-field on the D7-branes' world-volume. Fab is the 
world-volume gauge field, and Nj = 1 here. The resulting Lagrangian is: 



16(p2 + L(p)2)V [ ■ (64 + 4 (p2 + ^(^)2)2y 



(70) 

Note that at zero temperature the corresponding Lagrangian is given by 



which was analyzed in detail in ref. ||53l . For large p b, the Lagrangian asymptotes 
to: 



C^-p^^/l + lTW , (72) 
which suggests the following asymptotic behavior for the embedding function L{p): 



L{p) = m + 4 + • ■ ■ , (73) 
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where the parameters m (the asymptotic separation of the D3- and D7-branes) and c 
(the degree of transverse bending of the D7-brane in the (p, </)) plane) are related to the 
bare quark mass ruq = m/2Txa' and the quark condensate {ijjil)) oc — c respectively [|99l . 
To proceed, it is convenient to define the following dimensionless parameters: 



P 

m 

This leads to the Lagrangian: 



b ' 
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H 



m 



m 
J ' 



L{bp) 



c 



m + — + . . . 

p2 



(74) 



C 



-p' 



16 p2 + ^(p)2 



16 p2 + L(p)2 



1 + 4 p2 + ^(p)2 



1 + L'{PY . 



(75) 



For small values of i], the analysis of the second order, non-linear differential equation 



for L{p) derived from equation (75 ) follows closely that performed in refs. iTTl fTSlfTTOl . 
The solutions split into two classes: the first class are solutions corresponding to embed- 
dings that wrap a shrinking 5''^ in the 5*^ part of the geometry and (when the vanishes) 
closes at some finite radial distance r above the black hole's horizon which is located at 
r = h/\/2. These embeddings are referred to as 'Minkowski' embeddings. The second 
class of solutions correspond to embeddings falling into the black hole, since the of 
the Euclidean section, on which the D7-branes are wrapped, shrinks away there. These 
embeddings are referred to as 'black hole' embeddings. As we have demonstrated in the 
previous section, this gives rise to a first order melting transition at some critical mass 
rricT- We will now show that the effect of the magnetic field is to decrease this critical 
mass, and, at some critical magnitude of the parameter rjcr, the critical mass drops to 
zero. For r] > rjcr the phase transition disappears, and only the Minkowski embeddings 
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are stable states in the dual gauge theory, possessing a discrete spectrum of states corre- 
sponding to quarks and anti-quarks bound into mesons. Furthermore, at zero bare quark 
mass, we have a non-zero condensate and the chiral symmetry is spontaneously broken. 



3.4 Properties of the solution 
3.4.1 Exact results at large mass 

It is instructive to first study the properties of the solution for m ^ 1 . This approxima- 
tion holds for finite temperature, weak magnetic field, and large bare quark mass m, or, 
equivalently, finite bare quark mass m, low temperature, and weak magnetic field. 

In order to analyze the case m ^ 1, let us write L(p) = d + for a ^ 1 and 
linearize the equation of motion derived from equation ( [75] ), while leaving only the first 
two leading terms in (p^ + m^)^^. The result is: 

2(n'^ + 1)2-1 

dpip'C) - + ^^i-^m + 0(C) = . (76) 



Ignoring the 0(C) terms in equation (76), the general solution takes the form: 



4:p^[m^ + p'^) 9bp^{m^ + p'^y^ 

where we have taken C'(0) = C(0) = 0. By studying the asymptotic behavior of this 
solution, we can extract the following: 

ri^ 1 + 4r?2 + 2r7^ ^ / 1 

" = "-4^+ 32y +'^(y 

c = t- '^'tt'"\ o(^) . (78) 
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By inverting the expression for m, we can express c in terms of m: 

rf 1 + Arf + 8r]^ f ^ 



+ O — . (79) 



4m 96m^ \fh? 

Finally, after going back to dimensionful parameters, we can see that the theory has 
developed a fermionic condensate: 

oc -c = --i?^ H- . (80) 

The results of the above analysis can be trusted only for finite bare quark mass and 
sufficiently low temperature and weak magnetic field. As can be expected, the physically 
interesting properties of the system should be described by the full non-linear equation 
of motion of the D7-brane. To explore these we need to use numerical techniques. 

3.4.2 Numerical analysis 



We solve the differential equation derived from equation (75 ) numerically. It is conve- 
nient to use infrared initial conditions [El 13. For the Minkowski embeddings, based on 
symmetry arguments, the appropriate initial conditions are: 

L(p)|^=o = i^in, L'(p)|p=o = 0. (81) 

For the black hole embeddings, the following initial conditions: 

L 



L(p)\e.h. = Lin , -^'(P)le.h. 

P 



(82) 

e.h. 
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ensure regularity of the solution at the event horizon. After solving numerically for 
L{p) for fixed value of the parameter 7], we expand the solution at some numerically 



large pmax, and, using equation ( [73[ ), we generate the plot of — c vs m. 

In section 3.2, we have discussed the nature of the phase transition and the associated 
Maxwell construction for the thermodynamic free energy. Here we will follow the same 
recipe in the presence of the additional parameter t] which controls the strength of the 
magnetic field. 

Let us turn on a magnetic field. As one can see from Figure [3} the effect of the 
magnetic field is to decrease the magnitude of rhcr- In addition, the condensate now 
becomes negative for sufficiently large fh and approaches zero from below as m — )• oo. 



It is also interesting that equation (79 1 is still a good approximation for m > rhcr- 

For sufficiently strong magnetic field, the condensate has only negative values and 
the critical value of m continues to decrease, as is presented in Figure |3| If we further 



c(m) 




(a) 



(b) 



Figure 3: For strong magnetic field the condensate is negative. The value of rhcr continues to 
drop as we increase ry. 



increase the magnitude of the magnetic field, some states start having negative values of 
m, as shown in Figure |4} The negative values of m do not mean that we have negative 
bare quark masses; rather, it implies that the D7-brane embeddings have crossed L = 
at least once. It was argued in ref. [[TSlI that such embeddings are not consistent with a 
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holographic gauge theory interpretation and are therefore to be considered unphysical. 
We will adopt this interpretation here, therefore taking as physical only the m > 
branch of the — c vs m plots. However, the prescription for determining the value of nicr 
continues to be valid, as long as the obtained value of rhcr is positive. Therefore, we will 
continue to use it in order to determine the value ofr] = rjcr for which rhcr = 0. 

As one can see in Figure |4[ the value of t]ct that we obtain is r^cr ~ 7.89. Note also 
that, for this value of r], the Minkowski m = embedding has a non-zero fermionic 
condensate Ccr, and hence the chiral symmetry is spontaneously broken. For r] > rfcr, 
the stable solutions are purely Minkowski embeddings, and the first order phase tran- 
sition disappears; therefore, we have only one class of solutions (the blue curve) that 
exhibit spontaneous chiral symmetry breaking at zero bare quark mass. Some black 
hole embeddings remain meta-stable, but eventually all black hole embeddings become 
unstable for large enough rj. This is confirmed by our study of the meson spectrum 
which we present in later sections of the chapter. The above results can be summarized 




Figure 4: For -q = rjcr the critical parameter rricr vanishes. There are two m = states 
with equal energies, one of them has non-vanishing condensate — Ccr ~ —4.60 and therefore 
spontaneously breaks the chiral symmetry. 

in a single two dimensional phase diagram which we present in Figure |5] The curve 



59 



m 



1 



0.2 



0.8 



0.6 



0.4 




spontaneously broken CS 



2 



4 



6 



8 



10 



12 



Figure 5: The curve separates the two phases corresponding to discrete meson spectrum (light 
mesons) and continuous meson spectrum (melted mesons). 

separates the two phases corresponding to a discrete meson spectrum (light mesons) 
and a continuous meson spectrum (melted mesons) respectively. The crossing of the 
curve is associated with the first order phase transition corresponding to the melting of 
the mesons. If we cross the curve along the vertical axis, we have the phase transition 
described in refs. iTTl fTSlfTTOll . Crossing the curve along the horizontal axis corresponds 
to a transition from unbroken to spontaneously broken chiral symmetry, meaning the 
parameter c jumps from zero to Ccr ~ 4.60, resulting in non-zero quark condensate of 
the ground state. It is interesting to explore the dependence of the fermionic condensate 
at zero bare quark mass on the magnetic field. We can fix this from a simple dimensional 
analysis: 



In the T — )■ limit, we should recover the result from ref. [|53l: Ccr ~ 0.226R^H^/^, 
which implies that Cct{i]) ^ 0.226//^/^ for 77 ^ 1. This rather rich phase structure is 
accompanied by an equally rich thermodynamics. We refer the reader to the original 



Ccr = b^Ccriv) 



^3/2 



(83) 
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paper. ^ for an extensive account of the associated free energy, entropy and magneti- 
zation. 



3.5 Meson spectrum 

In this section, we calculate the meson spectrum of the gauge theory. The mesons we are 
considering are formed from quark-antiquark pairs, so the relevant objects to consider 
are 7-7 strings. In our supergravity description, these strings are described by fluctua- 
tions (to second order in a') of the probe branes' action about the classical embeddings 
we found in the previous sections ||98l . Studying the meson spectrum serves two pur- 
poses. First, tachyons in the meson spectrum from fluctuations of the classical embed- 
dings indicate the instability of the embedding. Second, a massless meson satisfying 
a Gell-Mann-Oakes-Renner (GMOR) relation ll58ll will confirm that spontaneous chiral 
symmetry breaking has occurred. 



3.5.1 The scalar fluctuation 



In order to solve for the meson spectrum given by equation (A. 5 ), we consider an ansatz 
for the field x of the form: 

X = h{z) exp {—iut) , with z = b~'^z , co = R~% Cj . (84) 

In these coordinates, the event horizon is located at z = 1. 

We have two different kinds of embeddings here, but we will consider only the 
massless black hole embedding. Our goal will be to demonstrate that the external mag- 
netic field destabilizes this classical embedding and therefore gives rise to the chiral 
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symmetry breaking phenomena observed in earlier sections. For the black hole embed- 
dings, the appropriate fluctuation modes are the in-falling modes [1731 . which yields the 
corresponding complex-valued quasinormal mode spectrum for a); the real part of a) 
corresponds to the mass of the meson before it melts, and the imaginary part of u is the 
inverse lifetime (to a factor of 2) [|73l . 



The equation of motion (A. 5) for the fluctuation corresponding to the massless 
embedding, 9o{z) = 0, is given by : 



h"i2) + 



2z 



1 



1 ^ (1 + Z^T]"^ 



h'(z) H ^ T^hh) 



0. (85) 



Upon using the change of variables outlined in appendix C, we can obtain the effective 
Schrodinger potential denoted by Vs{z). In this case, the exact form of this potential 
is not very illuminating, so we do not present it here. Alternatively, in figure |6] we 
pictorially demonstrate the effect of the magnetic field on this potential. 

Vsil) 



\ J1 = 



;, = 6 



0.9 



Figure 6: The effective Schrodinger potential for the scalar field fluctuation at two different 
values of the magnetic field. 



It is clear that a sufficiently large non-zero value of r] generates a locally negative 
potential well. This in turn suggests that the imaginary part of tu may become positive 
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signaling the existence of tachyonic modes. Such tachyonic modes are purely sourced 
by the magnetic field which tends to disfavor the black hole embeddings. We will not 
explore the quasinormal modes in any more details, but refer to ref. @ where more 
extensive numerical studies have been performed. 



3.5.2 The vector fluctuation 



In general the scalar fluctuation $ and the vector fluctuation Ai are coupled (as shown 



in equations (A.6) and (A.7)). However, for the massless embedding these coupled 
equations reduce to a single equation 



(2;^ — 1) (1 + Z'^f]^) 
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where we consider an ansatz (as before) of the form: 



Ai = A{z) exp {—iut) . 



Now we can again use the prescription in appendix C to obtain the effective Schrodinger 
potential for the gauge field fluctuation. Again, instead of explicitly writing down the 
potential we pictorially demonstrate the effect of the magnetic field in figure |7] 

The potential has remarkable similarity to the potential of the scalar field fluctuation 
in figure [6] and we again see the local development of a negative energy potential well 
for sufficiently large magnetic fields. This further conforms to the magnetic field's role 
in destabilizing the black hole embeddings. More extensive numerical study has been 
carried out in ref. [|9l. 

It is possible to numerically extract the meson spectrum for the Minkowski embed- 
dings (when the scalar field $ and the vector field A are coupled) and obtain the analogue 
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Figure 7: The effective Schrodinger potential for the gauge field fluctuation at two different 
values of the magnetic field. 

of Zeeman splitting and the Gell-Mann-Oakes-Renner relation. We will not elaborate 
on these aspects here, but point the interested reader to the relevant reference [|9|. 

3.6 Concluding remarks 

We have extended the holographic study of large gauge theory in an external mag- 
netic field described in ref. [53], to the case of finite temperature, allowing us to study 
the properties of the quark dynamics when the theory is in the deconfined plasma phase. 

The meson melting phase transition exists only below a critical value of the applied 
field. This is the critical value above which spontaneous chiral symmetry breaking is 
triggered (in the case of zero mass). Above this value, regardless of the quark mass 
(or for fixed quark mass, regardless of the temperature) the system remains in a phase 
with a discrete spectrum of stable masses. Evidently, for these values of the field, it 
is magnetically favorable for the quarks and anti-quarks to bind together, reducing the 
degrees of freedom of the system , as can be seen from our computation of the entropy. 
Meanwhile, the magnetization is greater in this un-melted phase. 
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There have been non-perturbative studies of fermionic models in background mag- 
netic field as briefly reviewed in chapter 2. Generally, those works use quite differ- 
ent methods to examine aspects of the physics — some primary non-perturbative tools 
are the Dyson-Schwinger equations in various truncations). Our results (and the zero 
temperature result obtained with these methods in the zero temperature case lISBl ) are 
consistent with the general expectations from those works, which is that strong mag- 
netic fields are generically expected to be a catalyst for spontaneous chiral symmetry 
breaking in a wide class of models (see e.g., refs. [[681I1421I113II for a discussion of the 
conjectured universality of this result). 
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Chapter 4: Phase structure of finite 



temperature large Nc flavoured 
Yang-Mills theory in an external 
electric field 

4.1 Introductory remarks 

The material presented in this chapter is based on the work in collaboration with 
Tameem Albash, Veselin Filev and Clifford Johnson[l()|. 

In this chapter, we study large SU(Ai'c) gauge theory with non-backreacting 
hypermultiplet quark flavours in the presence of a background electric field, at both 
zero and finite temperature. The holographic description of this theory at finite baryon 
chemical potential was first studied in ref. [86J, where the existence of a global elec- 
tric current induced was demonstrated. Further study of this set up was considered in 
ref. [I117II . The phase structure of the theory at zero baryon chemical potential was 
studied in ref. [fTOll. 

We may expect the behavior of the theory to be quite different from the magnetic 
case considered in the previous chapter. In a sense the electric field is pulling the quarks 
apart, decreasing the bound energy and thus we may expect that the onset of the decon- 
finement phase transition will happen at lower temperature, than in the absence of an 
electric field. Furthermore, at zero temperature we may expect that for sufficiently 
strong electric field the binding energy of the quarks will be completely overcome and 
the mesons will dissociate into their constituent quarks. 
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It is worth noting that the phase transition will be driven by the quantum fluctuations, 
which will result in a quantum phase transition. On the other side, the mesons are 
electrically neutral, while the dissociated quarks aren't. Hence the dissociation can be 
seen as a conductor/insulator phase transition [fTOl . 



4.2 General set-up 

We consider the form for the metric of the AdSs-SchwarzschildxS*^ background pre- 



sented in equation (59 1: 



+R^d9'^ + sin^ Odcj)'^ , (86) 



ds'^ = -j-^ ( -'-^dt^ + f2dx'^ ) + —dr'^ + cos^ Od^l 



In order to introduce our background electric field (and current), we consider an 
ansatz for the world- volume gauge field of the form ll86l : 



Ai{r) = -Et + F{r) . (87) 

This ensures a constant electric field Fqi = E m the gauge theory; F{r) sources a 
current on the world- volume [|86ll . For the D7-brane embedding, we consider an 
ansatz as before: = 6{r). 

The DBI action is then given by 

3, 



^D7 = -ToriV; / d^i m (1 + rH\rf) 



-4 {2Tia'f {-flf2F'{rf + ^flR'E^ (l + r''e\rf) ] ^'^ 



■ (88) 
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The equations of motion derived from equation (88 1 for F(r) introduces a constant 
of motion K. We can invert the equation of motion to write a first order differential 
equation for F{r) in terms of K: 



{2naf cos' Off f,F' 



K 



In the limit of r — )■ oo, the asymptotic solution for F{r) is given by: 



(89) 



lim F{r) 



Normalizability of the solution requires us to take /o = 0. The constant K is related to 
the (vacuum expectation value) vev of the current by the following relation [|86ll : 



(Ji) = (V;7V) = -^-n^a'VNfTjijK 



(90) 



where V is the spatial volume of the gauge theory. It is convenient to exchange the field 
F{r) in favour of the constant K by performing a Legendre transformation: 



'D7 



S^r-ld'^F'ir)§ 



r) 



-2Ti^VNfTjij I dt I dr 



16r5/i\/72 



(/2 - 16r4 {2TTa'f E^) {-6Ar^K^fl + fffi cos^ e{r)) 



(91) 



It is a simple check to show that SIbj/SK = F'{r). The expression under the square 



root in equation (91 1 must be positive in order to keep the action real. Either both terms 
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should be negative or both should be positive, and they must change sign at the same 
radial distance r^:. This results in two conditions: 



[Art - b^Y-^Grt {27caf = 



{2'KaE) + J {2Txa'EY + h'^ 



(Art - b^f i'irt + h^f cos^ e{n) 
64r6 (4r4 + b^f 



(92) 
(93) 



Note that equation (93 1 relates the current to the electric field in a form of Ohm's law, 
and this was used in ref. ri86i to determine the conductivity of the fundamental matter in 
the quark-gluon plasma. We refer to r^, as the "vanishing locus". From the action Idj 



defined in equation (91 1, we can derive the equation of motion for the embeddings 9(r). 
Its exact form is not illuminating, and so we do not display it here. In the limit of large 
r, the solution of the equation of motion asymptotes to the familiar form: 



e{r) 



m c 

7 + ^:3 



(94) 



The constants m and c are related to the bare quark mass and the condensate respectively, 
in a manner that was reviewed in Chapter 3. 



4.3 Properties of the solutions 
4.3.1 Exact results at large mass 

It is instructive to study the properties of the quark condensate as a function of the 
bare quark mass for large mass. This corresponds to m ^ (27ra'E) /bin terms of 
the dimensionfuU quantities. This limit constrains us to consider only the so-called 
"Minkowski" embeddings for the probe D7-brane for which the vev of the current 
(defined in equation ([90])) vanishes by virtue of equation (93 ). To extract the behaviour 
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of quark mass and condensate we linearize the equation of motion obtained from equa- 



tion (91) in the same way as described in ref. [9J. Rewriting the result in terms of 



dimensionful parameters, we find the following analytic behaviour for the condensate: 

- R^E^ ^ + Ab^R^E^ + SR'^E^ ^ ^ ^^^^ 



where E = 27ia'E. This may suggest that for high enough bare quark mass, the con- 
densate vanishes. However unlike the results for an external magnetic field [HI, the 
condensate approaches zero on the positive side of the condensate axis. We verify these 
results in the next section using numerical techniques. 



4.3.2 The case of vanishing temperature 

It is instructive to first study the zero temperature case, since a number of important 
features appear in this case. For this simple case, the energy scale is set by Rv^, and 
therefore it is convenient to introduce the dimensionless quantities f and m via: 

r = rVe f ; m = rVe rh . (96) 



We solve the equation of motion for 9{r) using a shooting technique by imposing 
infrared boundary conditions described in refs. [|71[8l[85l. For "Minkowski" embed- 
dings, in order to avoid a conical singularity at 6^ = n/2 [|85l , we impose the following 
boundary condition: 

^'(f)|,=^/2 = -oo. (97) 

For embeddings reaching the pseudo-horizon we find it numerically convenient to shoot 
forward (towards infinity) and backward (towards the origin) from the vanishing locus 
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f* = 1. The boundary condition ensuring the smoothness of solutions across the van- 
ishing locus can be determined from the equation of motion itself, and it is found to 
be: 




^0 = 0{n) , 



- tan 
du u 2 



de 1 e. 



(98) 



From the above boundary condition we can see that the embeddings reaching the van- 
ishing locus at = 7r/2 avoids the conical singularity by having a diverging derivative. 
Next we proceed to discuss our numerical findings. 

Numerical analysis shows that there are three different classes of embeddings, which 
can be classified by the topology of the probe D7-brane. First, we have the smooth 
"Minkowski" embeddings that close above the vanishing locus, due to the shrinking 
of the that the probe brane wraps. Second, we have the embeddings that reach the 
vanishing locus before the shrinks. These embeddings can in turn be classified into 
two different categories. The first are singular solutions that cross the vanishing locus 
and close before reaching the origin. These have a conical singularity at the closing point 
because the derivative does not diverge. We discuss these solutions further below. The 
second are smooth solutions that pass through the vanishing locus and reach the origin 
with no singular behaviour. These different embeddings are summarized in figures [8] 



We can explicitly see the conical singularity appearing for the brane (second curve in 



the magnitude of the electric field this class of singular solutions starts appearing after a 
certain value of (6'o)mm' close to 7r/2 (which corresponds to the critical embedding), is 
reached and persists until 6q = n/2. 



andUl 



Figure |9]) 



since it intersects the vertical axis at some non-zero value of 6. Depending on 
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rsin^ 




Figure 8: The solid curves starting far left (red) represent solutions smoothly closing before 
reaching the vanishing locus. The solid remaining curves (blue) correspond to embeddings pass- 
ing the vanishing locus (denoted by a semi-circular dashed black curve). 

rsind 
1.4 




Figure 9: Two kinds of embeddings that pass the vanishing locus corresponding to D7-brane 
closing at the origin (left), and forming a conical singularity (right) respectively. The semi- 
circular (red) curve corresponds to the vanishing locus. 

The existence of these embeddings suggests the possibility of a transition in topology 
of the probe brane as a function of the parameters, as happened for finite temperature. 
Here, it can happen as a result of the external electric field that we've applied. The 
Minkowski embeddings simply have a shrinking S^, while the embeddings reaching the 
origin are distinguished by having an shrinking as well as touching the AdS horizon. 
However, the presence of singular solutions makes the nature of the transition subtle, 
as we discuss later. We will later see that finite temperature case also reveals similar 
classification of embeddings. 
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From the asymptotic behavior of these embeddings we can extract condensate as a 



function of the bare quark mass. In Figure 10 we show this dependence. There are two 
important mass scales. One, rhcr, is where the phase transition between the two types 
of embedding occurs, and the other m*, is where the singular solutions first appear. 
These values are in turn determined by the parameters (^^o)cr and (^^o)min, respectively. If 
^Omin > (^o)cr, thcu the siugular solutions are never thermodynamic ally favored. We see 



m„=1.328 




m„= 1.328 




1.315 1.320 1.325 1.330 1.335 1.340 



Figure 10: The condensate as a function of bare quark mass rh. The curve segment coming 
in from the left (blue) corresponds to embeddings reaching the vanishing locus and the curve 
segments going out to the right (red) correspond to Minkowski embeddings. There is a family 
of solutions that contain a conical singularity between the vanishing locus and the origin. They 
exist for (^o)min < ^0 < (^o)cr- On the right is a magnification of the turn-around region where 
these segments join, showing multi-valuedness. An analysis of the free energy reveals a first 
order phase transition at rhcr where there is a jump from one type of embedding to another. 

the presence of a first order phase transition speculated based on the general arguments. 
Furthermore, we observe that there is no chiral symmetry breaking, since in all cases 
the condensate vanishes for vanishing mass. This is to be contrasted with the case of 
external magnetic field studied with these methods in the previous chapters. 

This latter observation fits our intuition that the mesons in this theory have a binding 
energy that grows with (it is proportional to) the constituent quark mass [|98l . For a 
given quark mass there should exist sufficiently high electric field that can reduce the 
binding energy the quarks. This allows two things to happen: First, this inhibits the 
formation of chiral condensate, and second, this ultimately will dissociate the mesons 
into its constituent quarks. This dissociation is in fact a transition from an insulating to 
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a conducting phase, mediated by the external electric field. On the dual gravity side this 
corresponds to a transition from Minkowski embeddings that do not reach the origin to 
embeddings that do. 

We must note that the appearance of the singular solutions (those that have a conical 
singularity) are not well understood at the moment. There is therefore the possibility that 
there is an as yet to be identified intermediate phase right after the meson dissociates. 
One simple possibility (but not the only one) is that stringy corrections smooth out the 
singularity in the interior, while preserving the asymptotic behavior. 

The finite temperature version of the story rather similar to the vanishing temperature 
case. We will not elaborate on this case further, but refer the reader to the original 
paper. M^. 

4.4 Concluding remarks 

It is very encouraging that non-trivial non-perturbative phenomena resulting from exter- 
nal fields such as those we have seen here (a dissociation phase transition, metal- 
insulator transition and the associated response current) can be so readily extracted in 
this kind of holographic study. We found that since the electric field works together 
with the presence of finite temperature, the resulting phase diagram which accounts for 
the effects of both is rather simple (ignoring the complication of the special conical 
solutions — see below). 

Clearly, the story is not quite complete in this electric case, since above a certain 
value of the electric field, some of our solutions develop conical singularities in the 
interior. We are not entirely sure about the nature of these solutions. As already stated, 
one possibility is that the solutions are locally (near the singularity) corrected by stringy 
physics, perhaps smoothing the conical points into throats that connect to the horizon. 
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We would expect in that case that our phase diagram would largely remain intact, since 
the values of the masses and condensates for each solution are read off at infinity, far 
from where the conical singularity develops. The results for the relative free energies of 
those solution (compared to the other solutions at the same values of the mass) would 
be the same, and so the complete story would be unaffected. A more drastic possibility 
is simply that the physics is considerably modified by instabilities associated with those 
conical solutions that our current study has not revealed — the appropriate part of the 
meson spectrum is particularly difficult to study numerically here, and so the tell-tale 
signs of tachyonic modes that represent an instability are hard to check for in this case. 
Such a modification could remove quite a significant part of the phase diagram for large 
enough values of the electric field, allowing the possibility of a richer structure than we 
have seen so far. 
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Chapter 5: External Fields and Chiral 



Symmetry Breaking in the 
Witten-Sakai-Sugimoto Model 

5.1 Introductory remarks 

We will change gears to discuss the dynamics of flavours in large Nc gauge theories in 
the holographic model in Type IIA framework. We have reviewed this set-up briefly 
in the introduction, but we will recollect some key features and build up from there. 
We will be interested in very similar physics as we have discussed in the previous two 
chapters. The material presented in this chapter is based on the work in collaboration 
with Clifford JohnsoniH. 

The Witten-Sakai-Sugimoto model, as described in ref. HI 3711 is one construction 
which cleanly realizes chiral symmetry breaking and confinement/deconfinement tran- 
sition. The supergravity background of this model is constructed of near-horizon geom- 
etry of Nc D4-branes, following ref. HI 591 . The study of Nj flavour D8 branes in this 
background when Nf <^ Nc reveals a nice geometric realization of chiral symmetry 
breaking. The flavour branes do not backreact on the background geometry in this 
probe limit and therefore studying their dynamics using the Dirac-Bom-Infeld (DBI) 
action (including a Wess-Zumino term, if necessary) suffices to capture the correspond- 
ing gauge theory dynamics of fundamental flavours in an analogue of the quenched 
approximation. 

In this (DBI) regime it is possible to capture general gauge theory features such as the 
phase diagram for temperature vs chemical potential by considering probe brane in finite 
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temperature supergravity background and exciting specific gauge field on the world- 
volume of the probe brane. Previous such studies including the non-zero chemical 
potential in this model have been carried out in e.g., refs. Il72l 11361 |22l |38l. Here, we 
will introduce an external magnetic and electric field. 

A clear method for introducing a background magnetic field has been previously 
discussed in the chapters 3 and 4. We considered pure gauge fi-field in the supergravity 
background, which is equivalent to exciting a gauge field on the world-volume of the 
flavour branes corresponding to a magnetic field. 

We find that the presence of magnetic field promotes the spontaneous breaking of 
chiral symmetry as found previously in chapter 3 in a different model. This is expected 
from the field theory perspective and is widely recognized as a sort of "magnetic catal- 
ysis" for chiral symmetry breaking (see e.g., ref. [111311 ). Further study of the phase 
structure of this model reveals the existence of a finite critical temperature (for restora- 
tion of chiral symmetry) for large magnetic field. We analyze a number of other physical 
quantities such as the latent heat and relative magnetisation associated to the phase tran- 
sition. 

We also find that our phase structure is rather generic for the Witten-Sakai-Sugimoto 
type holographic models where we consider the dynamics of probe Dp-brane in D4- 
brane background. We address some of the physics of an external electric field in the 
model. We find that in the symmetry-restored phase an external electric field drives a 
current in the gauge theory due to pair creation, and the symmetry-broken phase does 
not conduct. However we have not considered the presence of baryons in our set-up, 
which could give rise to a non-zero current in the phase where chiral symmetry is bro- 
ken. 
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5.2 The Witten-Sakai-Sugimoto Construction 

The Witten-Sakai-Sugimoto model H 13711 consists of near-horizon limit of Nc D4-branes 
wrapped on a circle of radius R in the x'^ direction with anti-periodic boundary condition 
for fermions. The D4-branes are intersected in the compact direction by Nf D8- 
branes at x^ = — ^ and Nf D8-branes at x^ = ^ (with the constraint that L < txR). 
This is dual to a (4 + 1) -dimensional SU(Ai'c) Yang-Mills theory with gauge coupling 
constant \ left and right handed quarks are introduced by the flavour D8 and D8-branes 
in the probe limit that share three spatial directions with the D4-branes. The flavour 
branes introduce a global flavour symmetry U(A^/)l x U(A^/)/j as seen from the (4 + 1)- 
dimensional D4-brane worldvolume gauge theory. This global symmetry is identified 
with the chiral symmetry (non-abelian) of the effective (3 + l)-dimensional gauge theory 
with chiral fermions. In the probe limit, the dynamics of the flavour branes is described 
by the DBI action in the background of Nc D4-brane geometry. The background metric 
of D4-brane is obtained from Type IIA supergravity and is given by 

ds^ = (1)'^' {-dt' + dx.dx^ + f{u){dx'f) + (1)"'^' (^-^ + u'dni^ , 

e =9s[^^) , i^(4) = -j^^Q , /(m) = 1-1 — 1 , R-' = 7rgsNJi. 

(99) 

Here are the flat 3-directions, t is the time coordinate, x^ is the spatial compact circle, 
^4 are the 5*^ directions and u is the radial direction, ig is the string length, gg is the 
string coupling; V4 and £4 are the volume and volume form of S'^ respectively. Also, $ 
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is the dilaton and F(4) is the RR four-form field strength. To avoid a conical singularity 
in the {x^, u} plane one should make periodic identification: 



47r / 



1/2 



27ri?4 . 



(100) 



This endows the background with a smooth cigar geometry in the {x^, u} plane. 
On the other hand, the finite temperature Euclidean background is given by 



RJ 



[dxidx' + f{u)dtl + {dx^f) 



U \ -3/2 / du 



R 



+ u^dVtl 



tE — tE 



47r/?3/2 



1/2 



T 



T 



1 _ UttR'/A 



fin) 



U' 



T 



U 



• (101) 



All the parameters are given by the same formula as equation (99). The dilaton, RR 



4-form, the radius R are also given by the same formula as equation (99 1. 

It is clear that at finite temperature we have two candidate supergravity solutions: the 
Euclidean version of ( 99 1 and the background in ( 101[ ). Depending on the temperature of 
the background one of these two possible backgrounds is energetically favoured. At low 



temperatures, <!/ {2ttR4), the Euclidean version of (99 ) is the favoured background 



whereas for high temperatures, > 1/ (27ri?4) the background in ( 101 1 is energetically 



favourable. At Td = 1/ (27ri?4) there is a first order phase transition, which corresponds 
to the confinement/deconfinement transition in the dual gauge theory. Henceforth we 
will sometimes refer (99) as the confining background and ( 101| ) as the deconfining 
background. 

Now one can introduce the flavour brane-anti-brane system in the probe limit, 
namely Nf <^ N^. In this limit the probe branes do not backreact on the geometry and 
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the classical profile of the probe is solely determined by the Dirac-Born-Infeld action. 
We will consider the following ansatz for the flavour D8-D8 branes: 

{t, Xj, x*^ = r, ^4, M = u(r)} , (102) 

where we note that the coordinates in the parenthesis should be understood as the world- 
volume coordinates of the D8/D8-brane. Before discussing the effect of the external 




Figure 1 1 : The possible shape of embeddings of the probe brane. The left diagram shows the 
only possible (U)-shape of the probe brane in the confining background. On the right, we have 
shown two candidate embeddings of the probe brane in the deconfined phase: the U-shaped (in 
red) and the straight (in blue). 

magnetic field, let us review the physics when temperature is the only external parame- 
ter. We can understand the main qualitative feature without resorting to actual calcula- 
tion. 
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The profile of the probe brane is clearly a function denoted by x'^{u) or m(x^). In 
the low temperature regime, the relevant background has a smooth cigar-shape in the 
{x^, u}-plane. Therefore the D8/D8-pair separated by a distance L at the boundary is 
bound to join together in the bulk. This is illustrated in the left of figure[TT] This implies 
that the chiral symmetry group U(A^/)l x U(A^/) r is broken down to a U(A^/)diag, which 
is a rather simple geometric realization of chiral symmetry breaking. Thus in this model, 
in the low temperature regime chiral symmetry is dynamically broken. 

The physics of chiral symmetry breaking is different in the high temperature regime. 
In this case, the relevant background does not have a cigar shape in {x"^, m} -plane (this is 
illustrated in the right of figure [TT]). Therefore for a given asymptotic brane-anti-brane 
separation L, it is possible to have two candidate solutions. The brane-anti-brane pair 
can join smoothly as before which realizes the spontaneous breaking of chiral symme- 
try. We will refer such embeddings as the "curved" or the "U-shaped" embeddings. The 
brane-anti-brane pair can also end separately on the bulk horizon which corresponds 
to restoration of the full asymptotic chiral symmetry group. We will refer such embed- 
dings as the "straight" embeddings. For a fixed background temperature, depending 
on the asymptotic separation L one of these two possible embeddings are energetically 
favourable. It can be shown that for "low" temperatures the "curved" embeddings are 
the favoured one, which means that the chiral symmetry is dynamically broken. At 
"high" temperature, the "straight" embeddings are favoured and therefore the broken 
symmetry is restored. This is a first order phase transition. We will make precise sense 
of what we mean by "low" and "high" temperatures next. 

We have two first order phase transitions taking place in this set-up. The background 
undergoes a confinement/deconfinement transition and the probe flavours undergo a chi- 
ral symmetry restoring transition. The first transition is controlled by the radius R4 of the 
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compact x^-direction and the second transition is controlled by varying L, the asymp- 
totic separation of the brane-anti-brane pair. Therefore the only relevant dimensionless 
quantity is the ratio L/R^. 

If the separation L obeys the bound that L/R4 < 0.97, then there exists a phase 
which is deconfined with broken chiral symmetry. We will sometimes refer this as the 
"intermediate" phase. In this case, the chiral symmetry restoring temperature is given 
by ^xSB — 0.154/L. Thus in this phase there is a separation of scales between the gluon 
deconfinement temperature and the chiral symmetry restoring temperature. If however 
L/i?4 > 0.97, then deconfinement and chiral symmetry restoration happens together. 
This has been discussed in refs. [|6l ll20L 

We can also set i?4 — )■ 00 so that the spatial circle direction now becomes a flat 
extended direction; this particular limit has been studied in detail in refs. [[TTlfTllfTSll . It 
turns out that in this particular limit the dual gauge theory is described by the Nambu- 
Jona-Lasinio model (when probing with D8-branes) or the Gross-Neveu model (when 
probing with D6-branes) with a non-local four-fermi interaction. The strength of the 
four-fermi interaction is set by the asymptotic separation L. In this case, the confining 
phase does not exist and therefore the chiral symmetry restoring transition is the only 
phase transition. 

5.3 The probe brane analysis 

Many aspects of the finite temperature physics have been studied before, e.g. in 
refs. [|6l, [|120|| and references therein. To introduce external magnetic field we excite 
a gauge field A3 = Hx^ on the worldvolume of the probe brane. 
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5.3.1 The low temperature background 

For notational convenience we rename coordinate to be r. In this case the relevant 
background is given by equation (99). It is straightforward to check that in this case 
the Wess-Zumino term does not contribute to the probe brane action, which is therefore 
completely determined by the DBI action: 

Sds = /is J d^^e-^^Jd^i{P[G^, + B^,]) = C J drC{u,u') , (103) 

where C = n^Vs^V^s / QgT and u' = du/dr. Recall that we put the flavour branes 
(D8/D8) with the asymptotic condition that as u — t- oo, r — )• ±L/2; where the ± 
corresponds to the D8 and D8-brane respectively. Also note that the coordinate r is 
restricted from —nR to +itR. 



One can immediately note from equation (103) that C{u,u') is independent of r, 
therefore the Hamiltonian corresponding to r will be a constant of motion. Carrying out 
the following Legendre transformation we should have 

ou' 



So the first integral of motion that follows from equation ( 104) is given by 



u' 



(l + H^SfY f(u) I / B\'\'' 

\ ^'''> '[! ^udi + H-'i^) VTM. (105) 



We have rewritten the constant in the right hand side in a convenient way (in terms 
of a new parameter Uq). Note that Uq is the minimum value of u that the probe brane can 
reach satisfying u'\u=Uo — 0- We assume that the brane-anti-brane pair join smoothly, 
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which implies that there is no resultant force present at the point where they meet. Typ- 
ically this would mean that there is no other source (e.g., a baryon vertex or a bunch of 
F-strings) present at this point. For zero background magnetic field this set up reduces 
to the low temperature case analyzed in ref. Let us now focus on the solution for the 
probe brane profile. 

We will compare the behaviour of the brane profile in the presence of magnetic field 
to the case when it is turned off. For notational convenience we define the following: 

y = TT ' ^KK = , R = Uod, L = UqI . 

With the above redefinitions we can obtain the difference in slope of the profile in pres- 
ence and in absence of magnetic field as 

with y E [1, oo]. So we get that < |m'//=oI e^ch value of y. This in turn means 
that the magnetic field bends the profile of the D8/D8 brane and therefore forces the 
brane-anti-brane pair to join closer to the boundary (and hence break chiral symmetry) 
for fixed asymptotic separation. 

We can study this explicitly as follows. The brane-anti-brane separation at the 
boundary (m — )■ oo) is given by 



2 J JuoiH) ^'h \/Uq(h) 



1/2 

\ -i- HI \ — \ r / \ 

f{y) 



/oo 
IiH)iy)dy. (106) 
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Clearly putting if = we get the corresponding separation when the background mag- 
netic field is switched off. 



L 
2 



^3/2 



y ^^"^dy 



Ji 



m {j^)y' - 1) 



1/2 



dyXj 



H=0 



(107) 



For the same asymptotic separation magnetic field changes the brane profile's point 



of closest approach Uo(^h)- We can compare f/o(//) and Uq. Equating equation ( 106) and 
equation ( |107| ) one gets 

lH=ody 



0{H) 



Uo J^X^H)dy 



(108) 



Some algebra shows that Ih=o > ^{h) for all y, so the ratio on the right hand side 
is greater than or equal to one, which also means that Uq(^h) > Uq. Therefore for the 
same asymptotic separation the magnetic field can only help to join the brane-anti-brane 
pair favouring chiral symmetry breaking. This is pictorially represented in figure [12] 
We can extract more qualitative features in appropriate limits. To do so, let us rewrite 



u 




Ukk 
T 



Figure 12: The dashed U-shaped curve represents a profile in vanishing background field and 
the solid U-shaped curve represents a profile when a non-zero magnetic field is present. These 
profiles are obtained by numerically solving the equation of motion for the probe brane. 
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equation ( 106 1 with the change to variable z = y ^. With this equation ( 106 ) becomes 



'o Jo 



z 



2 svl^v 7„ /'-'-nl'r""' -(i-i/?<K)(i + g-^») 

(109) 

Now small asymptotic separation corresponds to large values of Uq which means 
Ukk ^ 1. So for small L and weak magnetic field (1 / cP^"^ ^ H), the leading behaviour 



of the separation is given by (using equation (109)), L ~ -R^/^/i/fTo- This is same 



as the leading behaviour in zero magnetic field case in ref. [6J. However, for strong 



magnetic field (l/rf^/^ <^ H), the leading behaviour obtained from equation (109) is 
given by, L ~ R'^H/Uq. So for fixed value of Uq the asymptotic separation scales with 
the applied magnetic field strength H. This is however true only in the ^/kk ^ 1 limit. 

Note that in this model the bare quark mass is always zero as there is no separation 
between the flavour and the colour branes at the boundary; however since the branes 
join at some length scale Uq > f/xK in the core one can consider a string stretching 
from u = Ukk io u = Uq. The mass associated with the string can be identified to 
be the effective constituent quark mass as argued in ref. [6]. If we denote this mass by 



Mq, then Mq = J^° y/gttguu , where Uq has to be determined from equation ( 109 ) 



for given L. This turns out to be a self-consistency equation for Uq. This is analogous 
to the Gap equation in the field theory context (e.g., in ref. [llj). The self-consistency 
equation turns out to be 

^0 = ^f^ (1 - ^kk) {U^ + H'R') liUo, Hf , where 



I{Uo,H) 



9L2 



.3 .^-1.-5/6^, 



(1 - yl^z) (f/3 + Hm^) - (1 - Vl^) + iJ2^3))l/2 " 

(110) 
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Now equation (|1 10|) can be solved perturbatively, i.e., starting with an initial value 



for the parameter Uq we can determine the next order approximation to Uq using equa- 



tion (110); and we continue until the desired accuracy has been achieved. It is straight- 



forward to guess the initial value of f/o. Plugging in if = in equation ( 1 10 1 we should 
get the constituent mass for the low temperature case. This can serve as the initial guess 
for small magnetic fields. Once Uq is known for small magnetic fields, it can be used 
as the initial guess for successively higher values of magnetic fields. Thus we obtain 
the dependence which is shown in figure below. It is expected that the mass of the 
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Figure 13: The dependence of constituent quark mass (measured in units of (27ra') ^) on the 
external magnetic field in the low temperature phase. We have fixed Ukk = 0.4 and R = 1. 



vector and pseudoscalar meson would monotonically increase as the constituent quark 



mass increases. Figure 13 therefore should capture the behaviour of massive vector and 
pseudoscalar meson spectra in presence of magnetic field. 

Note that for this background disjoint brane pair do not exist: The constant r- 
solutions namely, r = — 7ri?4/2 and r = +7ri?4/2 join at u = Ukk- This is because of 
the cigar geometry of the background in the {t,u} submanifold and the fact that r-circle 
is wrapped by the probe branes. The brane pair must join together since there is no place 
in the geometry for them to end separately. So the only configuration possible in the low 
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temperature phase breaks chiral symmetry by reducing the global U(A^/)i x U(A^/)r 
to the diagonal U(A^/). This is geometrically understood as the joining of the flavour 
8-branes. 



5.3.2 The high temperature background 



Recall that the high temperature background is given by equation ( 101[ ); and we again 



use the same ansatz for the probe given by equation ( 102| ). As before the profile of the 



probe brane is completely determined by the DBI action, from which a first integral of 
motion can be readily obtained to be 




For convenience we use the dimensionless variables defined in equation (106) along 
with the new variable yx = Ut/Uq. Starting from the first integral of motion in equa- 
tion ( |1 1 1| ) it is easy to verify that the finite temperature analogue to equation ( 106 ) takes 
the following form 



/2 /2 



m^/M!,s H^d^ (112) 
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This also suggests that \u'jj\ < |m'^^q|, leading us to the same conclusion that the mag- 



netic field helps bending the branes. The analogue to equation ( 106 1 and ( 107 1 now take 
the following forms 



L _ i?3/2 
2 ^/Uo(H) J I 



^3/2 



y ^^'^dy 



f{y) 



y> „,8 



2 ^/iMmJi ^ 



/(I) 1+^2^3 y J- 



/M„8 _ 1 



1/2 



1/2 



(113) 



leading us to a similar conclusion as the low temperature case. This is pictorially repre- 



sented in figure 14 As before one can extract the dependence of the asymptotic separa- 




Ut 



Figure 14: The dashed U-shaped curve represents a profile in vanishing background field and 
the solid U-shaped curve represents a profile when a non-zero magnetic field is present. The 
straight (red) solution does not have any qualitative change in presence or absence of the external 
field. These profiles are obtained by numerically solving the equation of motion for the probe 
brane. 



tion in the small yr limit. Again with the change to variable z = y we get 

L H'/^ r — — — {\-y\zY^^'' z-^'lHz 



J{\-y\)(\^H^d?) [ 

'^y^Uo(H) ^ Jo 



(114) 
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So from equation (fTT4j) one can see that for small y-r and weak magnetic field L 



R^^"^ / \/Uo(H)', whereas for small yr and strong magnetic field we get L ~ B?H /Ul^jjy 
The joining of the flavour branes inside the core can be associated with the effective 
constituent quark mass. This corresponds to a self-consistency equation for Uq as in the 
low temperature case. The equation can again be solved using the same perturbative 



approach and the results are summarised in the figure 15 This behaviour of constituent 
mass is valid when the curved solutions are the lowest energy solutions (the chiral sym- 
metry broken phase). In the presence of finite temperature there will be a first order 
transition to chiral symmetry restored phase. We will study this transition later in the 
next section. 




Figure 15: The dependence of constituent quark mass (measured in units of {2'Ka') ^) on the 
external magnetic field in the high temperature phase. We have fixed Ut = 0.3 and R = 1. 



Now the trivial solution of equation ( |1 1 1[ ), which is given by r' = has a different 
physical meaning from the point of chiral symmetry breaking. Since the {t, u] subman- 
ifold has the cigar geometry, the solutions r = ±L/2 can remain disjoint and end at 
u = Ut separately. Therefore the trivial solutions in the high temperature case preserve 
the full \J{Nf)L X \]{Nf)ji symmetry by remaining disjoint. In order to determine the 
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true minimum energy configuration we need to compare the energies of the curved and 
the straight branes. We pursue this in the next section. 



5.3.3 Chiral symmetry breaking 

To determine the true vacuum we consider the difference between the energies of the 
curved and straight branes, which is given by 



1 

- 1 



f dyy {y' + H'd'f' . (115) 



Here AS" < would mean chiral symmetry breaking, AS* > would mean chiral sym- 
metry restoration and AS* = would characterize a transition from symmetry broken 
phase to a symmetry restored one. We employ numerical analysis to study this. 

It is known from, e.g., ref. [6| that for high temperature and zero magnetic field 
there exists a critical temperature beyond which the straight branes are energetically 
favoured, implying that in the dual gauge theory chiral symmetry is restored. Below the 
temperature, however, the symmetry is broken by energetically favoured curved brane 
pair that join together. 

From the zero of AS we can find out the critical value of yx for which the symmetry 
restoration occurs. Now to represent the phase diagram in terms of physical quantities, 
we recall that there is a length scale L corresponding to the separation of the brane-anti- 
brane pair at the boundary. So we express the chiral symmetry restoring temperature 



in units of 1/L using the critical value of yx in equation (114). The resulting phase 
diagram is shown in figure [T6j It is interesting to note that the presence of magnetic 
field increases the symmetry restoring temperature. In other words it promoted the 
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Figure 16: The phase diagram mT — H plane. The vertical axis is expressed in units of 1/L 
and the horizontal axis is expressed in units of R. 



spontaneous breaking of chiral symmetry. This fits with the general expectations from 
field theory (see e.g. refs. HllSID and the supergravity /probe intuition that introducing a 
magnetic field places more energy into the system; therefore in order to minimize the 
energy, condensates are formed (the branes bend more) resulting in more readily broken 
chiral symmetry. (It should be noted, however, that in this specific holographic model 
the identification of a quark condensate is a rather subtle issue (see e.g., ref. l!25l)). An 
interesting and compelling identification of this condensate was made in ref. [4J in terms 
of an open Wilson line (or in short, the OWL) operator. It was further demonstrated in 
ref. [fT4l that this OWL operator does have the right dependence with respect to the 
magnetic field which is consistent with the "magnetic catalysis" phenomenon. 

We can extract some more information about the transition by studying certain ther- 
modynamic quantities at the phase transition. To that end, let us note that the first order 
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phase transition is accompanied by entropy density that jumps aXT = T^. yielding a 
non-zero latent heat as reported in ref. [120], also a change in magnetization 

1 5 (iScurved >S'straight) ^ A 

S — —77 7^ 5 (-'latent — -l-c^S , 

^^~\^^ dH ■ ^^^^^ 

The absolute free energy and any thermodynamic quantity obtained from it (such 
as the absolute magnetization) for the two classes of embeddings (the straight and the 
curved branes respectively) are formally divergent quantities. Hence we compute the 
relative quantities which are finite. We studied the dependence of the change in entropy 
density and the relative magnetization numerically, and the results are shown in figures 



17(a) and 17(b) respectively. The relative magnetization also shows a similar saturation 
behaviour for high enough magnetic field. The straight branes correspond to the melted 
phase where quarks are free whereas the curved branes correspond to the mesonic phase 
where quarks exists in the form of bound states or chiral condensates. Therefore it is 
expected that the chiral symmetry restored phase (corresponding to the straight branes) 
is more ionized than the chiral symmetry broken phase (corresponding to the curved 
branes). This is in accord with our observation that the relative magnetization is negative 



in figure 17(b) 



5.4 A note on background electric field 

For this section we again consider the 8-branes as flavours. We can introduce a 
background electric field by considering the following form of world-volume gauge 
fieldiHISa: 

Ai{t,u) = {-Et + h{u)) . (117) 
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(a) (b) 

Figure 17: The behaviour of jump in entropy and change in magnetisation at the critical tem- 
perature where chiral symmetry is being restored. The vertical axis is evaluated in units of 
and we have set d = 1. 



This means we have a non-zero constant electric field along . The function h{u) 
encodes the response of the fundamental flavours to the external field, namely it encodes 
the information of the non-zero current when flavours are free to move and therefore 
conduct. Note that here we do not introduce any chemical potential, therefore in the 
gauge theory there is no a priori candidate for carrying the charge. However, there could 
still be current caused by pair creation in presence of the electric field. We comment on 
some observations and expectations henceforth. 

For simplicity we assume 2Txa' = 1. We choose the same ansatz for the probe brane 



profile as equation ( 102 1. The DBI action can be computed to be 

9uu = Guu + T , (118) 



where C = {Nffis)/(y^3Vs^), G^i, is the background metric and Quu is the induced 
metric component on the world-volume of the D8-brane along w-direction. 
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Now we will have two constants of motion corresponding to the function t{u) and 
h{u) as follows 



Ge-'^^/d^iG^^Guh' 

^jGuh'^ - {E^ - GuG^,){Guu + GrrT'^) 



B, (119) 
J . (120) 



The constants B and J are related to the minimum radial distance along the world- 
volume of the probe brane (Uq as before) and the gauge theory current (a similar identi- 
fication made in ref. [f86ll holds here also). 



Solving equation (119) numerically to look for all possible solutions is a difficult 



problem. Nevertheless, it can be shown from equation (119) that for confined phase 
there can be only the joined solution (curved as before), for which the current vanishes 
if we impose the condition that the brane-anti-brane pair join smoothly at Uq. This can 



be shown by formally solving equation ( 1 19 ) to obtain t{u) and h{u) as a function of u 



and the constants B and J. If we expand the solution around the joining point u = Uq 
and demand that r' — t- oo as n — t- Uq, we obtain that the current identically vanishes. 
However, if the branes do join but not smoothly, this is no longer necessarily true. One 
could, in this case, have cusp-like solutions for which the brane-anti-brane pair join at 
an angle and in order to stabilize the system it is necessary to consider the inclusion of a 
bunch of fundamental strings extending from the joining point u = U^iou = Ukk- This 
leads us to a construction much like in refs. [|136[ |22l |38l, where the effect of baryons 
was considered. Therefore the symmetry broken phase may have a non-zero current 
carried by the baryons. The chiral symmetry is always broken in this phase forced by 
the topology of the background. 

For the deconfined phase, it can also be shown in a similar way that there exists 
curved solutions joining smoothly at some Uq, which has zero current modulo the caveat 
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mentioned in the last paragraph. This is intuitive from the gauge theory point of view, 
since we are in a chiral symmetry broken phase therefore there is no charge carrier 
present to conduct (ignoring the possibility of a baryon current). The possible effect 
of pair creation is diminished by the existence of quark bound states in the symmetry 
broken phase. However, for straight branes we expect non-zero current to flow. 

A familiar fact from studying flavours in electric field tells us that the presence of 
electric field induces a so called "vanishing locus" for the probe DBI action. The "heal- 
ing" procedure (e.g., in ref. [,53,1 ) is to give a non-zero vev to the current. In practice this 



is obtained by substituting the functions t{u) and h{u) from equation ( 1 19 1 in favour of 



the constants B and J in the action in equation (118) and demanding the reality condi 



tion for the action for u G [t/^, oo). The condition leads to the following two equations 

{GttGxx — E"^) = , 

B^e^^Gu - {GuGxx - E') Grr {det S,)GuGl, - e'^^J^) = . (121) 
The two expressions are the terms in the on-shell action that go to zero in the numerator 



and the denominator respectively. It can be shown from equation ( 121 ) that in order to 
have J 7^ one has to have B = 0, which corresponds to the straight brane solutions. 
For a given electric field, we can determine the position of the vanishing locus u^h 



from the first condition in equation ( 121[ ). Knowing u^h we can then extract the current 



J, and therefore the conductivity using J = crE. To express the conductivity in terms 
of physical quantities let us recall that 

fis = (27r)-«a'-5 , T 



47ri?3/2 • 
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Now we will restore the factors of (27ra') and also set R!^ = \^a' /A = 1. Combining 
everything we get the following expression 



This conductivity is due to the melting of mesons at high temperature and pair creation 
mediated by the electric field. Since we get this formula using DBI action, it captures 
non-linear behaviour of the conductivity with respect to the electric field. 

The effect of the electric field on chiral symmetry breaking should be to reduce the 
symmetry restoring temperature by polarizing the bound states into constituent quarks. 
The analysis above points to the fact that electric field works as expected from field 
theory perspective. However, the energy consideration does not lead to the expected 
result, because as electric field increases Uch also increases, therefore the straight branes 
which extend all the way down to Ut has more DBI action energy as compared to their 
curved counterparts, which can only extend down to u^h- We hope to address this issue 
in future. 

5.6 Concluding remarks 

We have extended the study of Witten-Sakai-Sugimoto model to include the presence of 
external electric and magnetic fields, examining the dynamics of the flavour sector in (an 
analogue of) the "quenched" approximation. We have seen that external magnetic field 
helps in chiral symmetry breaking. This particular effect of an external magnetic field 
has been referred to as magnetic catalysis in field theory literature (see e.g., refs. II113II ). 
Our results and observations are consistent with results from those approaches. We 
found that the chiral symmetry restoring temperature increases with increasing magnetic 
field. We have further observed that such holographic models have an upper bound for 




(122) 
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the symmetry restoring temperature depending on the dimension of the gauge theory 
and the probe brane. 

We briefly studied the effect of external electric field, and though we have not 
explored all of the details, we expect the djmamics to be also consistent with the field 
theory intuition, although (as we did for magnetic field here) the precise details should 
be interesting to uncover. 

Another important avenue would be to study the meson spectra in the presence of 
these external fields, which we pursue in the next chapter focussing on the magnetic 
field case. 
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Chapter 6: The Mesons in the 
Magnetic Field 

6.1 Introductory remarks 

In this chapter we will investigate the effect of an external magnetic field on finite 
temperature the meson spectrum in the Witten-Sakai-Sugimoto model. The material 
presented in this chapter is based on the work done in collaboration with Clifford 
Johnson[80]. 

Following the field theory literature, e.g. ref. nil3ll . we have briefly reviewed in 
Chapter 2 that for (3 + 1) and (2 + 1) -dimensional field theories there is a universal 
phenomenon at work that is sometimes referred to as the magnetic catalysis of chiral 
symmetry breaking. 

As we have seen in the previous chapters, this physics (that of a magnetic field cat- 
alyzing chiral symmetry breaking) can now be studied using the alternative technique of 
gauge/gravity duals (the magnetic field maps to a background B-field in the string the- 
ory and causes the flavour branes in the geometry to bend, generically), and it has been 
observed fTQl to persist in all the Witten-Sakai-Sugimoto type models, where the flavour 
degrees of freedom are inserted as defects in the dual geometry for several (2 + 1)- and 
(3 + 1) -dimensional gauge theories. We encountered a similar effect in chapter 3 in a 
different model, where some of the structure of the meson spectrum was also studied, 
revealing phenomena such as Zeeman splitting, level mixing, etc in ref. [jH. Here we 
would like to study the meson spectrum in the presence of an external magnetic field in 
the Witten-Sakai-Sugimoto model supplementing the study of the phase diagram dis- 
cussed in the previous chapter. 



99 



Previous studies of the meson spectrum in this particular model have been carried 
out in e.g., ref. [123 for both low spin and high spin mesons. A comprehensive review 
on holographic studies of mesons can be found in ref. [|47l . We uncover the role of 
the magnetic field in these spectra in this particular Type IIA model. For the low spin 
mesons we restrict ourselves to the high temperature phase only and therefore study 
the quasinormal modes of the scalar and a subset of the vector mesons. This comes 
with a caveat. The precise identification between the quasinormal frequency associated 
with the supergravity fluctuations in this model and the meson spectrum in the dual field 
theory is rather subtle and poorly understood. For a detailed analysis regarding the issue 
of precise identification, see ref. UllSL In the Witten-Sakai-Sugimoto type models the 
mesons transform under an U(A^/)diagj whereas the quasinormal modes transform under 
either an U(A^/)l or an U(A^/)/j. This fact alone readily makes it difficult to precisely 
define what gauge theory quantities we would be studying by extracting the quasinormal 
modes. Nevertheless the quasinormal modes are naturally associated to the embeddings 
falling into the black hole of the gravity dual and therefore we study those in their own 
rights. 

Moreover, since the embedding function of the probe D8-brane is trivial (tq{u) = 
const, where to(u) denotes the profile function) in the high temperature phase, the scalar 
and the vector fluctuations remain decoupled even in the presence of a background anti- 
symmetric field such as a constant magnetic field. This is to be contrasted with the D3- 
D7 model in a background magnetic field lfTOl where the profile function of the probe 
D7-brane is non-trivial and therefore the scalar fluctuations couple with the vector fluc- 
tuations. Thus studying quasinormal modes in presence of a background anti-symmetric 
field in general is more difficult in the later model. 

The study of high spin mesons is another avenue that we pursue. We find that the 
presence of the magnetic field enhances the stability of such mesons by increasing their 
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angular momentum. This is turn increases the dissociation temperature at which they 
fall apart into their constituents. We also uncover a simple realization of the Zeeman 
effect in the presence of a background field. It can be readily seen that such effects exist 
in all the Witten-Sakai-Sugimoto type models with probe branes of diverse dimensions. 



6.2 Mesons with small spin 

Mesons with low spin in the dual gauge theory correspond to the small fluctuations of 
the classical profile of the probe brane. The fluctuation in the geometric shape (which is 
denoted by the function t(u)) of the probe branes corresponds to the scalar (or pseudo- 
scalar) meson mode and the fluctuation of the Maxwell field on the worldvolume of the 
D8/D8-brane corresponds to the vector (or axial vector) mesons. 



6.2.1 Scalar fluctuation 

Let us consider a fluctuation of the probe brane embedding function t(u) = tq(u) + x, 
where to{u) is the classical profile. Before proceeding further let us recall that the 
function to(u) is non-trivial only for the U-shaped configuration (for the straight branes, 
this function is just a constant). Here we will consider the high temperature phase only 
for which tq(u) = const, and the scalar fluctuation is always decoupled from the vector 
fluctuations. 

To recast the equation of motion for the fluctuation mode into an equivalent 
Schrodinger equation we define, f^^{u) = a{u)g{u), where a{u) is defined in eqn. 



(C.2) (the details of the coordinate changes are given in appendix C). The equivalent 
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Schrodinger equation is given by (to avoid notational clutter we choose R = 1 and 

Ut = 1) 



dlg{u) + uj^g{u) — Vs{u)g{u) = , where du 



u 



3/2 



-du 



and Vsiu) 



1) {hH^ {7u^ + 5) + 2H^u^ {71u^ + 7) + IQu^ {5u^ + 1)) 



(123) 



Here u is the "tortoise" coordinate defined in appendix C. For simplicity, we continue 
to express the Schrodinger potential in the original w-variable. Now we can understand 

Vs(n) 
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Figure 18: The Schrodinger potential for different values of external magnetic field. The dashed 
vertical black line represents the position of the horizon, the top-most curve represents the poten- 
tial at zero external field. As the field increases we find the potential grows more slowly with 
increasing u. 



qualitative features of imaginary part of quasinormal frequency by studying this effec- 
tive potential. We resort to numerics to study the potential Vs{u). This potential vanishes 
at M = 1 because of the presence of the horizon and diverges as n — )• oo reflecting infi- 
nite gravitational potential of the background geometry. 



The important feature that we find from figure 18 is that the Schrodinger poten- 



tial is always positive, therefore from eqn. (C.6) we can conclude that the imaginary 
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part of the quasinormal frequency is strictly negative. Moreover, the potential does not 
develop a negative energy well, therefore there exists no negative energy bound states 
in this phase. This is consistent with the fact that the probe brane configuration is sta- 
ble and tachyon free and chiral symmetry restoration is accompanied by meson melting 
transition. Tachyonic modes in the scalar field fluctuation are likely to appear in the 
overheated low temperature phase which is bypassed by the transition. Also we see that 
the potential vanishes at the horizon and therefore the spectrum is continuous, which is 
consistent with the melting transition. 

Now we can observe that the Schrodinger potential increases monotonically as we 
move away from the horizon. This implies that Re[a;^] > and hence Abs[tjR] > 
Absfcui]. By explicit numerical calculation we will find that this is indeed the case (we 
point to figure 1 1 9(a)] and [T9(b)] ) . 



To extract more of the physics, we now compute the quasinormal frequency corre- 
sponding to this scalar fluctuation. We find it extremely convenient to perform a change 
of coordinates as implemented in ref. [|52l . Let us introduce a new variable p = u/Ut 
and V = t + a{p), where a{p) is determined from the condition that 

dv = dt+ . dp . (124) 

This coordinate system {v, p} makes the numerics conveniently stable to find the quasi- 
normal modes. We refer the reader to ref. [l52ll for more details. In presence of an 
external field this equation does not take the general form described in ref. [|52l . This is 
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due to the breaking of the S0(3, 1) -> S0(1, 1) x S0(2) by the presence of the mag- 
netic field as a result of which the Laplacian in the {x^, x'^}-direction gets squashed. 
The equation of motion for the fluctuation now takes the form 



2,-3^1/2^/] = 0. (125) 



lUJ 



[d,p"\l + H'p--'f''f^ + 2p''\l + H'p 



It is now straightforward to check that this equation admits the existence of quasinormal 
modes. To find how it depends on the background magnetic field we use a shooting 
technique to solve the fluctuation equation and pick the appropriate value of uj for which 
the equation admits a normalizable solution (this is achieved by imposing the condition 
that — 7- as p — 7- oo). As boundary conditions we impose /^(l + e) = 1 and 
/^(l + e) equal to a value obtained by requiring the equation of motion to be regular 
near p = 1, the event horizon. Here the parameter e is a vanishingly small number in 



our numerical scheme. The result of this pursuit is summarised in figure 19(a) and 
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Figure 19: The real (wr = Re[w]) and imaginary (wi = Im[a;]) part of the quasinormal modes 
for the scalar fluctuation as a function of the external field. The quasinormal frequency is mea- 
sured in units of the background temperature T. The small wiggles are due to numerical errors; 
here we have plotted the polynomial which fits the data. As expected before, we confirm that 
AbsfwR] > Abs[wi], however they are roughly of the same order. 
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19(b)[ Since we are in the high temperature phase the quasinormal mode is comprised 



of both real and imaginary parts. It would be expected that the absolute value of the 
real part corresponds to the mass of the meson before it melts and the absolute value 
of the imaginary part corresponds to its inverse hfetime analogous to ref. [TTBll for the 



D3-D7 system. Figure 19(a) then tells us that the magnetic field increases the scalar 



meson mass and figure 19(b) hints that the inverse lifetime gets bigger as the field is 
dialed up. However the precise identification of the meson mass and inverse lifetime 
with the quasinormal modes of the scalar fluctuation is a subtle issue in this particular 
model (see ref. ifTTSll ). 



6.3.2 Vector fluctuation 

The vector meson spectrum in general is determined by a set of coupled equations of 
motion. It is however possible to completely decouple the transverse and the longitudi- 
nal modes if we restrict ourselves to the oscillations parallel to the magnetic field. If we 
restrict ourselves to the oscillation in one of the perpendicular directions (e.g. along the 
x^-direction) to the applied magnetic field, then it is also possible to decouple one of the 
transverse modes (in this case the A2-vector mode). However, the remaining transverse 
mode and the longitudinal mode remain coupled. The details are given in appendix B. 
Here we study the decoupled sectors only. 
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6.3.2a The transverse mode 

We begin analyzing the transverse vector meson spctra. To extract qualitative features 
of the spectra we study the Schrodinger potential. The equation of motion for the A2 
vector modes is given by 



du 



e 



det u^S^'S^^A^iu) = 



(126) 

where we chose an ansatz of the form A2 = A2{u)exp{—iut) . It is straightforward to 
turn this equation in the form of a Schrodinger equation with an effective potential given 
by 

, , , , (m3 - 1) (5 (7u^ + b)H^-2 [u^ - 43m3) H'^ + S iu^ + 2u^)) 

^3/2 

with du = — du , (127) 

u-^ — 1 

where u is the "tortoise" coordinate. Clearly the horizon is located at m — > 00 and the 
boundary is located at u 0. The corresponding potential is shown in figure |20j The 
main qualitative feature that emerges from this is the fact that there is no instability due 
to existence of bound states of positive imaginary frequency modes since the potential 
remains positive definite. We observe, however, that a non-vanishing magnetic field can 
create a dimple in the otherwise smooth and rather featureless potential and therefore 
create local (positive) potential well. Therefore it is possible to have long-lived bound 
states with Absf^n] ^ Absfwi]. In this case, Absl^R] is set by the depth of the well, 
which in turn is set by the background field. If we increase the magnetic field even 
higher then the dimple gets bigger and bigger and ultimately swallows the well. At even 
higher magnetic field the potential again takes a shape similar to the vanishing external 
field case. Therefore the imaginary part of the quasinormal frequency may have an 
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Figure 20: The Schrodinger potential for the vector meson for different values of external mag- 
netic field. Here u denotes the "tortoise" coordinate mentioned in appendix C. As the magnetic 
field increases we observe that a small well develops in the effective Schrodinger potential. The 
depth of this positive well is set by the background field. 



interesting behaviour with increasing magnetic field. We may expect that Abs[a;i] starts 
decreasing with increasing magnetic field, but after a critical value it starts increasing 
again. Next we turn to the numerics. 

The equation of motion (for the vector mode A2 with non-zero spatial momentum) 
in this case is given by 



da 



det (E(o))S22 {S'"'Fp2 + S^''F^2 + S'^'Fa) 



0, ze{l,3}, (128) 



where S'^^ are the analogue of 5"* in eqn. (|B.3|) when the coordinate change described 



in eqn. (124) has been performed. 

The background field breaks the full S0(3, 1) Lorentz symmetry and there is an 
unbroken SO (2) symmetry corresponding to the rotation in the plane perpendicular 
to the magnetic field. Therefore the transverse vector meson spectra with a non-zero 
momentum will have two distinct branches. Mesons can have momentum along the 
direction of the background field (corresponding to i = 1) or in the perpendicular plane 
(corresponding to i = 3). 
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Without any loss of generality we can consider studying the spectrum of the A2 
vector meson. We can give this meson a momentum along the x^-direction or along the 
x^-direction (but no momentum along s^-direction since this is the transverse mode). 
The spectrum would be entirely equivalent to the spectrum of the ^3 vector modes 
having momentum along the -direction or along the -direction. 

Now we follow the same numerical approach. We impose yl2(l + e) = 1 and fix 
^2 ( 1 + e) from the equation of motion near the horizon. With these boundary conditions 
we look for normalizable solutions for A2, which gives the quasinormal modes. In 




(a) 



(b) 



Figure 21: The real and imaginary part of the quasinormal modes for the vector fluctuation as 
a function of the external field. The frequency is expressed in units of background temperature 
T and the magnetic field is expressed in units of R. We find a local minima in Abs[u;i] in the 
parameter range we have explored. This is qualitatively similar to the behaviour anticipated from 
analyzing the effective Schrodinger potential. We also observe that in the most of the parameter 
space Abs[cjR] » Abs[a;i]. 



figure 21(a) and 21(b) we have shown the dependence of the real and imaginary part of 
quasinormal frequency for the particular choice of = 0. We observe that the inverse 
hfetime of the vector meson gets shorter as the external field is dialed up and therefore 
enhances the stability of the meson. This behaviour is opposite to what we observed for 
the scalar meson fluctuation. 

Next we obtain the dispersion relation when k ^ 0. In figure |22(a) we demon- 
strate this dispersion relation for a wide range of values for k and in figure |22(b)| we 
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Figure 22: The dispersion relation between the real and imaginary part of the quasinormal 
frequency for two different values of magnetic field. The red (solid) curves correspond to the 
dispersion relation for non-zero magnetic field. The black (dashed) curve corresponds to the 
dispersion relation when there is no magnetic field. The black dots represent the corresponding 
values of quasinormal frequency for vanishing spatial momentum. Figure |22(a)| and [22(b)|cor- 



respond to vector mesons having momenta along -direction and figure 22(c) corresponds to 
vector meson with momentum along -direction. 



have shown a magnified version of figure 22(a) to more clearly see the role an external 
magnetic field plays in this case. These two figures correspond to the dispersion rela- 
tion when the momentum is parallel (i.e. along a; ^-direction) to the magnetic field. We 
observe that for high enough spatial momentum the dispersion relation curve approx- 
imates the zero field dispersion curve quite well. We have shown in figure |22(c) the 
dispersion relation when the momentum is perpendicular (i.e. along -direction) to the 
background field. The qualitative behaviour between these two dispersion curves are 
clearly different; however for large enough spatial momenta these curves tend to become 
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insensitive to the background field. This is because the heavier the meson becomes the 
less sensitive it is to the background field. 



6.3.2b The longitudinal mode 

Now we study the longitudinal mode. Bearing in mind the change of variables in equa- 
tion ( 124| ), we define the gauge invariant longitudinal mode to be given by 



ujAi + A;A„ = E . (129) 

Now the equation of motion for this longitudinal mode is given by 

E"(p) + Zi(p)E'(p) + Z2(p)E(p) = , (130) 

where Zi and Z2 are known functions of p. The explicit expressions of Zi and Z2 are 
not very illuminating, thus we do not provide them here. 



Equation (130) now can be numerically solved by using the boundary conditions 
that E(l + e) = 1 and E'(l + e) equal to a value obtained from the equation of motion 
itself. The condition of normalizability then fixes the quasinormal modes. 



The numerical result is shown in figure 23(a) Note that here we numerically explore 



only for positive values of H. The equation of motion for the longitudinal modes have 
a linear term in H, and therefore the quasinormal modes may depend on the sign of H. 
However, we found numerically that the qualitative features of the quasinormal modes 
are similar also for negative values of H. The qualitative behaviour of the dispersion 
curve is similar to the transverse modes. It is evident that increasing magnetic field 
brings quasinormal frequency closer to the real axis in the complex plane. However, for 
a fixed magnetic field, the dispersion curves tend to approach the zero magnetic field 
dispersion curve as the momentum increases. 
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Figure 23: The dispersion relation of longitudinal oscillation in presence of a background mag- 
netic field. The quasinormal frequency is measured in units of the background temperature. In 



figure 23(a) the black solid curve represents the dispersion relation when the magnetic field is 
absent, the blue (circular) dots correspond to H = 5 x 10^^ and the purple (square) dots cor- 
respond to H = 10~^. H has been expressed in units of R. In figure 23(b) we have shown 
the dependence of the lowest quasinormal frequency corresponding to the hydrodynamic mode 
with the magnetic field for a fixed momentum k = 0.01. The imaginary part of the quasinormal 
frequency does not depend on the magnetic field in the hydrodynamic limit. 



6.4 Dispersion relation in the hydrodynamic Hmit 



In this section we wish to explore the dispersion relation of the longitudinal mode in the 
hydrodynamic limit. Recall that the longitudinal gauge field fluctuations correspond to 
the diffusive mode in the dual gauge theory. 

Going back to the original {t, u}-coordinate and defining the gauge invariant lon- 
gitudinal mode £{u) = cuAi + kA^, the equation of motion is given by (a la equation 



(B.9)) 



£"{u) + Zi{u)£\u) + Z2{u)8{u) = 



(131) 
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where 



%Hk(T (5=^=^ (8//A;(7 - (5")') + {SHau/k + £5"" (5^^)')) 

(132) 



with 

£(«)=^,e-*^-det(E(o)), a = (133) 
We also have the following constraint 

(cjA;^** - fcA'i^^^) + 8(7//£: = . (134) 

Now we venture for the dispersion relation in the hydrodynamic limit. In this par- 
ticular set-up, for simplicity, we will set the radius of the near-horizon D4-brane back- 
ground to unity (Ut — 1) and also set the horizon radius to unity {R — 1). We will 
restore these factors at the very end. To facilitate the hydrodynamic limit, we will 
assume a; <^ 1 and A; <^ 1 and further assume that H <^ 1. The last condition is 
imposed assuming that the length scale associated with the magnetic field is small com- 
pared to the plasma temperature. To further quantify this limit, we will parametrize 
H = kh holding h fixed in the hydrodynamic limit. 
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In this limit we will attempt to solve the equation ( 131 ) perturbatively in k. To this 
end, we write the following expansioij^ 



S{u) = {u- 1)-'"/=^ [J^miu) + kj^^i){u) + A;V(2)(n)] + O (P) . (135) 



Motivated from the hydrodynamic result in (E.22 ) we parametrize the dispersion relation 
in the following manner 



(136) 



where 71, 72 and 73 are three constants to be determined. 

We will attempt to determine these constants by perturbatively solving the differen- 



tial equation ( 131 1 subject to the following Dirichlet boundary conditions 



\mi£(u) = 1 , lim S(u) = . 



(137) 



These boundary conditions determine the so-called quasinormal modes of the corre- 
sponding fluctuation spectrum. Recall that when /i = 0, we have the known dispersion 
relation u = —i{2/3)k'^ [80] and hence we expect to obtain 71 = 2/3. 

Before proceeding further we have to fix the boundary conditions on the functions 



J^(^i). Since the differential equation in (131 1 is homogeneous, we can impose the fol- 
lowing boundary conditions without any loss of generality 



lim J^i 



n— >1 



(0) 



lim J'(i),(2) = 



lim J'(o),(i),(2) = . 

n— s>oo 



(138) 



^To directly compare to the notation of Appendix B of ref. fSOl, we need to identify Cj — uj/3. 
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Now we substitute the expressions in (|135[) and (|136[) in the fluctuation equation 



( 131 1 and analyze the corresponding equations at various orders of the expansion param- 
eter k. 

At the leading order in k, we get the following equation 

+ du [log J-('o)(n) = . (139) 

Note that there is no J^(o) term present in the above equation. The solution of this 



equation subject to the boundary condition in (155 1 is given by 



^(0) = ^ , (140) 

where we have used the specific functional forms of C, 5** and 5"". Clearly at this order 
we do not get any information about the dispersion relation. 
At the next order, we get the following equation 

J-(';)(iz) + [log J-('i)(m) = , (141) 

which is exactly similar to the equation satisfied by the zeroth order solution J^(o). How- 
ever, as we have arranged, satisfies a different boundary condition near the horizon 
as compared to the solution J-'{o). It is straightforward to check that in this case we get 
^-"(1) = 0. At this order, we do not get any non-trivial information. 

The differential equation at the next order is the most interesting one. The general 
structure of the equation is obtained to be 

J-("2)(w) + du [log {CS^'S^'')] F[,){u) + J (j-(o), J-('o),-F(i), = , (142) 
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where J sources an inhomogeneous piece constructed from the solutions obtained at 
lower orders]^ Unlike the lower orders, equation ( 142 1 is not exactly solvable, however 



we can solve it analytically near the horizon and near the boundary. 

Near the boundary, as u — t- oo, we can always obtain a solution for T{2) which 



is consistent with the boundary condition in (155). The non-trivial information comes 
from solving the above equation near the horizon. In this region the function T(2) takes 
the following form 

^■(2) = finite + (7i + /i(/i72 + 273)) (872/^^ + ^(373 + 16a)/i + 871 - 2) log(M - 1) , 

(143) 

where the "finite" piece consists of one undetermined constant of integration and can 
be arranged to be zero. Thus we need to separately impose the coefficient of the log- 
divergent term to be zero. This gives the following algebraic condition 

either 71 + /i {h^2 + ^3) = , 

or 872/1^ + i (873 + 16(t) /i + 871 - 2 = . (144) 

Interestingly, if we treat h to be an independent parameter, then the algebraic constraints 
impose numerical values for 7j's at different orders of h. Clearly the first algebraic 



condition implies 71 = = 72 = = 73 and hence gives a; = by virtue of (136) 
However, the second algebraic condition imposes 



2 16(T 
7i = o ' 72 = , 73 = — , (145) 



^The full expression of J is rather complicated, therefore we do not present it here. 



115 



which ultimately leads to the following dispersion relation 



CO = -i-P ± —kH . (146) 
3 3 



Comparing (146) with (E.22) we get 



V _2 _16a_ ^^^^^ 



e + P 3 ' 2 (e + P) 3 ' 

Thus we find that the magnetic field sources a propagating (sound) mode in the otherwise 
diffusive channel. In the hydrodynamic paradigm, such a mode is clearly sourced by 
the anomaly associated with a global U(l) symmetry as demonstrated in appendix E. 
It is also obvious from the presence of the Chern-Simons coefficient a in the above 
formulae. Finally restoring the factors of temperature, we get the following dispersion 
relation (after setting R = 1) 

1 72 

u = -i—;P± .kH . (148) 

2vrT {4:7iTf 

6.5 Mesons with large spin 

Mesons with large spin can be described by strings having both its end points on the 
flavour brane and rotating on a two plane along the directions where the gauge theory 
lives. For large enough angular momentum, the length of the string is much larger 
compared to is and we can use classical Nambu-Goto action to describe the dynamics 
of the rotating string. Mesons with large spin have previously been analyzed in, e.g., 
refs. [|122[|98l . Here we will follow the same framework, now including the effect of an 
external magnetic field. 
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In absence of magnetic field due to rotational symmetry, the spectrum does not 
depend on the direction of the angular momentum J. An external magnetic field breaks 
this symmetry and introduces a splitting of the energy levels (Zeeman effect). Here we 
will focus on the two cases, namely when the magnetic field is parallel or anti-parallel 
to the angular momentum. In this case there should be a Zeeman splitting in energy for 
mesons with non-zero magnetic moment. 

This idea has been implemented in ref. [75] in a different model. We adopt similar 
system and consider two flavour brane-anti-brane pair very close to each other so that in 
the bulk they join leaving two diagonal U(l)s (call this U(1)a and U(1)b corresponding 
to the two pair of flavour branes). We consider now turning on fields having equal 
opposite charge under these two U(l)s. Now AA and BB strings will have equal and 
opposite charges but AB and BA strings will have equal charges. In the first case, 
the total orbital magnetic moment of the meson will vanish whereas in the second case 
the meson will have a non-zero magnetic moment. The magnetic field will induce a 
quadratic correction to the meson energy for the former and a linear Zeeman splitting 
for the later. Here we will analyze the later case only. It is useful to note that this 
configuration is symmetric under reflection around the midpoint of the string (at p = 0). 
We will refer to this as the symmetric configuration. 

To proceed we write the relevant part of the background metric in the following form 
(rewriting the {x^, x'^}-two plane in polar coordinate {p, 0}) 



Here we have explicitly written down the metric corresponding to the high temperature 
phase only. The background metric corresponding to the low temperature phase can 




(149) 
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likewise be written in an analogous form. We recover the zero temperature background 
by setting Ut = 0. We will work with the Nambu-Goto action for the string with the 
following ansatz for the string profile 

t = T , p = p{a) , u = u{a) , (p = UT . (150) 

From now on the parameters {t, cr} will always refer to the worldsheet coordinates for 
the string. We also assume that cu is positive {i.e. only clockwise rotation). 

In the presence of an external constant magnetic field the Nambu-Goto action is 
accompanied by a boundary term coupled to the magnetic potential = Hp^ /2. This 
potential gives rise to a magnetic field B{2) = Hpdp A d(f) = Hdx^ A dx^. Taking 
this boundary term into account the action for the string is given by 




where represents the right and the left boundaries and the relative positive sign 
between the two boundary contributions coming from the magnetic potential is due to 
the symmetric configuration. For the non- symmetric configuration the there will be a 
relative negative sign between these two terms. We refer to ref. HTSll for further details 
on this construction. For convenience we introduce the rescaled variables 

p = puj , u = uco , R = Ruj , a = cru . (152) 
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The expression for the energy and the angular momentum for the spinning string in the 
rescaled variables can simply be obtained to give 



- T - ^ 



J 



AJ, 



1 



1 



dL 

1 



da 



?,2 



~3 I ~/2 I ~/2 
U \ p +U 



('1 






a+ ^ 




We observe that the external field can enhance angular momentum of the meson if it is 
aUgned in parallel (positive values of H)\ and can reduce the angular momentum if it is 
anti-parallel (negative values of U). 

To proceed we need to solve for the string embeddings. The equation of motion 
obtained from the Nambu-Goto action has to be supplemented by the boundary condi- 
tion 



dL 



TT 



d{X<) 



M 



sx 



M 







(154) 



as 



Since Su\qj. — and 5p\qj, is arbitrary, we have to impose the Neumann boundary 
condition (in absence of any external field) {dL/dp')\gj. = 0. This condition gives us 
the following constraint 



{u/RuY (1 - p 
p'2 + u'^ (ujR/u) 



1/2 



. 



(155) 
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So we need to impose the boundary condition p'\q-^ = 0, which means that the string is 
ending perpendicularly on the flavour brane-anti-brane pair. We will come back to the 
details of our numerical scheme in a later section. 

The presence of an external magnetic field changes this boundary condition. It can 
be easily seen that the boundary condition to be satisfied is now given by 



IT 



.11 



-H~p) = 



(156) 



where vr^l^^ is given in equation (155). Therefore we need to solve the equation of 



motion obtained from the action in eqn. ( 151 1 subject to the boundary condition given 
in eqn. ( 155 1 or (156]) in absence or presence of a magnetic field respectively. 

We will begin by obtaining some analytical results in the large angular frequency 
hmit in the zero temperature phase. We will choose a = p gauge and send a; — )■ oo; 



from the reality condition of the Nambu-Goto action in eqn. ( 151 ) this limit sets a bound 
IpI < l/w — 7- 0; and therefore equivalently we study the short string limit. 



6.5.1 Analytical results 

In this section we consider the zero temperature background and spinning strings with 
large angular frequency only. 

6.5.1a Vanishing magnetic field 

In this section we closely follow the approach adopted in ref. [|98l . In the large angular 
frequency limit the relevant string profile is well approximated by the local geometry 
very close to the point from where the string hangs. In this limit the string hangs from 
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the point where the brane-anti-brane pair join (the radial position denoted as Uq). We 
therefore pick the following ansatz to approximate the string profile 

u{p) =ujUo + ^ . (157) 

CO 

It is easy to check that this ansatz is consistent with the equation of motion for the string 
in a; — 7- oo limit. To obtain the leading order behaviour we have to solve for the function 
f{p) with the boundary conditions that /(O) = = /'(O). At £t = cjUq we require 
u'{p) = 0, which sets the latter boundary condition for the function /(p). Substituting 



the ansatz given in eqn. ( 157 1 in the equation of motion for the string and keeping only 
the leading order terms we obtain a differential equation for the function /(p) which can 
be integrated to obtain the following analytic form 

/(P) = ^(p' + arcsin2(p)) . (158) 



Plugging back /(p) from eqn. ( 158 1 to eqn. ( 153 ) and simplifying we get the following 
relations 

(159) 

From this we can see that a; — )■ oo limit is equivalent to J ^ A limit. Also we find 
that the mesons with large angular frequency follow a Regge trajectory with an effective 
tension 

Teff = £;V(2vrJ) = XU^^'/{2R'/') . (160) 

The radial position Uq is related to the asymptotic separation L between the flavours via 
[/J/2 ^ R^/^/{AL)B{9/16, 1/2), where B{a, h) is the Beta function. This length L sets 
the coupling strength of the non-local four Fermi interaction for the dual NJL model. 
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6.5.1b Non- vanishing Magnetic Field 



In presence of a magnetic field the leading order solution for u{p) remains the same as 



in eqn. ( 158 1; however the boundary condition needs to be modified as given by eqn. 



(156). Using the expression for tt^I^^ from eqn. (155) (in the gauge cr = p) in the 



modified boundary condition given in eqn. ( 156) we can determine the value of p at the 
boundary (meaning at n = f/o where the string ends on the flavour brane). In the limit 
of small magnetic field this critical value of p (denoted as pc) can be obtained to be 



PC 



IP 



R 



(161) 



With this pc, there will be two different string profiles which we have shown in figure 



24(a) and 24(b) respectively. 
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Figure 24: The two possible profiles are presented; Ug is the radial turnaround position for the 
string. These are obtained by solving the equations of motion numerically, although we come 
back to the numerical results in a later section. The profile in the left hand side has more energy 
than the one in the right hand side. 



These two different profiles result in two different energies (let us denote these ener- 
gies by E±). The one with the lower energy {i.e., with E_) is extended from p = to 
pc and the one with the higher energy {i.e., with E+) is extended from p = to p = 1 
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and then folds back io p = pc- Using the formulae for meson energy E and angular 



momentum J given in eqn. ( 153 1 we get the following linear order corrections (Zeeman 
splitting) introduced by the magnetic field 



E. 



E, 



3/4 



J19/4 



'\J - 



^9/4 VAJ+ — 



^ V A f/3/4 



(162) 



Recall that the parameter Uq is also related to the constituent quark mass. For a study of 
how the constituent quark mass behaves (for a fixed value of L) with the magnetic field 
see ref . [79,1 . 

It is easy to check that this splitting is consistent with the general formula for energy 
(obtained in ref. [IJ) of a Nambu-Goto string with an effective tension moving in a 
background magnetic field. The energy spectrum (for large J) is given by the following 
formula 



E^ = {27lT,s){l - k)J 



2 f H\ 

— arctan — 

Vr \TesJ 



(163) 



Recalling the expression for Tes from eqn. (160) and using the general result above the 
hnear order correction to the energy of the string can be obtained which matches exactly 



with what is obtained in eqn. ( 162). 



6.5.2 Numerical results 

In order to extract the relevant physics we need to solve for the string embedding. To 
that end we choose the parametrization a = p + u since it provides better numerical 
stability as discussed in e.g., ref. fTSl . 
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To solve the equation of motion first we fix the maximum radial position the string 
can attain (this is the position where the brane-anti-brane pair join, which is denoted 
by Uq). This in turn fixes the asymptotic separation between the brane-anti-brane pair. 
Fixing Uq is equivalent to fixing the constituent quark mass of the dual gauge theory 
and fixing L, the asymptotic separation of the brane-anti-brane pair, is equivalent to 
keeping the four fermi coupling of the dual NJL model. Now for a given Uq we look 
for a radial turnaround position for the string (call it Ug) such that shooting from Ug 
with the IR boundary conditions u(a)\jj^ = Us and u'{a)\jj^ = would satisfy the 
UV boundary condition u'{a)\jj^_^ = 1 (up to numerical accuracy). The existence of the 
turnaround position is guaranteed by the symmetric configuration that we consider here. 
Below we show two such representative profiles for to = 0.5 and a; = 3 respectively 
(corresponding to the blue and the red curve). 

uicr) 

___ 
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Figure 25: Two representative spinning string profiles. We have set i? = 1. 

Such profiles were discussed in ref. ["981 in a different but closely related model. In 
brief, we find that for a given value of u there are classes of solutions distinguished 
by the number of times they cross zero (referred to as "nodes" in ref. [.98.1 ) along the 
horizontal axis. Here we will constrain ourselves to the n = 1 case only. 
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We will now explore the behaviour of high spin mesons in presence of a background 
magnetic field. The key feature that we will observe for both zero and finite tempera- 
ture phases is the existence of a shift in physical quantities such as energy and angular 
momentum in presence of a non-zero magnetic field. We begin with our numerical 
results for the zero temperature phase. 

6.5.2a Low temperature phase 



E 




Figure 26: The effect of an external magnetic field in the zero temperature phase. The black 
dashed lines represent H = 0, blue dashed lines represent H = —0.1 and the red dashed lines 
represent H = 0.1; E and J have been evaluated in units of 1/(2110'). H have been evaluated 
in units of R. 



In figure 26 we have shown the dependence of the meson energy as a function 
of the angular momentum of the meson. For zero external magnetic field we recover 
the behaviour discussed in ref. [|122i |. We find that the external magnetic field shifts the 
energy and the angular momentum of the meson for generic values of the parameters 
of the system; however with increasing angular frequency these observables tend to 
become insensitive to the external field. In view of the approximate analytic results 



obtained in eqn. (162) it is straightforward to observe that an external field contributes 
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only at order H /u and hence has vanishingly small effects for large enough angular 
frequency. 

The existence of two energy branches for the same given value of J seems to persist 
(as in vanishing external field) even in the presence of a magnetic field, although as the 
magnetic field is increased we observe a tendency that this multi-valuedness in energy 
starts disappearing. Therefore any possible unstable upper branch in figure 26 can get 
promoted to a stable one by having sufficiently high angular momentum. 

Moreover there exists a maximum angular momentum beyond which the spinning 
meson dissociates. This dissociation is mediated by the acceleration of the spinning 
string. However the presence of magnetic field can again stabilize the mesons by rais- 
ing the maxima. Physically this simply corresponds to the fact that the magnetic field 
enhances the angular momentum of the meson and thus makes it stable. 
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Figure 27: Dependence of maximum spin with magnetic field for varying constituent quark 
mass. Jcrit and niq have been calculated in units of l/(27rQ'). We observe that even for a 
certain negative value of H, i.e. when the magnetic field is aligned anti -parallel to its angular 
momentum, the meson can be more stable than its zero field counterpart. 



To demonstrate this effect we observe the plot in figure 27 Clearly the H = —0.1 
branch is higher than the H = branch; an even higher value of H would inhibit any 
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dissociation at all. Also the heavier the constituent quark mass is the more J^it becomes 
making it less likely for the heavy mesons to dissociate. 
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Figure 28: Linear correction to energy for the equal charge case. From top to bottom the values 
of angular momenta are J = 0.937, 0.907, 0.862 respectively. The dashed straight lines are the 
best fit lines. E and J are both presented in units of 1/ {lira'). We observe that all the three best 
fit lines have same (up to numerical accuracy) slope, which we expect from the very nature of 
Zeeman splitting. 



A key feature of the high spin meson spectrum in presence of an external magnetic 



field is the linear Zeeman splitting. In figure 28 we have shown such linear correction 
in energy for fixed values of the meson angular momentum. 



6.5.2b High temperature phase 

We now perform a similar numerical study of high spin mesons in this intermediate 
temperature phase. 

In figure [29] we have shown the dependence of meson energy as a function of its 
angular momentum. The qualitative features of finite temperature physics are similar 
to that of the zero temperature physics. There is however one distinction as compared 
to the energy spectrum at zero temperature phase. The presence of magnetic field at 



finite temperature does not promote the possible unstable upper branch in figure 29 to a 
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Figure 29: The shifts in presence of an external magnetic field. E and J have been computed in 
units of 1/ (27ra'). H has been evaluated in units of R. The black curve corresponds to vanishing 
magnetic field, the blue curve corresponds to setting H = —1 and the red curve corresponds to 
setting H = I, where H has been expressed in units of R. 



stable one within the range of values for the magnetic field that we have explored using 
our numerical approach (we expect however that for sufficiently high values of H this 
unstable mode will be promoted to be a stable one). A possible thermally enhanced 
decay channel therefore remains open at finite temperature for a rather high value of the 
external field. The other familiar role of the magnetic field, we again find, is to introduce 
shifts in the physical quantities such as the energy and angular momentum of the meson. 



The meson dissociation due to spin should now be enhanced due to thermal fluctua- 
tion. We have previously observed that the magnetic field adds extra angular momentum 
to the system and stabilizes the mesons. We study in figure [30(a)] the effect of two com- 
peting parameters on the critical angular momentum of the meson. Finally we observe 



the familiar Zeeman splitting in figure 30(b) 
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Figure 30: Figure 30(a) sliows the dependence of maximum spin calculated in units of 1/ {lira') 
with the applied magnetic field. We can observe that there exists a certain range of values for 
Uq, which in turn fixes the value of the constituent quark mass, where the angular momentum is 



lowered by the magnetic field. Figure 30(b) shows the linear correction to meson energy for the 
equal charge case at finite temperature. From top to bottom the values of angular momenta are 
J = 0.630, 0.560 and 0.490 evaluated in units of l/(27ra') respectively. The dashed Unes are 
the best fit lines. 



6.6 Concluding remarks 

We have uncovered a sub-sector of the meson spectrum in the presence of an external 
magnetic field, complementing the phase structure that we obtained in ref. [|79l . Mesons 
with small spin that can be obtained by studying the quadratic fluctuations of the probe 
brane classical configuration are harder to study in absolute generality. In general the 
task is to solve a set of coupled differential equations, which we have here analysed 
in the situations where it is possible to decouple a subset of the modes. We have also 
analysed the large spin meson spectra in presence of a magnetic field to realize the 
well-known Zeeman effect. We also found that an external magnetic field enhances the 
stability of mesons and inhibits the dissociation. In appendix D we have presented a 
model calculation analyzing spinning strings in Rindler space to capture key features of 
the large spin meson dissociation. 
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There are several directions to pursue in future work. It is of interest to obtain the 
low spin meson spectra in the low temperature phase (i.e., when chiral symmetry is 
broken) and their dispersion relations. In the present work, we have focussed on the 
effect of a magnetic field. The analysis of the meson spectra and quasinormal modes in 
the presence of an electric field could be pursued in a similar spirit. The effect on the 
spectra of the presence of a chemical potential would also be of great interest. These 
external fields can all be realized as the non-normalizable modes of anti-symmetric 
fields on the world volume of the probe brane, and so we would expect that in all these 
cases in the high temperature chirally symmetric phase, the scalar and vector fluctuations 
always remain decoupled, as we have seen here. 
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Chapter 7: Conclusion 



We have extensively used the gauge-gravity duality to learn about the strong coupling 
dynamics of fundamental flavours. In particular we have used the framework of two 
different holographic models which share some rather generic features of QCD and are 
therefore expected to teach us about physics of the real world. In this thesis we have 
focussed entirely on the fundamental flavour sector. Our studies mainly included the 
effect of an external electric-magnetic field on the phase structure of the flavours at finite 
temperature, associated thermodynamics, chiral symmetry breaking and the associated 
meson spectrum. 

The two holographic models we have discussed are structurally quite different. In 
retrospect, let us make a comparative note. In the Type IIB model we studied the dynam- 
ics of a probe D7-brane in the background of the near-horizon limit of a stack of Nc 
D3-branes. At zero temperature this construction preserves N = 2 supersymmetry [[831 , 
however finite temperature breaks supersymmetry completely. In this model, we have 
identified the U(1)r symmetry corresponding to the rotation group in the transverse 2- 
plane of the probe D7-brane as the chiral symmetry. At zero temperature the presence 
of supersymmetry is enough to guarantee that the quark condensate must vanish at zero 
bare quark mass as mentioned in ref. [[TSl . Thus there is no spontaneous breaking of chi- 
ral symmetry in this model at zero temperature. Moreover, at zero temperature and zero 
bare quark mass this theory is still conformal in the limit Nf <^ N^. Therefore there 
is no scale in this system and thus no phase transition can occur. The zero temperature 
physics is somewhat trivial, although we can analyze many aspects such as the meson 
spectrum completely analvtically [|98l [99|[ . At zero temperature, the chiral symmetry can 
be broken explicitly by adding a mass term for the quark. In the gravity dual picture this 
mass term corresponds to a separation of the probe D7-branes from the D3-branes in 
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the transverse 2-plane, which can be introduced naturally in this model. This mass term 
corresponds to a non-normalizable mode in the profile of the probe brane and using 
the AdS/CFT dictionary we can read off the corresponding normalizable mode as the 
condensate. 

Let us contrast this with the Type IIA model. Unlike the Type IIB case, this is not 
a conformal theory. In this case, we study the dynamics of a probe D8/D8-brane in the 
background of the near-horizon limit of a stack of Nc D4-branes compactified on a spa- 
tial circle along the direction denoted by x^. By construction we impose anti-periodic 
boundary condition for the adjoint fermions which breaks supersymmetry. The global 
symmetry group is U(A^/) l x \][Nf)ii, which we continue to refer as the chiral symme- 
tr)Q The probe brane-anti-brane pair intersect the D4-branes at = ±L/2, therefore 
the quarks are massless by construction. Introducing mass term for the fundamental 
flavours and identifying the chiral condensate is rather subtle in this model. Work along 
this direction was initially carried out in refs. [|25ll42l and later in ref. flU. 

Let us recall that for all practical purposes the so called instructive limit of flattening 
out the a;^-direction (by sending R4 — t- 00) has almost no consequences in our gravity 
calculations. In this limit the zero temperature background does not break supersymme- 
try. However, introducing flavours by adding the brane-anti-brane probe immediately 
breaks all the supersymmetry. Therefore in this model the flavour physics is not influ- 
enced by supersymmetry at all. 

In the low temperature phase (i.e. the dual gauge theory is in the confining phase), 
when x^ is a compact direction the brane-anti-brane pair join together forced by the 
cigar topology of the background in the {x"^, wj-plane. This implies that the asymptotic 
global symmetry group U(A^/)l x \]{Nf)ji breaks down to a U(A^/)diag, which is very 

^Note that our definition of chiral symmetry in the Type IIB and the Type IIA set-up are quite different, 
however in an abuse of language we often do not make this distinction. 
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similar to the spontaneous chiral symmetry breaking in QCD. In the limit R4 — )• 00, 
the confining phase disappears. Furthermore, the background topology in the {x^, u}- 
submanifold is not cigar-shaped anymore. Therefore the flavour brane-anti-brane pair 
can either join together or end separately in the bulk. However it was shown in ref. [fTT| 
that it is energetically favourable for the probe brane-anti-brane pair to join together. 
Thus in this model, in the low temperature regime chiral symmetry is spontaneously 
broken. This is in sharp contrast with the Type IIB model. 

It was observed in ref. f53l| that in the Type IIB set-up an external magnetic field 
breaks all supersymmetry and therefore can induce spontaneous chiral symmetry break- 
ing. Furthermore this external magnetic field provides a natural scale for this system at 
zero temperature given by Ry/H, where R is the radius of the AdS-space and H is the 
strength of the magnetic field. Thus conformal invariance is completely lost even at zero 
temperature. A naive dimensional analysis can be used to fix the behaviour of the chiral 
condensate c ~ R^H^/"^, which we have discussed in the T — )• limit in chapter 3. 

The background temperature is another important physical scale. In the absence of 
any external magnetic field we discussed the temperature driven meson-melting transi- 
tion and associated meson spectrum in chapter 3. This transition can simply be described 
as follows. Below a certain critical temperature quarks exist in bound states such as 
mesons. Above this critical temperature, these bound states melt and form a plasma 
of the constituents. In chapter 3 we also investigated the case when both temperature 
and magnetic field are present. The phase structure is now controlled by two physical 
scales: the temperature and the chiral symmetry breaking scale induced by the external 
magnetic field. We have found that these two physical scales have competing inter- 
play and discussed the relevant phase structure in detail in chapter 3. The main feature 
that emerged from this study is the so called effect of "magnetic catalysis" of chiral 
symmetry breaking. Namely, for a given temperature sufficiently high magnetic field 
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induces chiral symmetry breaking by promoting the stability of quark bound states such 
as mesons or the chiral condensate. This means that for a given temperature we can 
increase the strength of the magnetic field so that the meson melting transition disap- 
pears altogether. In this model, this effect exists for arbitrarily large temperature. 

In the Type IIA model, the low temperature phase breaks chiral symmetry spon- 
taneously. Therefore introducing a magnetic field is expected to further enhance this 
effect. This is what we observed in chapter 4. However, the finite temperature physics is 
much richer in content. Sufficiently high temperature can trigger the restoration of chiral 
symmetry, which is similar to the meson-melting transition in the Type IIB set-up. Inter- 
estingly we found that an external magnetic field and temperature again appear as two 
competing parameters resulting in a non-trivial phase diagram. An external magnetic 
field promotes the breaking of chiral symmetry and avoids the melting transition. How- 
ever, in this model there seems to exist an upper bound of the chiral symmetry restoring 
temperature for arbitrarily large magnetic field. The existence of this upper bound in 
temperature is also in sharp contrast with what we observed in the Type IIB set-up. 

On a related note, in chapter 6 we have discussed the effect of the magnetic field on 
associated spectrum of mesons with both small and large spins. We have observed that 
the magnetic field increases the lifetime of the melting meson and enhances the angular 
momentum of mesons with large spin. Both these effects favour the bound states over 
the free constituent fundamental matter. 

There is an interesting point to note. As we have discussed, both these gravity duals 
are formally quite different. Even the breaking of the symmetry group that we have 
studied under the name of chiral symmetry is rather different. However, the qualitative 
feature of the "magnetic catalysis" seems to be a robust one in both these models. It 
is tempting to conjecture that this "magnetic catalysis" in chiral symmetry breaking is 
indeed an universal phenomenon. Similar conclusion has been arrived previously by 
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studying several field theory models in various non-perturbative frameworks (e.g. using 
the Dyson-Schwinger equation) in e.g. refs. [|68lll42llll3ll . It is perhaps not impossible 
to imagine probing some of these properties of the real quark-gluon plasma within an 
experimental set-up with an external magnetic field. 

Emboldened by the interesting physics observed in the presence of an external mag- 
netic field, it is natural to consider the effect of an external electric field. We have 
analyzed the effect of a constant electric field (in the Type IIB model) in chapter 4 and 
briefly towards the end of chapter 5 (in the Type IIA model). Intuitively, an electric field 
tries to dissociate any quark bound state into its constituents and drives a current. This 
is what we observed in chapter 4. The probe brane analysis in both these models seems 
to have an interesting feature. It is the existence of a vanishing locus or the so called 
pseudo-horizon (above the actual horizon) at which the probe brane action vanishes. In 
order for the probe brane to extend beyond this point, an additional gauge field is excited 
on the worldvolume. This additional gauge field gives the expectation value of the cur- 
rent that is caused by the electric field. We have computed the conductivity in both the 
Type IIB and the Type IIA models. 

However, much work is still needed to understand the curious features in presence 
of an electric field. In chapter 4, we found that depending on the value of the polar angle 
of entering the vanishing locus the probe brane develops a conical singularity. Similar 
conical singularity was observed in the Type IIA model in ref. [1231 . We have conjectured 
that stringy corrections may possibly cure this singularity. Another relevant avenue that 
we have not pursued in this thesis is the analysis of the meson spectrum in presence of 
the electric field. 

Following our general approach, we can further analyze the effect of a chemical 
potential in both these holographic models. The phase structure of thermal J\f = A 
Yang-Mills was analyzed in presence of a finite baryon density in ref. [|94| ; similar 
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study was carried out in the Type IIA model first in ref . [72J and later in refs. [|121[ |38l 
I221 11361 11191 121. It is natural to further introduce an external electric-magnetic field 
and investigate their effect. It turns out to have a rich phenomenology. Turning on both 
magnetic field and chemical potential has many interesting features, e.g. it induces an 
axial current in the Type IIA model. Work along this direction has been discussed in 
refs. [I24lll52[|l03lll04ll . This phenomenological richness is mainly sourced by the non- 
vanishing Chern-Simons term supported by the worldvolume of the probe branes. This 
non- vanishing Chem-Simons term is equivalent to the inclusion of the axial anomaly in 
the flavour sector of the dual field theory. It is further noteworthy that the effect of such 
axial anomaly terms within QCD has been recently studied in the context of the RHIC 
physics and new phenomena, e.g. the Chiral Magnetic effect [|56l , have been predicted. 
It is rather encouraging that similar effects are also observed within the Type IIA model. 

The role of finite chemical potential was studied within the Type IIB model in e.g. 
refs. [|94lll09ll and later generalized to isospin chemical potential in ref. [|48ll . Interest- 
ingly this model also shows remarkable phenomenological richness when a magnetic 
field is further introduced. This has recently been studied in detail in refs. [|50l 1771 . 
where upon a further analysis a holographic Berezinskii-Kosterlitz-Thouless transition 
has also been observed[l76|[5T]|. 

Undoubtedly both these models have remarkable rich features. Let us now go back 
and comment on the main theme of this thesis: the physics of chiral symmetry breaking. 
As we have noted earlier what we have called chiral symmetry in the models we dis- 
cussed are quite different. We also noted that the symmetry breaking pattern in the Type 
IIA model, i.e. the U{Nf)LxU(Nf)ji — )-U(A^/)diag5 is much closer to the chiral sym- 
metry breaking in QCD which we ultimately aim to understand better. This therefore 
poses a question: whether one can realize similar symmetry breaking pattern within 
some other brane-construction or supergravity background. To slightly rephrase this 
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question, we can focus on the fact that despite the beautiful realization of chiral symme- 
try breaking in the Witten-Sakai-Sugimoto model it has some serious drawbacks: e.g. 
the string coupling is not constant and diverges in the UV. Furthermore, the dual field 
theory is really a (4 + 1) -dimensional gauge theory compactified on a circle and not 
a true (3 + 1) -dimensional gauge theory. One is thus led to consider other possible 
supergravity models which might realize the chiral symmetry breaking pattern as the 
Witten-Sakai-Sugimoto model but do not suffer from its drawbacks. 

Recently in ref. HI 0011 a similar chiral symmetry breaking mechanism has 
been proposed by considering probe D7-anti-D7 branes in the Klebanov-Witten 
background[92j^ In ref. [|45l . this construction has been further generalized for the 
Klebanov-Strassler background[91]. Note that both these models are within the Type 
IIB supergravity framework. Evidently, these constructions provide us with the oppor- 
tunity to explore more. It would be very interesting to consider the effect of an external 
electric-magnetic field excited on the worldvolume of the probe brane-anti-brane sys- 
tem. Although the finite temperature Klebanov-Strassler background is not analytically 
known (for more details see e.g. ref. [|2|) and therefore hard to analyze, the zero temper- 
ature physics may still contain a rich physics. 

A crucial ingredient of our analysis is the large Nc and the Nc ^ Nf limit. In the real 
world we have Nc ~ Nf, thus a probe computation is not sufficient. Instead one needs 
to consider the effect of the backreaction of the probe flavour branes on the background. 
Early attempts to include backreaction in the Type IIB model has been discussed in 
e.g. ref. [|46ll . Some recent work has also been done to include the backreaction of the 
probe branes in the Witten-Sakai-Sugimoto model in ref. OTIl . However, this remains 
a much less controllable scenario. Similar attempts have also been made to include the 

^This background was originally found by Romans lll35l . 
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backreaction in e.g. refs. [[T9l |20l in the flavoured Klebanov-Witten and the Klebanov- 
Strassler model that we have not discussed here. 

Nonetheless, we hope we made a convincing case in favour of using the techniques 
of gauge-gravity duality to shed light on strong coupling dynamics. Let us now offer 
some comments from a broader perspective. String theory eventually emerged as a 
potential candidate theory for quantum gravity, however many aspects of this theory 
still remain poorly-understood. On the other hand, quantum field theories are much 
better understood. Thus, as proposed in ref. HIOSII , the gauge-gravity duality can itself 
serve as a definition of a theory of quantum gravity. The concrete examples of this 
gauge-gravity duality within string theory is therefore an exciting prospect. It is also 
natural to construct more such interesting examples. 

Continuing on this thought, it is tempting to ask whether we can realize the grav- 
ity duals of everyday systems. There are many condensed matter systems which are 
described by strongly coupled, scale invariant theories. These include high Tc super- 
conductors near their quantum critical point, fermions at unitarity etc. At the time of 
writing, constructing and analyzing gravity duals of various strongly coupled condensed 
matter phenomena is an extremely active and developing field of research. We refer to 
the refs. [[6911711 for a comprehensive review of the current progress. 

In conclusion we hope our approach of learning strong coupling dynamics by using 
the gauge-gravity duality, in one way or the other, will teach us interesting and useful 
lessons. 
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Appendices 



Appendix A: The fluctuation action of the probe in Type 
IIB model 

The relevant pieces of the action to second order in ol are: 

S = -Td7 I d^i^gab + Bab + 'iT^a'Fab + (27ra') /ly / F(2) A 5(2) A P [^(4) 

+ {27ra'ffijl [ F(2) A F(2) A P [C(4)] , (A.l) 

^ J Ms 

1 

C(4) = — —dt A dx^ A dx'^ A dx^ , 

^7(4) = (l — cos"^ 6) sin T/^ cos dip A d(j)2 A dcj)^ A dcpi , (A. 2) 

9s 

where P [C(4)] is the pull-back of the 4-form potential sourced by the stack of Nc D3- 



branes, P 



<^(4) 



is the pull-back of the 4-form magnetic dual to C(4), and F{2) is the 
Maxwell 2-form on the D7-brane world- volume. At this point, we resort to a different 
set of coordinates than we have been using. Instead of using the coordinates (p, L) 



introduced in equation (58), we return to the coordinates {z = 1/u ,6) because the 



analysis is simpler. We consider fluctuations of the form: 

e = 9o{z) + 2Wxr ) , 01 = 2W$(r) , (A.3) 

where the indices a,b = ... 7 run along the world-volume of the D7-brane. 9q{z) 
corresponds to the classical embedding from the classical equations of motion. Plugging 
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the ansatz in equations (A. 3 1 into the action and expanding to second order in {2na'), 
we get as second order terms in the lagrangian: 



C 



XF 



/-WZ 



wz 



(A.4) 



2% 
1 



--V^S'^'^R^ sin^ eoda<^db<^ 



-XF23 

1 ^ fr F mnop 

8 z'^R'^ ' 



-XF23f 



-^FoiHR'^smilJcos^ljdz (l - cos^6'o) = -^FqiHR'^ sin ijj cos ^/jd^K 



We have taken Eab = gfij + Bab to be the zeroth order contribution from the DBI action. 
In addition, we use that = 5"^ + J''^ where S""^ = S''" and J"^ = - J^". We use 
this notation for brevity. The indices m,n, o,p = A . . .9 run in the transverse directions 
to the D3-branes. 

To proceed, we impose Am = and dmAa = 0. We further assume that diX = 0, in 
which case it is possible to consistently set A2 = = A^. Finally we get the following 
equations to consider 



da 



1 + 



,4^4 



4^2 



dbX 



X 



(A.5) 








-^0 [y^S^^S^^d^Ai^ + d.KHdo^ - d, [y^S'' S^^d-Ai) , (A.6) 
da (y^S^-^R^ sin^ ^0^6$) - FmHR"^ sin ^ cos %l)d,K . (A.7) 
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Clearly we end up with one decoupled equation for % and a set of coupled equations for 
$ and A- fluctuations. 
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Appendix B: The fluctuation action of the probe in Type 
IIA model 



Here we provide the quadratic action for fluctuations of the probe brane. We consider 
the following fluctuation of the embedding function 



where ro{u) is the classical embedding of the probe brane and {^°} refer to its world- 
volume coordinates where a = 0, . . . , 8. 

Now we can determine the effective Lagrangian corresponding to small fluctuations 
by expanding the DBI action up to quadratic orders. We use {27ra') as a bookkeeping 
parameter for this expansion. The contribution at order (27rQ;')^ is summarised below 



r(u) = ro(u) + (2W)x(r) , 



(B.l) 




{daX){dbX) , 



(B.2) 



153 



where 



-1 



S""^ = diaK <! -G", G 



Gx 



Gl^ + H^' Gi+m'- 



X m 



411 5 



ab 



H 



[S^Sl - 52^1) , Quu = Guu + TQ{ufGrT , 



Glx + 

n^W =A:mg\^{ulRY^'^u\{ulRy'''^u',{ulRy'''^u\{u/R)-^'^u^^ , 
-det = ^/GuG^xQuu {Gl, + H^) (det 11^411) , (B.3) 



where S and A are the symmetric and anti-symmetric part respectively and 
denotes the diagonal metric for the 5*4. 

The term C^^p is the interaction term that couples the scalar and the vector modes. 
Clearly if Tq(u) = then the coupling vanishes. Therefore for the high temperature 
phase (i.e., when tq(u) = const.) the scalar and vector meson modes are always decou- 
pled. 

Now we determine the contribution coming from the Wess-Zumino term. The total 
contribution coming from the Wess-Zumino term at order {27ra')'^ is given by 



'5'wz 

^4 



'M9 

dC3 F2 = dA 



P[C5] AF2AF2--F4AAAF2AB2 
o 



with 



(B.4) 



In what follows we will focus on the special case of Tq(m) = in which case the first 



term in eqn. (B.4) does not contribute at the leading order. 



The equations of motion for the gauge fields then decouple from the scalar fluctua- 
tion and is given by (using R = 1) 



da 



(B.5) 
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where now a, a', b, b' E {M}''^ |J{^}' the probe brane worldvolume. It is now possi- 
ble to consistently set all the = where a E S^. This choice imposes a constraint 
of the following form 



This constraint can be trivially satisfied by looking at the gauge field fluctuations inde- 
pendent of the spherical directions (in other words focussing only on the SO (5) singlet 
states). 

It is straightforward to show that with the following ansatz for the gauge field fluc- 
tuations decouple into two sectors 



(B.6) 



At = At{u)e 



Ai = Ai{u)e 



— 



. 



(B.7) 



The choice A. 



■u 



imposes the following constraint 



e-*^-det(£;(o))5"" [S'\ujA'^) - S^\kA[)] + Ag;^H{uA^ + kAt) 



. (B.8) 



The equations of motion are given by 




(B.IO) 



(B.9) 



155 



and 



ft, 



y-det (£;(o))5225"M^ 3 



- e-*-^-det (E(o))522 (A;25" + cu^^") ^2,3 = . (B.U) 



In this case, the longitudinal modes (At, Ai components oscillating along {t, x^} plane) 
and the transverse modes (^42, A3 components oscillating along {t,x^} plane) clearly 
decouple. It is useful to note that A2, A^ transform as vectors and {At, Ai} transforms 
as a scalar under the SO (2) rotation group in the plane perpendicular to the magnetic 
field. 
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Appendix C: To obtain the effective Schrodinger equa- 
tion 

Here we outline the general variable changes to obtain the corresponding Schrodinger 
equation from the original equation for fluctuation modes. Without any loss of generality 
we can write the equation for any fluctuation mode as follows 



where, ai{u), 02 (m) and 03 (m) are known functions and prime denotes derivative with 
respect to u, fxiu) denotes the radial profile of the fluctuation mode and u is the cor- 
responding oscillatory frequency (in the finite temperature case this is the quasinormal 
mode). Let us rewrite fx{u) = a{u)g{u) with 



ai{u)f'^{u) + a2{u)f'^{u) + u:'^az{u)fx{u) = , 



(C.l) 



1 a2{u) 1 / ai{u) 



) 



9, 



( 




(C.2) 



a{u) 



2 ai{u) 2 \a3(M) 



u 



The Schrodinger equation is then obtained to be 




(C.3) 



where 




(C.4) 



(C.5) 
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A more conventional form of the Schrodinger equation can be obtained after changing 
variables to "tortoise" coordinate du = B{u)du\ the horizon is then located at £t — )■ oo. 
The main goal of this exercise is to obtain the potential Vs{u) and extract qualitative 
features of the meson spectrum. 

To determine the quasinormal modes we consider only in-falling modes at the hori- 
zon; such a solution can be written as g{u) = exi)(—iujt)ilj{u). Now multiplying equa- 



tion in (C.3 ) by ip*{u), integrating by parts from the horizon to the boundary and using 



the equation of motion we get 

Note that the left hand side is positive definite as long as Vs is positive. This implies 
that lm[uj] < and therefore the corresponding mode is stable. On the other hand, 
a sufficiently strong negative support of Vs would indicate the presence of tachyonic 
modes in the fluctuation spectrum and hence an instability. 

Let us now represent the quasinormal frequency as w = + iui. It can be noted 
from the effective Schrodinger equation that there is a Z2 symmetry in (namely, 
if g{u) is a solution with a frequency tu then g*{u) is also a solution with frequency 

Here we have explicitly assumed that the fluctuation does not have any momentum 
mode. The fluctuation equation for the meson having a momentum can also be recast 
in the form of a Schrodinger equation. The changes of variables are exactly similar as 



already mentioned in eqn. (C.3 1, but the effective potential receives a positive contribu- 



tion coming from the momentum 



Vsiu,k) = Vsiu,0) + e{^-^] , (C.7) 
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where the left hand side of the equation ( |C.1[ ) is now accompanied by another term of 
the form — \a4{u)\k'^f^, the relative negative sign is due to the fact that /c-corresponds to 



a spatial oscillation. 
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Appendix D: Lowest hydrodynamic mode in vanishing 
magnetic field 

In holographic contexts the study of the lowest-lying (hydrodynamic) quasinormal 
modes of various fluctuation fields in a given background has become a vast literature 
by now. Initial work along this direction was done in e.g. ref. l|149[ll30lll31[|T46ll and 
later reviewed in e.g. ref. II150L These studies were done for the M = A SYM the- 
ory and hence within the Type IIB supergravity. Similar analysis was also done for the 
Witten-Sakai-Sugimoto model in e.g. ref. ifTSlfTTll . 

Here we analyze the gauge field fluctuations when there is no magnetic field. We 
provide an analytical derivation for the dispersion relation of the lowest hydrodynamic 
mode of the longitudinal oscillation, which has been obtained by other methods in refs. 
[|115[ l52l . Similar analytical study of a U(l) fluctuation field in AdSp+2 background 
was done in refs. [|87lll50i 

The longitudinal mode is defined to be £{u) = tuAi + kAf. Using this definition, the 



constraint equation in (B.8) and the equations of motion in (B.9) (after setting B = 0) 



we obtain the equation of motion for the longitudinal mode to be 

d'^S 5 (1 - - P (1 - x-^)) + Gu^x'^ dS 

dx"^ 2x {1 — x^^) {uj"^ — kP' {1 — x^'^)) dx 

+ ^ ^ = , (D.l) 

where x = u/Ut- 

The near-horizon hmit u — )■ f/^ is now achieved by taking the x — j- 1 limit. In this 
hmit, the equation of motion takes the form 
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The general solution of this equation is given by a linear combination of the incoming 
and the outgoing modes 



£{x) = Cicos(£l;log(x - 1)) + C2 sin(a; log(a; - 1)) , oj = \l = -(D-S) 



The incoming boundary condition singles out the solution with Ci = —C2. Near the 
horizon £(x) is therefore obtained to be 



S{x) = C^Exp {—iulog{x — 1)) 



(D.4) 



where Cl is an yet undetermined constant. 

Let us also define k = k/{ATxT). Now, for small enough u and k (such that a) ^ 1 



and k ^ 1), we can ignore the last term in equation (D.l I. It turns out that in this limit 
the equation of motion for the longitudinal mode is also exactly solvable. 



1 / 2k'^ 



uP' log 



X 



3/2 + 1 



3,3/2 _ 1 



(D.5) 



where Mi and M2 are constants of integration. Near the boundary, i.e., x — 00, this 
reduces to 



(D.6) 



Therefore normalizability oiE{x) forces us to impose M2 = 0. The quasi-normal modes 
are therefore obtained to be the solution of this constraint. 
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On the other hand, in the vicinity of the horizon the solution in ( |D.5[ ) takes the 
following form 



S{x) = ^ (^-2MiP + 3M2 + Miu^ log (^^^ ^ - ^Miw^ log(x - 1) + 

-Ml ( 6P + c:;^) (x - 1) - ^7^1 ( + c:;^) (a; - 1)2 + . . . 
6 V / 144 V / 



(D.7) 



Now, for sufficiently small values of u, we can expand the solution obtained in (D.4 1 
to get 



S{x) =Cl- iCLu\og{x - 1) - hoLU^ (log(x - 1))' 



(D.8) 



Comparing equation (D.7 1 and (D.8 1 we get 



^ f -2MiP + 3M2 + Miw^ log (/^ 



Cr 



(D.9) 



Setting M2 = 0, the solution of this equation is given by 



(D.IO) 



Restating this result in dimensionful parameters, we get the lowest hydrodynamic quasi- 
normal frequency 



u) = -iDnk'^ , Dr 



1 

2^ 



(D.ll) 
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where Dj^ is the R-charge diffusion constant. This value for the diffusion constant 
was obtained numerically by studying the spectral functions of holographic flavours 
in ref. [I115II and was also confirmed in ref. [[52ll by numerically obtaining the lowest 
quasinormal mode of the longitudinal oscillation. Not surprisingly we reproduce the 
(3 + 1) -dimensional result of ref. HI SOU since the flavour branes describe an effective 
(3 + 1) -dimensional field theory. 

On the other hand, it is easy to show that the transverse mode does not have any 
solution compatible with the hydrodynamic limit. See, e.g. ref. [1971 . 
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Appendix E: Global anomalies and dispersions from first 



Here we will briefly review the dispersion relations resulting from ideal hydrodynamics 
in the presence of a quantum anomaly [I147II . Our goal here will be to demonstrate the 
existence of the "anomalous sound" mode (observed in section 6.4) within the hydro- 
dynamic equations. Let us consider a (3 + 1) -dimensional relativistic fluid with one 
(classically) conserved U(l) which results in an U(l)'^ anomaly (after quantization). 
One example of such will be the QCD triangle anomaly coming from the chiral U(l) 
which is a true symmetry of the system at the classical level. 

We also turn on a slowly varying non-dynamical background gauge field which 
couples to the current J^. We assume that ~ 0{k'^), where k is the momentum, and 
therefore F^^ ~ 0{k). We will work within the framework of first order hydrodynamics 
and therefore keep terms of order 0{k) in the constitutive equations for T'^" and J^. We 
will closely follow the set-up and conventions of ref. [I147II . 

In the presence of external fields (such as electric -magnetic fields) the hydrodynamic 
equations take the following form 




A , 



(E.l) 



where C is the anomaly coefficient and we have also defined 



1 



(E.2) 



U 5 



2 



where is the fluid velocity vector normalized as u^u^ 



-1. 
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Note that B'^' is identically zero for lower dimensional theories. The non-zero right 



hand sides of the hydrodynamic equations in (E.l I reflect the work done by the external 
field on the system and on the anomaly [|147L 

The constitutive equations for T^" and are given by 

JA* = puf' + u^', (E.3) 



where 

A^"" = 7]^"" + uV , r^^^ = diag(-l, 1,1,1) , (E.4) 



and e, P and p denote the energy density, the pressure and the charge density respec- 
tively. The terms U^'^ and are of 0{k) and capture the dissipative effects. These 
dissipative terms are obtained from the requirement of the existence of an entropy cur- 
rent with non-negative derivative. The result of such an exercise is given by lll47ll 

U>"' = -ria"" - CA^"" {d^u") , 

and 

= aQA^'[-d,p + F,^u'' + ^d,T)+Cv.w^ + CBB^, 

= ^e^^'^'^u^d.u^ , (E.6) 

where d is the spacetime dimension (which in this case is 4), r] is the shear viscosity, ( is 
the bulk viscosity and (^u), Cb are the new transport coefficients sourced by the anomaly. 
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We will consider linear fluctuations around the equilibrium configuration. The fluc- 
tuation variables will be the velocity vector 5u^^ as well as the temperature 5T and the 
chemical potential Sji. The normalization condition of the fluid velocity vector requires 
Su^u^ = 0, i.e., the fluctuations be orthogonal to the equilibrium configuration. 

We choose the following equilibrium configuration 



Therefore our fluctuation variables are: Su^, Suy, Suz, ST and S/j, and all fluctuations are 
of the plane-wave form exp{—iujt + ikx). The equation for the fluctuations modes are 

given by 



F2Z — H — — F32 , F12 — — F31 , 



(1,0,0,0) , T = const. , /X = const. 



(E.7) 



A , 



d^5J^ = C5 {E^B^) 



(E.8) 



where 



(E.9) 



From these we can obtain the equations for the fluctuations to be given by 




= 0, 



0, 



0, 



(E.IO) 



(E.ll) 



(E.12) 



166 



and 



w (e + P) + ik'^T] + iagH^ + -(^kH 



6uy — ipHdUz 
6uz + ipH6uy 



, (E.13) 
. (E.14) 



Clearly the fluctuation sectors {Sp, 5T, Su^} and {5uy, Su^} are decoupled. 

Let us first we consider the case when H = 0, p = and p = 0. The equations of 
motion for the fluctuation modes take the following form 



u5e — k {e + P) Su^ 



4 



(e + P) 6u^ - k6P + ik^ [ C + ] = , 



and 



0, 



(E.15) 



LoSp — ik'^aqSp = 
[uj (e + P) + ik'^T]] 6uy^z = 



(E.16) 



The dispersion relation resulting from the equations in (E. 15 ) is given by 



u = ±Cc:k — ik"^ ^ 



2(6 + P) 



+ . 



dP 



(E.17) 



where Cg is the speed of sound (more generally the speed of the propagating degree of 
freedom). On the other hand, the dispersion relation resulting from the last equation in 
( |E.16| ) is given by 



CO = -ik^^^ + ... 
e + P 



(E.18) 
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whereas the first equation in (|E.16[) describes a diffusive mode for 5/i given by 



(9/i 
dp 



+ ... 



(E.19) 



Therefore equation ( E. 1 5 ) is known as the sound mode equation and ( E. 1 6 1 is known as 



the shear mode equation. 

Let us now focus on a bit more general case when H ^ 0, but p = = p. In 
this case, the diffusive mode corresponding to 5p couples to the sound mode equation 
and thus it does not make sense to address them separately. The other diffusive mode, 
corresponding to Suy or Suz, however does remain decoupled. Let us write down the 
equations below: 

1. The coupled mode: 



coSe — q{e + P) 6ua 



4 



oo{e + P) 6ux - q6P + iq' \^ + -r^J du^ 
uSp + u(bH6ux — iq^aqdp 



0, 
0, 
. 



(E.20) 



2. The shear mode: 



(e + P) + iq^T] + laqH^ + -Cn,qH 



u{e + P)+ iq^T] + mqH^ + -Cu,qH 



Su,. 



0, 
. 



(E.21) 



From the shear channel, we readily obtain the following dispersion relation 



-ik 



2 V 



e + P e + P 2 (e + P) 



(E.22) 
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Assuming that H is small in the hydrodynamic limit (specifically if we focus on the 
regime when H ~ k), we get a new propagating mode whose velocity is now given by 

^anomalous jj 

Therefore the presence of an external magnetic field can source a propagating mode 
in the otherwise diffusive channel. Note that this effect is completely universal, i.e., 
present for both conformal and non-conformal systems obeying arbitrary equation of 
state (since we did not use the equation of state anywhere). The magnetic field also 
modifies the diffusive mode. 

On the other hand, the coupled modes have a simple behaviour. The corresponding 
characteristic equation factorizes in the following diffusive and sound channels 

^ = ^k^^Q_^ ^ = ±Csk-ik^^y^i^. (E.24) 
e + P' ' 2(e + P) ^ ^ 

In other words, the magnetic field does not affect the already existing shear and sound 
channels. 

Let us now consider the case when all the external parameters are present i.e. H ^ 
and p 7^ 0. As alluded to earlier we will consider only {6T, Sji, 5ux} and {5uy, 5uz} as 
the independent fluctuation modes, which incidentally do remain decoupled. 

The shear channel dispersion relation now takes the form 

u} = —ik — — iH — ^ ± —. zzrkH ± — . (E.25) 

e + P e + P 2 (e + P) e + P ^ ^ 
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This has the same qualitative features as before, namely the propagating mode still 
shows up. The only difference now is the new term proportional to the charge den- 
sity and Unear in H. Note that this dispersion relation should also be rather general, 
since we have not assumed any equation of state specific to any system. 

The other set of coupled mode can also be solved to obtain the corresponding dis- 
persion relation, but this is cumbersome for the most general case. We need to know the 
equation of state to obtain the specific dispersion relation of a particular system. 
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Appendix F: Spinning strings and background field in 
Rindler space, a model calculation 

Here we would like to solve a toy problem to study the effect of the magnetic field on 
meson dissociation. We implicitly assume the framework where we can have two equal 
charges at the two ends of the string and therefore analyze the symmetric configuration 
only. Rindler spacetime provides an useful arena where many qualitative (if not quan- 
titative) features of an event horizon can be realized in a simple set up. Such study has 
been previously carried out in ref. [I123II without any background field. The main result 
of this exercise is that the high spin meson dissociates once a critical value of the accel- 
eration is reached. This acceleration is in turn determined by the angular momentum of 
the meson. 

There are two ways to physically interpret the results of such pursuits (e.g., see 
ref. [I123II and ref. [1211 ). We consider accelerating the rotating string in one of the space- 
like directions and therefore a Rindler horizon forms in the four dimensional spacetime. 
Alternatively, we can consider accelerating the string in the holographic direction and 
therefore a bulk spacetime horizon forms which sets the temperature of the dual gauge 
theory. To begin with we take the former point of view. 

The metric for the Rindler space can be written as 

rfs^ = -i'^K^d^f + di^ + dp^ + p^dct? , (F.l) 

where {p, 0}-plane represents the plane where we would rotate the string. To add a 
background magnetic field we consider similar potential given before. Now the 
ansatz for the string in terms of its worldvolume coordinates {r, cr} is given by 

^ = ^ , i = , P = p{(y) , (p = 0JT . (F.2) 
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This ansatz implies that the direction of acceleration of the string is orthogonal to the 
direction of its rotation. We can readily see that there are two conserved charges asso- 
ciated with the string, namely the angular momentum J, which is associated with the 
rotation along 0-direction and the boost charge, which is associated with the translation 
in ?7-direction. 

The Nambu-Goto action for the string now will be accompanied by a boundary term 
exactly similar to the term ASb given in eqn. ( |151| ). Therefore the full action is given 
by 

S=^ [ rfrrfav/(e«^2_^2p2)(^/2 + p/2) + /s,Sb . (F.3) 
2na' J 

Varying this action we can obtain the equation of motion for the string profile and the 
boundary conditions we should impose. In absence of any external field this exercise 
has been carried out in details in ref. [I123II . which we do not repeat here. The presence 
of the magnetic field enforces us to impose the following boundary condition 



^pl±£/2 



dL 



= -Hpu , (F.4) 

±e/2 



where £ is the distance between the two end points of the string. The other boundary 
term corresponding to the variation of ^ is satisfied imposing the Dirichlet boundary 
condition 5^ = 0. This corresponds to a situation where the end points move with a 
constant acceleration a = ^(±£/2)^^. 
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Now choosing a gauge p = a we can solvqj the equation of motion to obtain the 
profile for the string given by 



^(p) = —cosh — arcsin f — 



with C2 = — (l + i/^^ 



4 

1 f „ox / 1 \ „ 1 



As the notation suggests the parameter J represents the angular momentum of the 
meson. 

As we have seen before, the presence of the background magnetic field does not 
change the equations of motion for the Nambu-Goto string. Therefore the functional 



form of the string profile given in equation (F.5 1 is entirely determined by the Rindler 
geometry HI 231 . We notice that the presence of the magnetic field does induce an effec- 
tive length £eff = ^^(1 + by changing the boundary condition of the string end- 
points. 

A non-zero value of H brings about two key changes in the evaluation of the angular 
momentum. The profile of the string modifies in order to satisfy the boundary condi- 



tion in eqn. (F.4) and therefore the effective length 4ff arises in the formula in (F.5). 



There is also an additional explicit contribution AJ^ of the form shown in eqn. ( 153), 
which finally gives a functional factor of tan^^(l/iJ) in the expression for the angular 
momentum. 

It can be shown that the explicit contribution of the of the magnetic field to the 
angular momentum, denoted by A J^, is actually cancelled by a term coming from inte- 
grating over the string profile. The net resulting angular momentum therefore can be 



'The analytical solution of the equations of motion obtained from the Nambu-Goto action remains the 
same as in ref. II123I . We refer to ref. 11231 for further details about solving the equations of motion. We 
use the same analytical solution and impose the boundary condition given in ( |F.4| l to obtain the results in 
eqn. (|E3i. 
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lower than the angular momentum of the string in absence of the field. This is pictori- 



aUy summarised in figure 3 1 




Figure 3 1 : The total angular momentum of the string in units of (27ra') in the background field. 
The horizontal dashed line represents the value of angular momentum of the string when there is 
no background magnetic field. We observe that the S-field initially reduces and then increases 
the angular momentum creating a local minima. However, we will observe that the stability of 
the meson is always enhanced by the S-field. 



To find the maximum acceleration we identify oT^ = and look for minima of 

the right hand side of the equation as a function of the parameter C. For now we content 
ourselves with only positive values of H, i.e. when the string angular momentum and the 
background field are parallel. It is a straightforward exercise to show that such minima 
corresponds to values of C which satisfy the following equation (obtained from taking 
the first derivative of ^ and setting it equal to zero at p = ±£/2) 



X = mtanh 



m 



X 



C 



X 



where 



m = arcsm 



arctan 



H 



-1/2 



(F.6) 



Clearly the parameter H generates a set of such values of C which satisfy the above 
equation and therefore promotes the maximum value of acceleration to a function of the 
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magnetic field. Differentiating the eqn. in the first line of ( |F.6[ ) and using identities for 
the trigonometric hyperbolic functions we obtain that the family of roots of the equation 
can be simply given by the relation x = am where a = 0.834 is a constant. 

Using these roots the maximum acceleration can be obtained as a function of the 
background magnetic field which is given by 



Omax oc \ / aictan f^j ^ ^ . (F.7) 

V /arcsin(^^^ 

Alternatively we could imagine the Rindler space to be the near-horizon approxima- 
tion of a blackhole background and the acceleration is along the holographic coordinate 
(we make the identification that ^ = u, where u is the radial coordinate in the holo- 



graphic set up). The parameter k (in eqn. (F.l)) in this case represents the surface 
gravity and sets the temperature of the dual gauge theory to be T = k/{2tt). 

Now we imagine the string end points to sit on the flavour brane at some value of 
Co = Uq, which fixes the constituent quark mass; and we consider a similar ansatz for 



the string profile as presented in eqn. (F.2). The interpretation of the set up is different 



from the earlier one; however this brings about no change as far as the mathematics is 



concerned. We therefore again obtain the string profile as given by eqn. (F.5 ). 



However in this case the physical parameter that captures the dissociation of the 



meson is the size i (which sets the angular momentum J via the relation in eqn. (F.5 1) 
of it. The length is obtained to be 

2(:7arccosh(K[/o/C) /^—eh\]'^ (pg) 



arcsm 



We would hope to see that for a given value of Uq, there exists a maximum admissible 
size of the spinning meson beyond which it dissociates. 
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We can treat as before the parameter C to be independent with respect to which 



we will consider maximizing the length function £. From eqn. (F.8) it is clear that the 
external magnetic field does not play any role to determine the value of C corresponding 
to the maximum of t, however it determines the the function imaxiH) and therefore also 
^max(i^) to be given by 



arcsm 



■/max C)C arctan ^— ^ arcsin + H"^ 



(F.9) 
(F.IO) 



where the constants of proportionality depends on the background temperature T and 
the parameter Uq. 

For the background considered in the main text the Rindler spacetime does emerge 
when we zoom in the region close to the horizon li55l IllOH . The toy model here can 
therefore be identified with a corresponding study of high spin mesons in the overheated 
phase. So if we want to make a direct connection between the Rindler background 
dynamics to the specific system studied in the text, we should remember that such a 
background appears precisely in the over-heated phase of the system, where any small 
fluctuation is likely to destroy the meson and drive the over-heated phase to the melted 
phase. Our emphasis here however is not to make this connection, we view this exercise 
as a more tractable analysis of some of the physics we studied in the more thorough 
holographic calculations. 
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